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SUMMARY 


This  Final  Report  under  Contract  No.  FI9628-78-C-0200  summarizes  the  results  of  a  data  reduction  and  analysis 
program  addressing  optical  atmospheric  properties  in  the  western  European  environment.  The  data  were  collected  in  a 
cooperative,  but  independent  effort  associated  with  the  NATO  Research  Study  Group  8  of  Panel  IV,  AC243  and  their  Pro¬ 
gram  OPAQUE. 

The  Aircraft  Data  Reduction  effort,  task  one  under  this  two  part  research  contract  has  concentrated  on  the  measure¬ 
ments  of  visible  spectrum  atmospheric  optical  properties,  primarily  scattering  coefficient  profiles  plus  sky  and  terrain  radi¬ 
ances,  that  were  collected  during  the  Winter  and  Summer  seasons  of  1978.  During  this  interval  there  were  approximately 
45  data  flights  conducted  throughout  western  and  northern  Europe.  Over  b70  aerosol  si/e  distributions  covering  particle 
radii  from  0.2  to  5.9  micrometers  were  measured  contemporaneously  with  the  measurements  of  sky  and  terrain  radiances 
during  these  flights.  Both  the  optical  and  aerosol  data  are  reviewed  in  this  report. 

Several  general  characteristics  of  the  volume  scattering  coefficient  profiles  are  illustrated  showing  the  Trapani.  Sicily 
profiles  as  clearly  different  from  those  measured  in  northern  Europe.  The  Trapani  profiles  show  significantly  less  low  alti¬ 
tude  structure,  often  with  a  complete  lack  of  any  abrupt  haze  top.  Summer  hazes  in  this  Mediterranean  area  remain  heavy 
throughout  the  entire  6  kilometer  altitude  regime. 

Very  low  altitude  scattering  coefficient  profiles  measured  between  ground  level  and  approximately  700m  AGL. 
representing  all  four  temporal  seasons  and  all  northern  and  southern  sites,  indicate  that  in  lorty-seven  out  of  fifty-four 
(87%)  instances,  the  profile  is  essentially  constant  in  value  within  this  very  low  altitude  regime  during  the  middle  to  late 
afternoon  measurement  intervals. 

A  high  degree  of  radiance  symmetry  and  uniformity  is  illustrated  for  sky  and  terrain  radiances  measured  under  full 
overcast  conditions.  Conversely,  under  conditions  of  highly  variable  cloud  cover,  the  measured  sky  and  terrain  radiances 
tend  to  fluctuate  both  above  and  below  the  radiance  distribution  of  a  clear  sky  by  a  factor  ol  two  or  more  in  some 
instances,  but  in  general  remain  within  about  ±50  percent  of  the  clear  sky  radiance  for  corresponding  paths  of  sight 

Several  experimental  uncertainties  associated  with  the  measurement  of  high  altitude  scattering  coefficients,  and  near 
sun  sky  radiances  have  been  addressed,  and  corrective  procedures  have  been  devised  to  reduce  the  errors  existing  in  the 
original  measurements.  The  early  evidences  of  stray  light  errors  in  the  nephelometer  data  resulted  in  several  mechanical 
modifications  to  the  instrument,  but  continuing  analysis  indicated  an  additional  residual  error  in  the  measurements.  These 
analyses  indicated  that  the  error  was  most  probably  the  result  of  incomplete  aerosol  exchange  within  the  nephelometer 
sampling  volume  as  the  measurement  altitude  increased.  It  is  clear  that  an  incomplete  flushing  of  only  a  fraction  of  one 
percent  of  the  low  altitude  aerosol  concentration  can  result  in  a  substantial  erroneous  signal  when  the  true  high  altitude 
concentration  has  dropped  over  an  order  of  magnitude.  The  evidences  from  these  analyses  are  consistent  in  their  predic¬ 
tions  of  the  true  optical  scattering  ratios  0(z)  being  1.3  to  2.6  at  6  kilometers  altitude,  roughly  half  the  measured  values. 
Sky  radiance  measurements  made  by  the  automatic  2n  scanner  system  were  found  to  contain  stray  light  errors  which  con¬ 
sistently  increased  markedly  in  magnitude  both  as  the  instruments  field  of  view  approached  within  25°  of  the  solar  disc, 
and  as  the  atmospheric  aerosol  loading  decreased  An  approach  to  the  correction  of  these  near-sun  radiance  values  has 
been  presented  earlier  in  Johnson  (1981). 

The  aerosol  size  distributions  were  analyzed  using  the  dV/d  log  r  volumetric  format  resulting  in  characterization  by 
three  additive  log-normal  distributions,  one  for  each  of  the  three  volume  modes  observed  in  the  data 

The  accumulation  mode  (r  =  0.25  to  0.47^m.  <r  =  0.1 1)  is  found  almost  always  confined  to  the  mixing  layer,  and 
appears  more  frequently  in  the  data  collected  at  the  northern  stations  than  those  in  the  south.  The  middle  mode  <7  =  0.45 
to  \ u  =  0.3)  occurs  only  in  the  mixing  layer,  appearing  most  often  in  Winter  and  when  the  relative  humidity  is 
greater  than  70  percent.  The  course  particle  mode  (7  =  5fA.ni,  <r  ■»  0.2)  can  exist  above  or  within  the  mixing  layer,  and 
shows  no  preference  for  season,  but  appears  most  often  in  the  data  from  southern  stations  and  when  relative  humidity  is 
less  than  70  percent. 

The  data  and  comparisons  with  aerosol  models  imply  that  the  volume  modes  are  associated  with  specific  aerosol 
sources,  and  therefore,  with  the  history  of  the  airmass  in  which  the  aerosol  sample  was  measured. 
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A  REVIEW  OF  MEASURED  ATMOSPHERIC  OPTICAL  PROPERTIES 
AND  THEIR  CONTEMPORARY  AEROSOL  SIZE  DISTRIBUTIONS 


Richard  W  Johnson 
Bruce  W.  Filch 


1.  INTRODUCTION 

This  Final  Report  has  been  prepared  under  AFGL 
Coniraci  No.  F19628-78-C-0200  and  addresses  the  first  of 
two  tasks  enumerated  under  the  terms  of  the  contract,  the 
reduction  and  publication  of  previously  unprocessed  raw 
data  collected  by  an  instrumented  C-130  aircraft  during 
the  life  of  a  flight  program  which  terminated  in  Sep¬ 
tember,  1978.  The  flight  program,  related  to  the  measure¬ 
ment  of  atmospheric  optical  properties,  was  conducted  by 
the  Visibility  Laboratory  of  the  University  of  California, 
San  Diego  under  the  sponsorship  of,  and  in  cooperation 
with  the  Air  Force  Geophysics  Laboratory,  Hanscom 
AFB,  Massachusetts. 

The  instrumentation  and  experimental  procedures 
used  in  the  collection  of  the  raw  data  processed  during  this 
contract  interval  were  reviewed  and  summarized  in  Dunt- 
ley  ei  at.  (1978),  and  thus  will  receive  only  minimal  addi¬ 
tional  comment  at  this  time.  The  airborne  instrument  sys¬ 
tem  is  illustrated  in  Fig.  1-1. 

The  optical  measurements,  which  consisted  pri¬ 
marily  of  profiles  of  atmospheric  volume  scattering 
coefficient  and  47r  sky  and  terrain  radiances,  were  made  in 
several  spectral  bands  throughout  the  visible  spectrum  as 


illustrated  by  the  curves  in  Fig.  1-2.  The  experimental 
flights,  which  were  associated  with  the  NATO  program 
OPAQUE,  were  made  along  flight  tracks  throughout 
several  European  regions  as  illustrated  in  Duntley  ei  at. 
(1978)  and  as  summarized  in  Fig.  1-3  and  Tables  1.1  and 
1.2  which  have  been  abstracted  from  that  earlier  report. 

In  addition  to  the  radiometric  measurements  made 
by  the  various  radiometer  systems  illustrated  in  Fig.  1-1, 
there  were  airborne  measurements  of  particle  size  distribu¬ 
tions  at  each  of  several  altitudes  during  each  data  flight 
These  distributions  were  measured  over  the  range  0.2  to 
5.9  micrometers  (inn)  in  particle  radius  using  a  Royco 
model  220  particle  counter.  An  initial  discussion  of  these 
aerosol  measurements  was  presented  by  Cress  (1980),  and 
by  Fitch  and  Cress  (1981). 

In  the  following  sections  of  this  Final  Report  on  the 
Aircraft  Data  Reduction  task,  we  will  summarize  and 
review  the  dua  processed  and  published  subsequent  to  the 
status  reported  in  Duntley  et  al  (1978).  Section  II  will 
present  a  discussion  of  the  radiometric  data,  and  Section 
HI  will  contain  an  analysis  of  additional  measurements  of 
aerosol  size  distributions. 


Hr.  1-1.  C- 1. to  airborne  instrument  sysiem 
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Table  1.1.  Project  OPAQUE  Flight  Track  &  Ground  Sites 
Location  &  Ground  Elevations. 


Geographical  Reference 

Center  of  Right  Track 

Ground  Site 

Approx 

Sue  Identification  Code 

Latitude 

Longitude 

Latitude 

Longitude 

along  track 

Ground  Sues 

Catania.  Sicily 

(CT) 

. 

37*24*N 

14'55'E 

Trapani,  Stctly 

(TR) 

37*55*N 

1 2°29'E 

Bruz, France 

(BR) 

4g°01’N 

1*451* 

Btrkhof,  Germany 

(BKJ 

. 

48°I3*N 

nrE 

Yeovillon.  England 

(YO) 

- 

snorts 

2*371* 

Soesterberg.  Netherlands 

(SO) 

52°08'N 

5®17'E 

Meppen.  Germany 

(MP) 

52°52'N 

7®23’E 

Right  Tracks 

Sigonella.  Sicily 

ISO) 

37°24‘N 

IS’20'E 

- 

. 

Sea  Level 

Trapani.  Sicily 

(TR) 

37’33'N 

12°30’E 

Sea  Level 

Bruz,  France 

(BRI 

48’01'N 

l°41*W 

50  m 

Btrkhof.  Germany 

<BK> 

48“15'N 

9*0$ ’E 

762  m 

Yeovil.  England 

•  YO) 

50°56'N 

2°27’W 

60  m 

Soesterberg.  Netherlands 

(SO) 

5P56*N 

5°35'E 

6  m 

Mtldenhall.  England 

(Ml.) 

S2°24'N 

1‘41’E 

Sea  Level 

Ahlhorn  Germany 

(AL) 

52°53'N 

7°51'F. 

18  m 

Meppen.  tier  many 

(MP) 

53*00'N 

,037'E 

18  m 

Rodbv  Denmark 

IRB) 

54'41'N 

trosi: 

Sea  Level 

Table  1.2.  Geographical  Distribution  of  Protect  OPAQUE  Data  Flights 

Attempted  (and  reported  I  Data  Sequences 


Flight  Track  Idem 

OPAQl.'E  1 

OPAQl/i:  II  , 

OPAQUE  III 

OPAQUE  IV 

OP  AQUE  V 

Spring  '’6 

Fall  ‘76 

Summer  ‘77 

Winter  ‘78 

Summer  ‘78 

Sigonella.  Sicily 

n 

0 

0 

2  (2) 

0 

Trapani.  Sicily 

0 

1) 

0 

3  (3) 

4  (4) 

Bru/.  France 

0 

5  (<) 

3  <J> 

0 

0 

Birkhol.  Germans 

0 

0 

0 

8  (7) 

5  (41 

Yeovil.  England 

4  (3» 

0 

0 

5  (5) 

3  (3) 

Soesterberg.  Netherlands 

t  M) 

() 

1  (0) 

2  <2) 

3  (3) 

Mtldenhall,  England 

0 

0 

0 

2  (2) 

0 

Ahlhorn.  Germany 

11 

n 

2  12) 

D 

n 

Meppen.  Germany 

<  12) 

4  M) 

3  (3) 

3  (3) 

3  (3) 

Rodby  Denmark 

2  <2) 

4  14) 

4  (4) 

2  (2) 

2  (2) 

TO  l  AL 

13  (8) 

13  ( 1 2 ) 

13  02) 

27  (261 

20(19) 

Technical  Report  No 

AFGl  FR 

AFGl  TR 

AFGl.  TR- 

AFGLTR 

AFGLTR- 
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2.  AIRBORNE  RADIOMETRIC  MEASUREMENTS 

During  the  three  year  interval  1975-1978,  there 
were  five  major  overseas  deployments  of  the  airborne 
instrument  system  illustrated  in  Fig  1-1  as  an  adjunct  to 
the  program  OPAQUE.  These  five  trips  were  planned  to 
provide  daytime  atmospheric  measurements  as  a  function 
of  altitude  in  the  immediate  vicinity  of  several  different 
OPAQUE  data  stations  The  plan  was  developed  to 
ensure  a  reasonably  diverse  geographical  sampling  of  verti¬ 
cal  profile  data  during  each  of  the  four  seasons  of  the 
year.  Spring,  Summer.  Fall,  and  Winter  Data  for  the  first 
three  of  these  deployments  was  summarized  and  reviewed 
in  Duntley  era/.  (1978)  In  a  similar  manner,  the  data 
from  the  last  two  deployments.  Winter  and  Summer  1978 
is  summarized  in  the  following  paragraphs 

Table  I  I  lists  the  general  location  of  each  of  the 


different  flight  tracks  used  during  the  OPAQUE  deploy¬ 
ments.  and  also  the  locations  where  the  air  transportable 
ground  unit  was  established.  The  general  locations  are 
illustrated  in  Fig.  1-3.  where  the  short,  solid  lines  indicate 
the  approximate  orientation  of  each  flight  track  For 
further  more  detailed  site  location,  the  reader  is  referred 
to  the  site  maps  in  Section  6  of  Duntley  er  al.  (1978). 


2.1  Description  of  Winter  and  Summer  Deployments 

As  a  convenience  to  the  reader,  the  deployment 
descriptions  for  the  Winter  1978  and  Summer  1978  field 
trips  have  been  abstracted  from  Duntley  el  al.  (1978). 
and  reproduced  below 
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OPAQUE  IV,  WINTER  1978, 

FLIGHT  NOS,  C-430  THROUGH  C-456 
AFGL-TR-79-0159 

The  Winter  deployment  began  with  the  Visibility 
Laboratory  team  departing  San  Diego  on  Wednesday,  25 
January  1978.  The  team  returned  on  Monday,  10  April 
1978.  During  this  72  day  deployment,  the  longest  of 
those  up  to  this  time,  the  joint  Air  Force  and  Visibility 
Laboratory  team  operated  the  aircraft  out  of  four  different 
bases  in  three  different  countries,  not  counting  the  initial 
staging  into  Rhein-Main  Air  Base.  Two  new  bases,  never 
before  visited  during  the  OPAQUE  related  missions,  were 
added  to  the  itinerary,  one  in  Sicily  and  one  in  southern 
Germany.  The  chronological  listing  of  the  operational 
bases  utilized  during  the  OPAQUE  IV  deployment  is 
shown  below: 


Order 

Visited  Staging  Base  Idem  Geographical  Location 

1  Rhein- Mam  Air  Base  near  Frankfurt.  Germans 

2  Sigonella  Naval  Air  Facility  near  Catania.  Sicily 

3  Memmingen  Air  Base  S-W  of  Munich.  Germany 

4  I  Mtldenhall  Air  Base  N-W  of  Cambridge.  England 

5  |  Wunstorf  Air  Base  near  Hannover.  Germany 


The  only  major  instrumentation  modification  imple¬ 
mented  for  the  OPAQUE  IV  deployment  was  with  respect 
to  the  MEMORY  cutoff  system  This  modification  was 
flown  in  prototype  during  OPAQUE  III.  The  improved 
performance  of  the  dual  irradiomeier  system  upon  which 
this  modification  had  been  retrofitted  was  so  impressive 
that  the  decision  to  extend  the  modification  to  all  systems 
was  made  immediately  upon  return  to  the  Laboratory  in 
September,  1977.  The  modification  was  accomplished  on 
all  of  the  airborne  and  ground  based  radiometer  systems 
used  by  the  Visibility  Laboratory  team  during  OPAQUE 
IV. 

The  fly-away  ground  station,  which  again  included 
both  the  CRM  and  integrating  nephelometer  systems,  was 
operational  in  the  vicinity  of  all  four  of  the  primarj 
OPAQUE  flight  tracks  While  the  aircraft  operated  out  of 
Sigonella  Naval  Air  Facility  the  ground  station  was  located 
at  Birgi  Air  Base  near  Trapani.  Sicily  This  base  was  the 
location  for  the  Italian  OPAQUE  data  site  When  airborne 
operations  moved  to  Memmingen  Air  Base,  the  ground 
station  was  relocated  to  Birkhof,  Germany  which  is  the 
location  of  the  German  OPAQUE  data  site  Only  in  the 
United  Kingdom  was  the  fly-away  ground  station  not 
located  on  the  official  OPAQUE  site  While  the  aircraft 
was  operating  out  of  Mildenhall  Air  Base  the  ground  sta¬ 
tion  was  located  at  HMS  Heron,  Yeovilton  Royal  Naval 
Air  Station,  Yeovilton,  England  This  site,  approximately 
35  miles  northwest  from  the  British  OPAQUE  site  at 
Christchurch,  was  chosen  for  its  proximity  to  the  actual 
flight  track  and  for  its  suitability  for  air  transport  of  the 
ground-based  equipment  to  and  from  the  site.  When  air¬ 
borne  operations  returned  to  Germany  out  of  Wunstorf 
Air  Base,  the  ground  fly-away  ground  station  was  relocated 


to  the  Meppen  test  range  which  was  the  site  of  the 
United  States  supported  OPAQUE  data  site. 

In  spite  of  the  bad  weather  conditions  that  were 
anticipated  for  this  Winter  deployment,  OPAQUE  IV  was 
the  most  productive  of  all  the  OPAQUE  series.  The 
weather  was  exceptionally  permissive  as  is  indicated  by  the 
completion  of  twenty  five  missions.  The  spectrum  of 
weather  conditions  encountered  was  exceptionally  broad, 
ranging  from  reasonably  cloud  free  with  moderate  haze  to 
low  ceilings  with  rain  and  snow. 

OPAQUE  V.  SUMMER  1978 
FLIGHT  NOS.  C-460  THROUGH  C-^9 
AFGL-TR- 80-0207 

This,  the  second  Summer  deployment,  similar  in 
season  to  OPAQUE  111.  began  with  the  Visibility  Labora¬ 
tory  team  departing  San  Diego  on  Friday,  28  July  1978 
The  team  returned  on  Saturday.  30  September  1978  Dur¬ 
ing  the  64  day  intervening  period,  the  joint  Air  Force  and 
Visibility  Laboratory  team  operated  the  aircraft  out  of  six 
different  air  bases  in  four  different  countries  This 
represents  the  largest  number  of  bases  ever  utilized  during 
a  single  deployment  However,  since  all  of  these  bases 
had  been  visited  before,  there  were  no  operational 
difficulties  and  the  exceptional  cooperation  of  the  host 
nation  personnel  made  the  exercise  a  remarkably  smooth 
and  efficient  experience  The  chronological  listing  of  the 
operational  bases  utilized  during  this  final  OPAQUE 
deployment  is  shown  below: 


Ofilci 

Visited  Slatting  Base  Idem  Geographical  Location 

1  Rhcm  Main  A't  Base  near  Franklurt.  Germans 

2  Sigonella  Naval  Air  Facility  near  Catania.  Sicils 

3  Wunstorl  Air  Base  near  Hannover.  Germans 

4  Memmingen  Air  Base  S  of  Munich.  Germans 

5  Mildenhall  Air  Base  N  W  ol  Cambridge.  England 

b  Vaerlose  An  Base  neat  Copenhagen  L>enmark 


There  were  no  instrumentation  modifications 
imposed  during  the  preparations  for  the  OPAQUE  S' 
deployment,  only  normal  maintenance  and  overhaul.  The 
fly-away  ground  station  was  ojverational  at  or  near  tour  of 
tne  five  primary  OPAf/.'E  flight  track--  As  during  ihe 
preceding  OPAQUE  IV  deployment,  it  was  in  operation  on 
the  official  OPAQUE  data  site  ai  Trapani/ Birgi  Air  Base. 
Sicily,  and  at  both  the  Meppen  and  Birkhof  sites  in  Ger¬ 
many.  Similarly,  ihe  ground  station  was  located  at  Yeovil¬ 
ton  while  the  aircraft  was  operating  oui  of  Mildenhall. 
England 

The  weather  was  much  more  cooperative  during  this 
second  Summer  deployment  than  it  had  been  during  ihe 
first.  OPAQUE  III  Even  though  there  were  lew  truly 
clear  days,  there  were  many  marred  only  by  minor  cloud 
and  ha/c  conditions  Consequently,  of  the  nineteen  total 
missions,  ten  of  them  were  classified  as  "full  profile"  mis¬ 
sions.  i.c.  four  straight  and  level  altitudes  in  a  total  of 
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four  spectral  bands.  The  remaining  nine  were  assorted 
flights  of  lesser  quality  ranging  all  the  way  down  to  5000 
ft.  ceilings  with  haze  and  light  rain. 

In  a  manner  similar  to  that  used  in  Duntley  ei  al. 
(1978).  a  tabular  summary  of  the  pertinent  descriptive 
information  related  to  the  OPAQUE  IV  and  V  data  flights 
has  been  abstracted  from  the  parent  data  report  covering 
each  of  these  two  deployments  (AFGL-TR-79-0159  and 
AFGL-TR-80-0207).  and  reproduced  as  Tables  2  1  and 


2.2.  The  times  listed  under  the  "Total  Time  of  Data  Tak¬ 
ing"  column  are  Greenwich  Mean  Time  (GMT)  and  Local 
Civil  Time  (LCT).  The  sun  zenith  angles  are  tabulated 
for  the  time  when  data  taking  began,  at  the  time  of  sun 
transit  (minimum  sun  zenith  angle),  and  at  the  conclusion 
of  the  last  data  taking  event.  The  maximt  n  and 
minimum  flight  altitudes  in  meters  AGL  (Above  Ground 
Level)  are  noted  in  columns  12  and  13.  The  filter 
identification  numbers  refer  to  the  spectral  responses  illus¬ 
trated  in  Fig.  1-2. 
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Table  2.2  (eon*!.).  Right  Data  Summary  for  OPAQUE  V 
(AFGL-TR-80-0207) 


|  Total  Time  of  Data  Taking  | 


- 

1  Sun  2emlh  Angle 
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End 
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Table  2..1.  OPAQUE  Related  Aircraft  Data  Refxrrls 


Report  No 


Deployment 


Data  Type 


AFGL 

TR-77-0078  1 

l 

AFGL 

TR-77-0239  i 

II 

AFGL 

TR-78  0168  | 

111 

AFGL 

TR- 79-01 19  | 

IV 

AFGL 

TR-80-0207  j 

V 

AFGL 

TR-80-0192  i 

I. 11.111 

AFGL 

TR8I0154 

IV  V 

AFGL 

TR-81  0237  j 

l.ll. ill 

AFGL 

TR-81  0271  ! 

II. III.  v 

AFGL 

ril-8l-0.lt  7 

l-V 

Spring  76 
i  Fall  76 
i  Summer  '77 
J  Winter  '78 
I  Summer  '78 
I  See  Above 

i 


Scattering  Coefficient  Profiles  &  Related  Meteorology 
Scattering  Coefficient  Profiles  &  Related  Meteorology 
Scattering  Coefficient  Profiles  A  Related  Meteorology 
Scattering  Coefficient  Profiles  A  Related  Meieoroi'-ev 
Scattering  Coefficient  Profiles  A  Related  Meteorology 
Aerosol  St/e  Ihsirtbuiions  and  Analysis 
Scattering  Coefficient  Very  Low  Altitude  Profiles 
Scattering  Coefficient  Very  Low  Altitude  Profiles 
Sks  and  Terrain  Radiances 
Variations  m  Sky  and  Terrain  Radiances 


2.2  Selected  Sample  Data 

Flight  data  resulting  from  the  five  OPAQUE  deploy¬ 
ments  have  been  presented  in  a  series  of  AFGL  Technical 
Reports  which  are  identified  in  Table  2.3.  In  general  the 
data  presentation  reports  are  of  three  types.  The  first  con¬ 
taining  primarily  scattering  coefficient  profiles,  the  second 
sky  and  terrain  radiances,  and  the  third  containing  aerosol 
size  distributions 

In  order  to  provide  a  quick  overview  of  the  data 
processed  subsequent  to  that  presented  in  Dumley  el  at. 
(1978),  portions  of  the  original  OPAQUE  data  reports 
have  been  reproduced  and  presented  herein  with  addi¬ 
tional  but  typical  descriptive  and  analytic  comment.  For 
this  summary  presentation,  only  composite  plots  of 
pseudo-photopic  volume  scattering  coefficient  and  tem¬ 
perature  which  have  been  sorted  by  geographical  location 
are  included  as  representative  of  the  profile  type  reports. 
In  a  similar  manner,  several  composite  representations  of 
the  sky  and  terrain  radiance  plots  are  also  included  For 
moa  complete  information  regarding  these  data,  the 
reader  is  referred  to  the  parent  reports  as  listed  in  Table 
2.3. 

Temperature 

Temperature  is  generally  a  function  of  latitude  and 
season  with  an  expected  diurnal  and  weather  variability 
during  any  one  season  at  a  given  latitude  The  tempera¬ 
tures  measured  during  each  (light's  several  profile  ele¬ 
ments  are  averaged  by  altitude  into  a  composite  plot  which 
typifies  that  individual  flight  These  resultant  average 
temperature  profiles  for  the  OPAQUE  IV  and  V  deploy¬ 
ments  have  been  sorted  into  sets  representing  similar  lati¬ 
tudes  and  presented  in  Figs  2-1  and  2-2  These  in-flight 


temperatures  can  be  profitably  compared  to  temperatures 
from  the  U  S  Standard  Atmosphere  Supplements  (1966) 
which  are  shown  as  solid  lines  on  the  plots  in  Figs.  2-1 
and  2-2.  It  should  be  noted  that  the  altitude  scale  in  these 
figures  is  in  meters  above  mean  sea  level  (m  MSL).  and 
the  ground  elevations  at  the  test  sites  range  from  0  meters 
at  Trapani.  Sigonella.  Mildenhall  and  Rodby,  to  762  m  at 
Birkhof 

As  noted  in  the  parent  reports,  the  temperature 
profiles  are  reasonable  for  the  season  and  latitude  of  the 
measurement,  although  there  are  significant  examples  of 
low  and  mid-level  inversions  at  all  sites 

Total  Volume  Scattering  Coefficient 

Total  volume  scattering  coefficients,  as  measured 
with  the  filter  4  pseudo-photopic  spectral  response  for  the 
OPAQUE  IV  and  V  deployments,  are  presented  in  Figs. 
2-3  and  2-4  These  plots,  which  have  been  abstracted 
from  the  previously  issued  reports  listed  in  Table  1.2, 
have  been  sorted  in  a  manner  similar  to  the  preceding 
temperature  plots  i.e.  by  approximate  latitude  and  season 
One  should  note  that  the  altitude  scales  in  Figs.  2-3  and 
2-4  are  in  meters  above  ground  level  (m  AGL)  rather 
than  on  the  MSL  scale  used  for  the  temperature  plots. 
Also  note  that  the  horizontal  scale  spans  three  decades  for 
all  of  the  plots  except  2-3d,  which  spans  only  two. 

As  discussed  in  the  OPAQUE  IV  parent  report,  the 
degree  of  fine  structure  in  the  nephelometer  data  that  can 
be  attributed  to  instrumentation  and  system  noise  is  not 
more  than  ±5%  of  the  reading.  However,  as  illustrated 
by  the  profiles  representing  flight  C-432  in  Fig.  2-3a,  a 
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E  SEA  LEVEL  (KM)  ALTITUDE  ABOVE  SEA  LEVEL  (KM) 


OPAQliE  V  SUMMER  1978 


YEOVILTON  50°  56'  N 


RODBY  54°  51'  N 


TEMPERATURE  (DEG  C) 


Fig.  2-2.  (can’t.)  Temperature  for  OPAQUE  V  flights  2  August  to  2h  September  1978.  compared  to  temperature  from  U  S  Standard 

Atmosphere  Supplements  (l%6)  (AF<il_-TR-80-0207). 


semi-square  wave  biasing  voltage  was  imposed,  during 
some  flights,  upon  the  data  signal.  This  distortion  of  the 
profile  data  is  of  such  a  nature  that  the  erroneous  values 
are  those  representing  the  smaller  value  of  scattering 
coefficient.  These  distortions  appear  to  affect  only  the 
flights  in  Sicily  and  no  determination  of  their  cause  has 
been  completely  established.  Users  should  be  aware  of 
these  circumstances  during  their  evaluation  of  the  suitabil¬ 
ity  of  these  date  for  their  specific  application. 

Users  should  also  be  advised  that  the  potential  stray 
light  error  that  was  discussed  in  each  of  the  earlier  reports 
listed  in  Table  1.2  has  not  been  corrected  in  the  data 
displays  of  Figs.  2-3  and  2-4,  except  for  the  Trapani  data 
shown  in  Fig.  2-3a.  For  a  complete  discussion  of  this  cir¬ 
cumstance  the  reader  is  referred  to  Johnson  and  Gordon 
(1979). 

Several  general  characteristics  of  the  profiles  illus¬ 
trated  in  Figs.  2-3  and  2-4.  briefly  alluded  to  in  the  parent 
reports,  are  repeated  here.  In  both  the  Winter  and  Sum¬ 
mer  data  sets,  the  Trapani  profiles  are  clearly  different 
from  their  companion  profiles  measured  in  northern 
Europe.  In  the  Winter  set,  the  Trapani  profiles,  though 
noisy,  show  significantly  less  low  altitude  structure  than 
any  of  the  other  sets,  and  have  a  complete  lack  of  any 
abrupt  haze  top.  Similarly  the  summer  Trapani  profiles 
show  little  structure  and  additionally  maintain  relatively 
high  scattering  coefficient  values  throughout  the  entire  6 
km  altitude  band. 


It  appears,  in  general,  that  there  is  substantially 
more  structural  irregularity  in  the  scattering  coefficient 
profiles  during  the  Winter  set  than  during  the  Summer. 
This  is  particularly  true  for  the  northern  European  sta¬ 
tions,  Soesterberg,  Meppen,  and  Rodby  as  illustrated  in 
Fig.  2-3d  and  Figs.  2-4b,  d  and  f. 

A  necessary  but  unfortunate  artifact  of  the  profile 
data  discussed  in  the  preceding  paragraphs  is  that  the 
measurements  all  terminate  at  some  significant  altitude 
above  ground  level.  This  being  a  necessary  condition 
imposed  by  the  safety  of  flight  regulations  which  apply  to 
all  civil  air  space,  and  an  unfortunate  experimental  con¬ 
straint  in  view  of  the  extreme  sensitivity  of  slant  path  con¬ 
trast  transmittances  to  variations  in  the  near  surface  haze 
conditions.  The  data  contained  in  two  of  the  reports  listed 
in  Table  2.3  (AFGL-TR-81-0154  and  -0237),  were  pro¬ 
vided  to  reduce  the  uncertainties  in  our  knowledge  of  the 
scattering  coefficients  near-surface  profile. 

Several  samples  of  these  near  surface  profiles  are 
presented  as  Figs.  2-5  and  2-6.  There  are,  in  general,  only 
occasional  variations  of  any  substance  within  these  low 
altitude  profiles.  In  this  sample  of  54  low  altitude  daytime 
profiles,  acquired  in  clear  air  away  from  clouds  (i.e.  under 
VFR  conditions),  47  (87%)  are  represented  by  nearly  con¬ 
stant  scattering  coefficients  as  a  function  of  changing  alti¬ 
tude,  These  stable  conditions  are  clearly  illustrated  in 
Figs.  2-5a  and  2-6b.  Four  of  the  seven  anomalous  profiles 
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OPAQUE  V  SUMMER  1978 


TRAPANI  37°  33'  N  SOESTERBERG  51°  56'  N 


MEPPF.N  53°  00'  N  BfRKHOFF  4S°  15'  N 


Fig.  2-4.  Airborne  nephelonieter  pseuilo-phoiopiv  il.n.i  lor  OPAQUE  V  flight-.  Irom  2  AuguM  m 
26  September  I67H  I  AKil  -  TR -80-112071 
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OPAQUE  V  SUMMER  1978 


YE0V1LT0N  50°  56'  N 


TOTAL  VOLUME  SCATTERING  COEFFICIENT  (PER  MI 


RODBY  54°  51'  N 


Fig.  2-4.  (con't.)  Airborne  nephelomeier  pseudo-phot  opic  data  lor  OPAQUE  V  flights 
From  2  August  to  26  September  1978  (AEGL-TR-80-0207) 


that  were  measured  are  contained  in  the  illustrations  of 
Fig.  2-5b  and  2-6a.  As  noted  in  the  earlier  reports,  the 
identification  of  the  conditions  resulting  in  the  seven 
profiles  showing  variations  within  their  vertical  structure 
will  require  additional  analysis.  In  this  regard,  a  supple¬ 
mentary  set  of  precision  local  meteorological  observations 
including  local  trajectories  would  (y*  beneficial 

Sky  and  Terrain  Radiances 

As  noted  in  Table  2.3,  there  were  two  general 
classes  of  radiometric  data  collected  during  a  typical  air¬ 
borne  measurement  sequence  The  first  being  primarily 
scattering  coefficient  profile  data  as  illustrated  in  the 
immediately  preceding  figures  and  comments,  and  the 
second  being  measurements  of  sky  and  terrain  radiances  at 
selected  assorted  altitudes  This  second  complementary 
data  set,  introduced  in  Johnson  (1981).  includes  angularly 
precise  spectral  measurements  of  the  4n ■  radiance  field 
surrounding  the  aircraft  at  several  different  altitudes 
throughout  the  zero  to  six  kilometer  altitude  regime 

These  radiance  data  enable  one  to  characterize  a 
broad  variety  of  the  environmental  conditions  extant  dur¬ 
ing  the  flight  episode  and  thus  lead  directly  to  the  develop¬ 
ment  of  operationally  useful  predictive  models  The  radi¬ 
ance  data,  in  conjunction  with  their  companion  scattering 


coefficient  data,  are  readily  applicable  to  the  determination 
of  slant  path  contrast  transmittances,  atmospheric  optical 
depths,  aerosol  directional  scattering  characteristics,  flux 
divergences  and  their  attendant  determinations  of  turbid 
atmosphere  single  scattering  albedos. 

For  this  brief  overview,  only  a  small  sample  of  the 
nearly  500  radiance  arrays  that  have  been  processed  have 
been  selected  for  presentation  here.  As  with  the  preceding 
temperature  and  scattering  coefficient  data,  these  sky  and 
terrain  radiance  displays  have  been  abstracted  from  a  pre¬ 
viously  published  parent  reports.  Johnson  (1981),  and 
Johnson  and  Hering  (1981) 

The  data  displays  presented  in  Figs  2-8  through 
2-11  are  in  the  Composite  Sky-Terrain  Radiance  format 
used  in  several  earlier  reports  In  these  composite  plots, 
the  radiance  variations  throughout  the  sun-zenith  plane,  in 
both  the  up-sun  and  down-sun  directions,  are  plotted  as  a 
single  curve,  and  labeled  "PHI  —  0  &  -180"  Similarly,  the 
data  in  the  cross-sun  plane  are  plotted  continuously  and 
labeled  "PHI  -  90  &  -270".  In  all  cases  a  zenith  view 
angle  of  0°  represents  the  point  in  the  sky  directly  over¬ 
head,  and  a  zenith  view  angle  of  180°  represents  the  point 
on  the  underlying  terrain  directly  below.  The  negative 
signs  associated  with  the  angular  indices  are  an  artifact  of 
the  computer  sorting  technique,  and  should  not  be 
rigorously  interpreted. 
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As  an  aid  to  the  interpretation  of  the  radiance  distri¬ 
butions  illustrated  in  Figs  2-8  through  2-11.  a  pictorial 
representation  of  the  measurement  field  is  shown  in  Fig 
2-7  This  artist's  conception  of  the  4rr  radiance  field 


which  surrounds  the  aircraft  is  annotated  with  the  angular 
notations  and  terminologies  used  in  the  data  displays  and 
discussions  which  follow. 
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Six  radiance  plots  representing  poriions  ol  three 
different  data  (lights  have  been  absiracied  from  Ihe  parent 
reports  referenced  above  and  reproduced  as  Fig  2-X 
Ihese  si\  displays  are  the  standard  computer  output  lor 
d.ita  ol  this  type  and  lhu«  present  one  data  point  at  evert 
lire  degree  increment  in  zenith  angle  Radiances  meas¬ 
ured  in  lire  near-sun  regions  which  eo  off-scale,  as  in  lag 
2-  be.  have  been  auiomaiicallv  truncated  at  the  maximum 
plot  boundary  i.t  111  »  If  m  urn  .  and  thus  should  not 


be  imerpreicd  ;is  reliable  measurements  Neither  have 
ihese  data  been  corrected  lor  the  near-sun  stray  light 
errors  discussed  in  Johnson  (1981)  Therefore  radiance 
measurements  made  within  ±25°  from  the  solar  disc  wall 
contain  a  stray  light  component  which  is  dependent  upon 
the  aerosol  loading  ol  the  atmosphere  above  the  measure¬ 
ment  altitude  As  noted  in  Johnson  (1981 1.  the  near-sun 
stray  lighi  effects  can  range  from  zero,  as  in  the  full  over¬ 
cast  case  ol  Fig  2-8a.  to  as  much  as  a  factor  of  28  in  cases 
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similar  to  the  clear  high  altitude  situations  shown  in  Fig 
2-8d. 

It  is  clear  however,  that  even  with  the  above 
caveats,  there  are  several  readily  distinguishable  charac¬ 
teristics  inherent  in  the  data  plots  of  Fig  2-8. 

The  data  in  Fig.  2-8a  represents  one  of  the  most 
symmetric  radiance  distributions  retrieved  nom  the 
OPAQUE  series.  This  example,  one  of  eight  in  the  parent 
report,  illustrates  radiances  in  all  four  cardinal  azimuths 
from  the  sun  showing  nearly  identical  gradients  as  a  func¬ 
tion  of  zenith  angle  with  very  little  sky  structure  One 
should  note  the  relatively  large  radiance  difference 
between  the  overhead  sky  and  the  underlying  ocean  sur¬ 
face,  and  contrast  this  difference  with  the  data  in  Fig. 
2-8b.  In  Fig.  2-8b  the  data  represent  a  similar  o\ercas> 
sky  condition  over  a  completely  different  underlying  ter¬ 
rain.  The  relatively  high  reflectance  of  the  snow  covered 
plateau  underlying  the  Birkhof  flight  track  has  produced 
the  most  directionally  constant  radiance  levels  retrieved 
from  the  OPAQUE  series. 

The  characteristics  of  radiance  uniformity  and  sym¬ 
metry  illustrated  by  the  data  from  flights  C-422  and  C-435 
are  relatively  unique  and  can  be  considered  to  represent 
one  extreme  in  the  overall  spectrum  of  sky  and  terrain 
radiance  distributions. 

Just  as  the  fully  overcast  conditions  illustrated  in 
Figs.  2-8a&b  illustrate  one  end  of  the  meteorological  spec¬ 
trum,  the  clear  sky  conditions  of  flight  C-379,  Figs. 
2-8cAd,  represent  the  other.  These  plots,  two  of  the  four¬ 
teen  which  appear  in  the  parent  report,  represent  clear  sky 
conditions  whose  radiance  patterns  are  typically  well 
behaved  and  tell  a  consistent  story  with  respect  to  upper 
and  lower  hemisphere  characteristics,  both  as  a  function  of 


observational  wavelength  and  altitude.  The  bright  hor¬ 
izons  and  highly  asymmetric  up-sun  and  down-sun  radi¬ 
ance  levels  are  consistent  and  clearly  discernible.  For  this 
over  water  track,  the  uniformity  of  the  down-sun  and 
cross-sun  sea  surface  radiances  is  well  illustrated  as  is  the 
obvious  glitter  pattern  which  appears  as  a  broad  high  radi¬ 
ance  peak  in  the  up-sun  data  of  Fig.  2-8c  and  narrows  as 
the  flight  altitude  increases  to  6  km  in  Fig.  2-8d.  The  rela¬ 
tively  high  radiance  of  the  nearby  shoreline  gradually 
appears  in  the  down-sun  lower  hemisphere  measurements 
as  the  altitude  of  observation  increases,  also  shown  clearly 
in  Fig.  2-8d. 

Right  C-468  provided  fifteen  data  plots,  two  of 
which  are  shown  here  as  Figs  2-8e&f.  These  data  were 
measured  on  a  nominally  clear  day  that  was,  however, 
dominated  by  high,  thin  cirrus  The  subtlety  and  difficulty 
of  precisely  identifying  cirrus  effects  are  illustrated  by 
comparing  plots  c&d  with  plots  e&f  In  this  example,  both 
of  the  low  altitude  plots  c&e  show  very  similar  cross-sun 
and  down-sun  patterns  although  the  Rodby  data  (c)  was 
taken  under  clear  sky  and  the  Meppen  (e)  under  thin 
cirrus.  In  the  high  altitude  pair,  the  Rodby  data  (d)  show 
the  irregular  cross-sun  and  down-sun  characteristics  of  the 
low  altitude  data,  apparently  cirrus  contaminated,  while 
the  Meppen  plots  (f)  show  the  regular  and  well  matched 
characteristic  normally  anticipated  from  high  altitude  clear 
sky  measurements 

Several  additional  comparisons  of  sky  and  terrain 
radiance  characteristics  are  illustrated  in  Figs.  2-9.  2-10 
and  2-11  which  have  been  abstracted  from  Johnson  and 
Hering  (1981).  In  these  composite  plots,  several  flight 
sequences  have  been  combined  on  a  single  display  in 
order  to  more  clearly  illustrate  the  characteristic  of 
interest  The  three  characteristics  illustrated  in  this  partic¬ 
ular  senes  are  radiance  variations  as  a  function  of  spectral 
passband.  variations  as  a  function  of  altitude  and  varia¬ 
tions  as  a  function  of  cloud  cover. 

Comparative  measurements  of  the  sky  and  terrain 
radiance  fields  as  a  function  of  wavelength  for  flight  C-466 
are  shown  in  Fig.  2-9.  Since  the  coefficient  for  Rayleigh 
U.e.  molecular)  scattering  varies  very  nearly  as  the  -4th 
power  of  the  wavelength,  A.  and  the  Mie  (i.e.  aerosol) 
scattering  coefficient  varies  roughly  in  the  range  A  0 '  to 
A  ",  the  radiance  field  has  a  strong  spectral  dependence. 
This  dependence  is  illustrated  clearly  in  Figs.  2-9a&b  by 
the  very  pronounced  decrease  in  sky  radiance  with  increas¬ 
ing  wavelength  through  this  visible  spectrum  interval. 
Note  also  that  the  spectral  variations  of  the  measured 
apparent  terrain  radiance,  while  heavily  cluttered  by  the 
influence  of  small  scale  fluctuations  in  surface  reflectivity, 
still  illustrate  significant  spectral  separations.  A  prominent 
feature  in  this  regard  is  the  relatively  large  apparent  terrain 
radiance  measured  in  the  750nm  spectral  band.  These 
high  radiances  result  from  the  characteristically  large  sur¬ 
face  reflectance  for  rural  farm  and  wooded  areas  in  the 
near  infrared  portions  of  the  spectrum. 

The  data  from  flight  C-466  also  can  be  used  to  illus¬ 
trate  the  variation  in  sky  and  terrain  radiances  as  a  func¬ 
tion  of  altitude.  These  altitude  dependent  characteristics 
are  presented  in  Fig.  2-10. 
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Fig.  2-8.  Selected  sky  and  terrain  radiance  plots  from  flights  C'-.VM.  C-422.  C-43S  and  C -468  Plots  a.  c.  and  c  are  low  altitude 
(  -  27m  AGL)  Plots  b.  d.  and  f  are  high  altitude  ( >  2S00  m  AGL) 
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The  influence  of  atmospheric  path  radiance  upon 
the  apparent  radiance  of  both  the  sky  and  terrain  is  clearly 
demonstrated  in  both  Figs.  2-10a&b.  At  the  wavelength 
chosen  for  illustration.  A  ~  475nm,  as  well  as  for  those 
not  shown  here  but  included  in  the  parent  report  Johnson 
and  Hering  (1981).  the  systematic  decrease  in  apparent 
sky  radiance  as  the  altitude  of  measurement  increases  is 
obvious,  as  is  the  corollary  increase  in  apparent  terrain 
radiance.  One  should  note  the  tendency  for  the  percent 
change  in  sky  radiance  between  specific  altitude  levels  to 
be  much  the  same  over  a  wide  range  of  viewing  angles  as 
is  discussed  from  an  analytic  point  of  view  in  the  parent 
report 

To  illustrate  the  variation  of  sky  and  terrain  radi¬ 
ances  as  a  function  of  variable  cloud  cover,  a  composite 
plot  containing  data  from  five  separate  flights  is  presented 
as  Fig  2-11  All  of  the  measurements  selected  for  this 
comparison  were  made  in  the  pseudo-photopic  passband 
i  f.  A  —  550nm.  see  Fig  1-2.  These  five  flight  episodes 
represent  a  broad  variety  of  conditions,  ranging  from 
clear.  to  overcast.  <$.  as  measured  from  an  intermedi¬ 
ate  altiu.de  of  approximately  1 150  m  AGL. 

One  notes  that  the  measured  sky  and  terrain  radi¬ 
ances  for  these  cloudy  cases  tend  to  fluctuate  both  above 
and  below  the  radiance  distribution  for  the  clear  sky.  see 
Johnson  and  Hering  (1981).  depending  upon  the  scatter¬ 
ing  and  reflectance  properties  of  the  clouds,  ha/e  and  ter¬ 
rain  The  measured  cloudy  sky  radiance  is  found  to  be  a 
factor  of  two  or  more  larger  in  some  instances,  but  in  gen¬ 
eral  remains  within  about  ±50  percent  of  the  clear  sky 
radiance  for  corresponding  paths  of  sight  It  is  appropriate 
to  keep  in  mind  the  significant  differences  in  terrain 
characteristics  for  the  flights  shown  in  Fig.  2- 1 1  a&b.  The 
Meppen  (MP).  Soesterberg  (SO),  and  Yeovilton  (YO) 
flights  were  made  over  cultivated  farmland,  and  the  two 
Birkhof  flights  were  made  over  a  snow  covered,  wooded 
plateau  The  influence  of  these  different  terrain 
reflectances  upon  their  associated  path  radiances  can  be 
significant  and  is  appropriate  for  study  using  fast  opera¬ 
tional  models  such  as  that  described  by  Hering  (|98|) 

2.5  Data  Evaluations 

In  the  preceding  sections,  a  variety  of  profile  and 
fixed  altitude  radiometric  measurements  of  total  volume 
scattering  coefficient  and  sky  radiances  have  been 
presented  as  an  overview  of  the  data  reductions  accom¬ 
plished  as  part  one  of  this  two  part  contract  Over  a 
period  of  several  years,  as  the  various  experimental  meas¬ 
urements  and  data  reduction  activities  progressed,  there 
were  intermittent  but  continuing  evidences  of  inst rumen- 
tally  induced  error,  in  one  system  or  another,  that  often 
were  persistently  unyielding  in  identifying  the  source  of 
the  suspected  error  In  most  cases  adequate  identification 
of  the  system  malfunction  was  accomplished  through  a 
study  of  the  automatic  diagnostic  tests  applied  to  the  data 
in  the  early  stages  of  processing  However,  there  were 
two  instances  where  an  uncertainly  in  the  measurements 
continued  throughout  the  program  Since  these  uncertain¬ 
ties  relate  (o  the  primary  instrument  systems  and  their 


data,  i.e.  the  integrating  nephelometer,  and  the  upper 
hemisphere  scanner,  some  further  comment  seems 
appropriate. 

Total  Volume  Scattering  Coefficient 

As  discussed  in  several  earlier  reports,  among  them 
Duntley  el  at.  (1975),  and  Duntley  era/.  (1978),  the  air¬ 
borne  integrating  nephelometer  underwent  a  major 
reconfiguration  shortly  before  the  initial  OPAQUE  deploy¬ 
ment  in  the  Spring  of  1976  As  an  outgrowth  of  this 
reconfiguration  the  essentially  in-line  flow  through  device, 
as  illustrated  in  Fig.  3-1  of  Duntley  era/.  (1975),  was 
converted  into  the  folded  path,  compact  design  illustrated 
in  Fig  3-4  of  Duntley  era/.  (1978)  There  were  several 
stray  light  problems  associated  with  this  new  design  which 
became  apparent  during  the  analysis  of  the  resultant  data, 
particularly  with  respect  to  the  high  altitude  data  collected 
well  above  the  primary  ha/e  layer  A  discussion  of  this 
problem  and  the  development  of  a  corrective  procedure  is 
outlined  in  Duntley  era/.  (1977).  A  similar  discussion 
related  to  the  data  from  each  of  the  five  OPAQUE  deploy¬ 
ments  was  included  in  each  of  the  five  OPAQUE  data 
reports  listed  at  the  bottom  of  Table  12 

Continuing  analysis,  as  referred  to  in  Johnson  and 
Gordon  (1980),  of  both  the  total  volume  scattering 
coefficient  measurements,  as  well  as  the  directional 
scattering  function  measurements  at  scattering  angles  /)  — 
30°  and  150°.  led  to  the  conclusion  that  in  most  cases,  the 
measurements  ol  the  upper  altitude  total  volume  scatter¬ 
ing  coefficients  were  too  high  Thus,  several  additional 
procedures  were  employed  to  try  and  specify  the  character 
of  the  error  and  to  develop  a  reasonable  corrective  pro¬ 
cedure 

Two  separate  data  sets  were  employed  in  the  evalua¬ 
tion  of  the  probable  error  in  the  total  volume  scattering 
coefficient  data.  The  first  was  that  containing  the 
nephelometers  contemporary  measurements  of  directional 
scattering  function  at  scattering  angles.  £f  =  30°  and  150° 
The  second  was  that  containing  the  sky  radiance  distribu¬ 
tions  illustrated  in  Fig  2-8.  which  also  were  made  simul¬ 
taneously  with  the  nephelometer  sets.  Two  procedures 
were  used  with  each  set.  and  were  essentially  unanimous 
in  their  predictions,  which  are  summarized  in  Table  2.4. 
and  discussed  briefly  in  the  following  paragraphs  In  both 
cases  the  procedures  were  conducted  in  terms  of  the  opti¬ 
cal  scattering  ratio  Q  where  £)(:)  is  defined  by  the  ratio  of 
the  total  volume  scattering  coefficient  at  altitude  r,  to  the 
Rayleigh  U.e.  Molecular)  volume  scattering  coefficient  for 
a  given  spectral  band  at  the  same  altitude  (Hering  (1981)) 

Predictions  Based  Upon  Directional  Scattering 
Function  Measurements 

Comparisons  of  directional  scattering  functions 
measured  using  the  Visibility  Laboratory  nephelometer. 
with  those  measured  by  Barteneva  (1960)  indicate  a  strong 
correlation  between  the  two  sets  of  measurements.  This 
relationship  has  been  illustrated  in  each  of  the  reports 
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Table  2.4.  Comparison  or  Optical  Scattering  Ratio 
Determinations  for  the  Region  6-10  km  AGL. 


ship  may  be  established  either  by  using  the  Barteneva  data 
whose  median  values  are  shown  in  Fig.  2-12,  or  by  estab¬ 
lishing  another  analytic  representation  such  as  the 
Henyey-Greenstein  functions  (Irvine,  1968)  which  were 
used  by  Hering  (1981).  Estimates  of  (Hz)  were  made 
using  each  of  these  representations  with  the  measured 
values  of  3(z,/}),  where  2(z,/3)  is  defined  by  the  ratio  of 
the  total  volume  scattering  function  at  altitude  z  and 
scattering  angle  /3,  to  the  molecular  (or  Rayleigh)  volume 
scattering  function  at  the  same  altitude  and  scattering 
angle,  and  are  listed  in  Table  2.4. 


listed  in  Table  2.1,  and  is  repeated  here  as  Fig  2-12, 
which  has  been  abstracted  from  Johnson  and  Gordon 
(1980).  Based  upon  these  data,  one  may  readily  accept 
the  premise  that  the  optical  scattering  ratio  Q(z)  can  be 
predicted  on  the  basis  of  measurements  of  the  volume 
scattering  function  ratio  (}(:.(3)  The  functional  relation- 


Predictions  Based  Upon  Sky  Radiance  Measurements 

Proceeding  from  the  development  by  Livshits  and 
Pavlov  (1970).  it  can  be  shown  that  aerosol  optical  thick¬ 
ness  may  be  obtained  from  measurements  of  solar  almu- 
cantar  radiances.  It  is  clear  therefore,  that  if  one  were  to 
use  the  high  altitude  i.c.  (6  km  AGL)  sky  radiance  data  as 
measured  by  the  airborne  scanner  system  described  in 
Johnson  (1981)  to  determine  the  aerosol  optical  thickness 
above,  and  if  one  makes  reasonable  assumptions  regarding 
the  composition  of  the  regions  above  10  km.  then  one 
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might  deduce  an  accurate  estimate  of  the  optical  scattering 
ratio  Q  within  the  6  km  to  10  km  regime.  Thus,  using 
measured  solar  almucantar  data  at  f) 55°  and  estimates  of 
the  almucantar  radiance  at  /3125°  where  necessary,  esti¬ 
mates  of  aerosol  optical  thickness  and  thence  estimates  for 
optical  scattering  ratio  were  determined  for  flights  C-466 
and  C-468.  The  estimated  values  also  appear  in  Table  2.4. 

It  is  a  straight  forward  procedure  to  determine 
atmospheric  transmittance  due  to  scattering  from  measure¬ 
ments  of  sky  radiance  using  the  radiance  ratio  technique 
of  Kushpil  and  Petrova  (1971).  This  technique  has  been 
used  often  by  the  Visibility  Laboratory,  in  its  computer¬ 
ized  form,  using  the  2 it  radiance  arrays  produced  by  the 
airborne  scanner  system.  These  transmittances,  if  deter¬ 
mined  from  the  sky  radiances  measured  at  6  km  AGL  and 
corrected  for  nominal  ozone  absorption  can  also  be  con¬ 
verted  to  optical  depth  (/.<>.  optical  thickness)  and  used  to 
estimate  the  optical  scattering  ratio  Q  much  as  in  the 
preceding  paragraph  This  alternate  technique  was  also 
used  with  the  data  from  flights  C-466  and  C-468  resulting 
in  the  additional  estimates  shown  in  Table  2  4 

It  is  clear  that  while  each  of  the  four  procedures 
outlined  above  are  derived  from  a  valid  and  fairly  rigorous 
theoretical  analysis,  each  also  has  in  its  application  a 
necessary  uncertainty  due  to  assumptions  primarily  regard¬ 
ing  the  upper  atmospheres  constitution  and  absorption 
properties,  and  the  radiometric  measurements  themselves 
It  is  however  equally  clear  that  the  consistency  in  their 
estimated  values  of  optical  scattering  ratio  give  strong 
credence  to  the  contention  that  the  measured  values  (bold 
face  m  Table  2.4)  are  consistently  too  high. 

An  attractive  and  technically  profitable  program 
would  result  from  an  increased  effort  in  the  development 
of  these  techniques  and  their  application  to  the  existing 
data  bases  representing  the  altitude  regime  defined  by  the 
6  km  AGL  measurements. 

Near-Sun  Sky  Radiances 

The  measurement  of  the  4n-  radiance  field  which 
surrounded  the  airborne  instrument  system  at  each  flight 
altitude  has  provided  a  data  base  of  wide  utility.  The  gen¬ 
eral  scope  of  the  4n-  radiance  data  base  has  been  discussed 
and  illustrated  bv  Johnson  and  Hering  (1981),  and  a  basic- 
fault  existing  in  each  of  the  arrays  due  to  stray  light 
influences,  has  been  addressed  in  Johnson  (1981 ) 

As  discussed  at  length  in  Johnson  (1981).  both  the 
comparison  of  irradiances  calculated  from  scanner  radiance 
arrays  with  directly  measured  irradiances,  and  comparison 
of  measured  near-sun  radiance  gradients  with  model  pred¬ 
ictions.  indicated  conclusively  that  the  near-sun  radiance 
measurements  were  consistently  too  high  Subsequent 
laboratory  measurements  indicated  that  off-axis  light 
entering  the  scanner's  optical  system  was  substantial  at  all 
viewing  angles  near  the  sun.  and  would  consistently 
increase  markedly  both  as  the  instrument's  field  of  view 
approached  the  solar  disc  and  as  the  atmospheric  aerosol 
loading  decreased 

The  implications  to  be  derived  from  the  measure¬ 


ments  just  discussed  were  that  high  altitude  or  clear  day 
measurements  of  sky  radiances  within  10°  of  the  solar  disc 
will  have  stray  light  ratios  of  as  much  as  20  to  35,  whereas 
low  altitude  near-sun  measurements  through  more  turbid 
atmospheres  will  have  substantially  lower  stray  light  ratios 
of  about  3  to  6  at  the  same  angular  distance  from  the  disc. 
There  were  no  detectable  stray  light  influences  upon  meas¬ 
urements  made  beyond  25°  from  the  disc  under  any  of  the 
conditions  tested.  Thus,  using  sky  radiances  to  predict 
optical  scattering  ratios  as  discussed  in  the  preceding  para¬ 
graphs  will  not  be  influenced  by  these  effects  since  most  of 
the  radiances  used  were  at  scattering  angles  of  55°  or 
more. 

Recommendation 

The  application  of  the  measured  atmospheric  optical 
properties  discussed  in  the  preceding  paragraphs  has  led  to 
the  development  of  a  reliable  and  fast  operational  model. 
Hering  (1981).  It  is  quite  appropriate  at  this  point  to  use 
the  model  in  conjunction  with  the  procedures  outlined  in 
Section  2.3  to  provide  a  fully  corrected  set  of  selected 
scattering  coefficient  profiles  and  their  associated  sky  and 
terrain  radiances.  These  optimized  data  sets  would  pro¬ 
vide  a  valuable  tool  in  the  development  of  next  generation 
models,  as  well  as  accurate  specifications  for  the 
identification  of  typical  meteorological  boundary  condi¬ 
tions 

3.  AIRBORNE  AEROSOL  MEASUREMENTS 
3.1  Introduction 

Aerosols  play  an  important  role  in  determining  the 
characteristics  of  skylight,  the  transmittance  of  light,  and 
other  optical  properties  of  the  atmosphere.  The  more 
rigorous  radiative  transfer  calculations  used  to  predict 
these  atmospheric  optical  properties  need  as  input  a  model 
of  the  aerosol  size  distribution.  Some  of  the  models  com¬ 
monly  used  have  been  the  modified  gamma  function,  used 
by  Deirmendjian  (1969)  to  construct  his  haze  models;  the 
power-law  size  distribution  function  of  Junge  (1963);  and 
the  log-normal  distribution  as  described  by  Aitchison  and 
Brown  (1957),  Raabe  (1971).  Davies  (1974).  and  Jaenicke 
and  Davies  (1976) 

To  define  an  aerosol  distribution  for  a  specific  situa¬ 
tion,  the  model  fit  parameters  must  be  expressed  in  terms 
of  predictable  physical  observables  such  as  altitude,  loca¬ 
tion,  and  the  meteorological  variables.  Many  research 
investigations  have  been  devoted  to  just  such  a  description 
of  a  maritime  aerosol  in  the  marine  boundry  layer.  The 
Munn-Katz  marine  aerosol  model,  presented  in  Hughes 
and  Richter  (1980),  and  the  parent  model  of  Wells  et  al. 
(1977)  are  examples  The  fit  parameters  of  the  Munn- 
Katz  model  are  expressed  as  functions  of  altitude,  relative 
humidity,  and  the  wind  speed  at  the  sea  surface.  The 
expense  and  difficulty  in  making  aerosol  measurements, 
however,  makes  it  difficult  to  investigate  the  behavior  of 
the  model  parameters  for  the  wide  variety  of  atmospheric 
conditions  that  can  occur  in  a  continental  or  maritime 
environment. 
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The  following  sections  present  a  discussion  of  an 
extensive  set  of  measurements  of  the  aerosol  size  distribu¬ 
tion  in  the  lower  troposphere  taken  during  the  spring  and 
fall  of  1976,  the  summer  of  1977,  and  the  summer  and 
winter  of  1978  at  eight  sites  in  Europe  and  Great  Britain 
The  data  are  used  to  explore  the  characteristics  of  the 
measured  aerosol  as  a  function  of  altitude,  season,  loca¬ 
tion,  and  relative  humidity.  Comparisons  of  the  measured 
distributions  are  made  to  some  of  the  more  commonly 
used  model  distributions. 

In  a  previous  report  by  Fitch  and  Cress  (1981), 
hereafter  termed  report  1,  the  characteristics  and  behavior 
of  the  aerosol  volume  distribution  were  explored  for  part 
of  the  almost  600  total  measured  distributions.  The 
analysis  of  report  1  concentrated  on  the  characteristics  of 
the  submicrometer  volume  mode,  showing  that  it  was  well 
fit  by  a  log-normal  curve  and  that  the  behavior  of  the  fit 
parameters  was  regular.  In  this  report,  the  entire  data  set 
is  analyzed  and  the  characteristics  of  the  entire  measured 
volume  distribution  explored. 

3.2  Method  of  Measurement 

Aircraft  measurements  of  the  number  and  size  dis¬ 
tribution  of  the  atmospheric  aerosol  were  conducted  along 
tracks  of  constant  altitude  above  the  ground  near  Yeovil 
and  RAF  Mildenhall  (Great  Britain).  Rodby  (Denmark). 
Meppen.  Ahlhorn.  and  Birkhof  (W  Germany),  Soester- 
berg  (Netherlands),  Bruz  (France),  and  Trapani  (Italy) 
The  Ahlhorn  data  was  included  in  the  Meppen  data  for 
this  study  because  of  the  close  proximity  of  the  two  air¬ 
craft  tracks  to  one  another,  approximately  15  kilometers 
All  of  the  50  km  long  sampling  tracks  the  locations  of 
which  are  shown  in  Fig.  1-3  and  tabulated  in  Table  3  2. 
were  over  land,  except  for  the  Rodby  track  which  was 
over  the  Baltic  Sea.  Trapani  which  was  over  the 
Mediterranean  Sea.  and  the  Mddenhall  (rack  which  was 
over  the  North  Sea  The  procedure  and  the  instrumenta¬ 
tion  used  in  the  measurement  program  have  been 
described  in  detail  in  report  1  and  by  Cress  (1980).  and  so 
will  only  be  briefly  reviewed  here 

The  aerosol  number  distribution  was  measured  over 
particle  radii  from  0  2  to  5.9/im  using  an  89  channel 
Royco  220  single  particle  counter  carried  on  board  the 
C-130  aircraft  illustrated  in  Fig  I- 1  Samples  were  taken 
in  cloud  free  air  at  altitudes  ranging  from  150  to  6100 
meters  (m)  above  the  surface.  The  sampling  time  of  the 
Royco  was  typically  four  minutes,  but  occasionally  samples 
were  taken  for  one  minute  or  as  long  as  ten  minutes, 
depending  on  the  situation  and  the  density  of  particles 
being  counted.  Sample  air  for  the  Royco  was  provided  at 
about  47  cubic  centimeters  per  second  by  a  ram  air. 
approximately  isokinetic  probe  located  on  top  of  the  air¬ 
craft  fuselage  33  cm  above  the  aircraft  skin  The  entire 
Royco  system  consisted  of  the  Royco  counter,  a  Technical 
Measurement  Corporation  Model  102  Gamma  Scope  II 
Pulse  Height  Analyzer,  and  a  Hewlett-Packard 
Clock/ Printer 

The  Royco  is  a  right-angle  scattering  counter  that 


was  calibrated  with  submicrometer  size  polystyrene 
spheres.  Since  the  instrument  is  a  right-angle  scattering 
device,  a  major  source  of  error  is  an  inadequate  knowledge 
of  the  index  of  refraction  of  the  ambient  aerosol.  The 
refractive  index  is,  unfortunately,  difficult  to  measure  and 
so  is  rarely  known  in  measurements  of  size  distributions. 
The  error  in  measured  particle  size  associated  with  the 
effects  of  refractive  index  can  be  estimated  from  labora¬ 
tory  measurements  for  an  average  range  of  probable 
indices.  In  such  an  investigation  Quenzel  (1969)  found  a 
maximum  sizing  error  of  a  factor  of  2  to  3.  A  ground 
based  intercomparison  study  between  a  Royco  220  and 
active  and  classical  Particle  Measurement  Systems  aerosol 
counters  showed  agreement  to  within  a  factor  of  2  to  5  in 
particle  concentration  (Cress,  1980) 

Another  source  of  error  may  have  been  the 
modification  of  the  physical  properties  of  the  ambient 
aerosol  sample  as  it  passed  from  outside  the  aircraft  to  the 
Royco  scattering  chamber.  Possible  modifications  to  the 
aerosol  sample  could  have  been  the  result  of  a  drying  of 
the  particles  or  a  selective  loss  of  particles.  The  small 
number  of  changes  in  the  direction  of  flow  within  the 
plumbing  probably  minimized  the  effects  of  particle  loss. 
According  to  the  investigations  of  report  I  and  Cress 
(1980),  the  characteristics  of  the  measured  aerosol  prob¬ 
ably  represent  the  influence  of  a  relative  humidity  of  lower 
value  than  the  measured  ambient  value.  Although  this 
could  have  reduced  the  sizes  overall,  the  shape  of  the 
major  features  of  the  distribution  were  likely  preserved 
due  to  the  small  size  range  measured.  Comparison  of 
consecutive  one  minute  data  samples  indicate  that  the 
repeatability  of  the  airborne  Royco  was  good. 


3.3  Results  of  Measurements 

For  many  years  size  distribution  data  has  been 
presented  in  figures  of  log(rf,V/</  log  r )  versus  log  r  where 
r  is  the  particle  radius  and  N  is  the  total  number  density 
of  particles  smaller  than  r.  The  measured  distributions 
are  generally  compared  to  a  common  size  distribution  for¬ 
mula  such  as  the  one  developed  by  Junge  (1963): 
dSid  log  r  =  it  "  where  <  is  a  constant  depending  on  the 
total  number  of  particles  and  d  is  estimated  to  be  about  3 
for  particle  radii  between  0.3  and  10/xm,  (Junge  and  Jaen- 
icke.  1971)  This  power-law  size  distribution  was  fit  to  the 
data  collected  in  Great  Britain  and  F.urope  and  the  results 
presented  in  Cress  (1980).  The  investigation  shows  that, 
overall,  the  value  of  /3  is  in  the  range  3  to  6  for  data  col¬ 
lected  within  the  mixing  layer  and  has  a  value  close  to  2 
for  data  collected  above  the  layer 

Departures  of  the  measured  size  distribution  from 
the  power-law  distribution  appear  small  on  the  four  or  five 
log  cycles  traditionally  used  in  figures  of  aerosol  number 
distribution.  These  small  departures  were  found,  how¬ 
ever.  to  yield  multi-mode  curves  when  the  data  were  plot¬ 
ted  as  volume  (or  mass)  versus  particle  size,  as  pointed 
out  by  Whitby  el  al.  (1972)  and  Shettle  (1975).  The 
appearance  of  a  volume  mode  is  demonstrated  in  Fig.  3-1 
which  presents  an  aerosol  distribution  measured  at  Mep- 
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fin-  3-1.  A  comparison  ol  the  aerosol  number  and  solumc  disin- 
hulions  as  a  I  unci  ion  of  particle  radius  lor  dala  taken  al  Meppen 
during  August  at  an  altitude  of  457m  above  the  ground  The  lit 
of  a  power-law  distribution  with  a  slope  of  -3  to  the  number  distri¬ 
bution  Idashed  line)  and  the  lit  of  Iq  I.)  It  to  the  volume  disiri- 
button  I  +  >  are.  also,  shown 


pen  during  ihe  month  of  August  al  an  altitude  of  457m 
above  the  ground  The  data  is  displayed  in  terms  of  a 
number  distribution  i/.V/t/  logr  and  a  volume  distribution 
c H’/d  log  r  where  dl'.UI  logr  =  4/3  n  r'  dN/d  logr.  Only 
the  data  from  every  other  Royco  channel,  as  shown  by 
dots,  is  displayed  The  data  at  particle  radii  greater  than 
I nm  was  smoothed  using  an  equally  weighted  seven  point 
average  The  smoothed  data  is  shown  as  a  solid  line  in 
the  volume  display 

The  curves  of  Fig  3-1  demonstrate  that  departures 
of  a  measured  number  distribution  from  a  fitted  Junge 
distribution  with  /J-3  (dashed  line)  can  translate  into 
modes  when  the  data  are  plotted  as  a  volume  distribution 
A  Junge  distribution  with  /J=3  would  appear  as  a  straight 
horizontal  line  in  the  volume  display.  The  mode  in  the 
submicron  region  was  termed  the  accumulation  mode  by 
Whitby  (1978).  who  slated  that  its  main  source  of  mass  is 
from  the  coagulation  of  the  smaller  size  aerosols  and  from 
gas  to  particle  conversion  of  generally  anthropogenic  ori¬ 
gin  A  similarly  shaped  mode  centered  at  radii  greater 
than  3gm,  shown  later  in  Fig  3-2.  was  termed  the  coarse 
mode  by  Whitby  (1978)  and  is  considered  to  have  sources 
such  as  wind  blown  dust  and  sea  spray,  which  tend  to 
make  its  characteristics  independent  of  those  of  the  accu¬ 
mulation  mode  A  fit  of  the  volume  distribution  with  a 
log-normal  curve  (  +  )  is  shown  in  Fig  3-1  for  comparison 
and  will  be  discussed  in  the  next  section 

The  limited  number  of  investigations  into  the  nature 
of  an  aerosol's  volume  distribution  and  the  large  spread  in 
the  values  of  d  from  the  power-law  fit  to  the  data  (Cress. 
1980)  prompted  an  analysis  of  the  data  using  the  volume 
size  display  This  analysis  docs  not  represent  an  attempt 
to  define  a  new  aerosol  model  but  is  intended  to  demon¬ 
strate  the  detail  with  which  an  aerosol  distribution  can  be 
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Fig.  3-2.  The  volume  distribution  sampled  at  Soeslerherg  during 
August  at  an  aim ude  of  152m  (a)  and  at  Birkhof  during  February 
at  an  altitude  of  1219m  (b)  The  til  of  Eq  (I)  to  ihe  data  (  +  )  is 
shown  for  comparison 


described  using  the  volume  display.  The  first  part  presents 
a  mathematical  description  of  the  volume  distribution. 
Subsequent  sections  discuss  the  relationships  between  the 
fitting  parameters  and  how  the  aerosol  distribution 
behaves  as  a  function  of  altitude,  humidity,  season,  and 
site  of  measurement. 

3.4  Fitting  Function 

An  examination  of  volume  distributions  for  the 
entire  data  set  revealed  that  a  majority  of  the  modes  fall 
into  the  accumulation  or  coarse  particle  mode  categories. 
There  are  a  small  but  significant  number  of  cases,  how¬ 
ever.  when  the  volume  display  showed  the  presence  of  a 
mode  centered  between  the  regions  generally  assigned  to 
the  accumulation  and  coarse  particle  modes.  The  presence 
of  this  middle  mode,  necessitated  a  model  fit  consisting  of 
three  log-normal  distributions.  By  using  three  separate 
log-normal  curves  to  fit  the  data,  the  features  of  the  accu¬ 
mulation  and  coarse  particle  modes,  together  with  the 
characteristics  of  the  curious  middle  mode,  could  be  stu¬ 
died  individually.  The  model  used  to  fit  the  data  was 

=  I  if,  expl  -0  5 1 ( log /■— logr, )/«r, I J  I  +  k  (3.1) 

where  5  is  a  constant,  <r  the  standard  deviation,  f  the 
mode  radius,  and  if  the  value  of  </!  '</ logr  occurring  at 
the  mode  radius,  as  demonstrated  for  the  accumulation 
mode  in  Fig.  3-1  The  values  i  *•  1-3  refer  to  the  accu¬ 
mulation,  middle,  and  coarse  modes,  respectively.  The 
quantity  if  may  also  be  expressed  as  W-  l'/<r(2ir)*,  where 
I  is  the  total  particle  volume  of  the  mode.  To  more  easily 
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visualize  the  character  of  the  entire  volume  distribution 
the  value  M  was  used  instead  of  the  equivalent  expression 
containing  V.  The  least-squares  fit  of  Eq.  (3.1)  to  the 
volume  distribution  was  performed  using  the  algorithm  of 
Marquardt  (Bevington,  1969). 

The  result  of  fitting  Eq.  (3.1)  to  a  volume  distribu¬ 
tion  containing  an  accumulation  mode  was  shown  in  Fig. 
3-1.  To  obtain  a  good  fit  for  this  example,  only  the  con¬ 
stant  and  one  log-normal  curve  were  needed.  A  bimodal 
distribution  is  shown  in  Fig.  3-2a  for  data  collected  at 
Soesterberg  in  August  at  an  altitude  of  152m  above  the 
ground.  The  data  is  well  fit  using  the  constant  of 
Eq.  (3.1)  and  two  log-normal  curves,  one  to  fit  the  accu¬ 
mulation  mode  and  the  other  to  fit  the  coarse  particle 
mode.  Cases  were  found  not  only  of  the  accumulation 
mode  existing  solo,  as  shown  in  Fig.  3-1,  but  also  of  the 
coarse  particle  mode  as  the  only  mode  in  the  measured 
volume  distribution.  Figure  3-2b  shows  a  coarse  particle 
mode  measured  at  Birkhof  in  February  at  an  altitude  of 
1219m.  The  data  are  well  fit  in  this  case  with  only  one 
log-normal  distribution  and  the  constant.  Two  measured 
aerosol  distributions  showing  the  middle  mode  are  dep¬ 
icted  in  Fig.  3-3  for  Yeovil  in  March  at  152m,  Fig.  3-3a, 
and  for  Rodby  in  March  at  152m,  Fig  3-3b.  In  addition 
to  existing  solo  or  with  an  accumulation  mode,  as  shown 
in  Fig.  3-3,  the  middle  mode  was  found  in  the  presence  of 
the  accumulation  and  coarse  particle  modes  together  (not 
shown). 

The  fit  of  a  log-normal  distribution  to  the  accumula¬ 
tion  modes  presented  in  Figs.  3-1,  3-2a,  and  3-3b  shows 
that  Eq.  (3.1)  consistently  overestimates  the  data  in  the 
first  eight  to  ten  Royco  channels.  Remember  that  only 
every  other  data  channel  is  shown  in  the  figures.  This 
over  estimation  may  be  the  result  of  a  modification  to  the 
aerosol  distribution  in  the  Royco  plumbing  or  a  departure 
of  the  actual  aerosol  distribution  from  a  log-normal  distri¬ 
bution. 

To  judge  how  well  Eq  (3.1)  fit  the  entire  data  set. 
the  value  chi-square  was  calculated,  but  it  is  only  a  relative 
indication  of  how  good  the  fit  is  because  the  uncertainty  of 
the  data  is  unknown.  For  purposes  of  comparison,  the 
values  of  chi-square  for  Figs  3-3a  and  3-3b.  which  are 
almost  equal,  were  exceeded  by  20  to  25%  of  the  fitted 
data  and  for  twice  the  values  of  chi-square  by  about  10% 
of  the  data.  A  majority  of  the  large  values  of  chi-square 
are  associated  with  the  presence  of  a  coarse  particle  mode 
Even  though  a  coarse  particle  mode  appears  large  com¬ 
pared  to  one  of  the  other  volume  modes  it  can  be  heavily 
influenced  by  random  errors  in  counting  since  it  may 
represent  a  much  smaller  number  of  counted  particles  per 
channel.  For  example,  a  value  of  20 ^m’/cm'  represents 
about  3000  times  the  number  of  counted  particles  for  the 
Royco  channel  with  a  mean  radius  of  0.3nm  than  for  the 
channel  with  a  mean  radius  of  4.0 nm. 

In  the  next  section  the  characteristics  of  the  three 
modes  and  the  constant  will  be  discussed  individually 
After  which,  the  frequency  of  occurrence  of  the  modes  in 
different  combinations  as  a  function  of  altitude,  season, 
relative  humidity,  and  site  of  measurement  will  be 
presented 


Hit.  3-3.  The  volume  distribution  sampled  at  Yeovil  <a»  and  at 
Rodbv  (h)  during  March  at  an  altitude  of  152m  The  lit  of 
N  ( 1 1  to  the  dai.i  t  4- 1  is  shown  for  comparison 


3.5  Results  of  the  Fit 

Three  parameters  (M.  r,  and  .r)  are  required  to 
define  each  of  the  three  log-normal  curves  used  to  fit  the 
volume  modes.  The  behavior  of  these  three  parameters 
and  the  constant,  defined  in  Eq.  (3.1).  is  of  interest  in 
ascertaining  the  predictability  of  a  mode’s  characteristics. 
The  accumulation  mode  is  heavily  influenced  by  anthropo¬ 
genic  sources  of  aerosol,  and  dominates  scattering  of  radi¬ 
ation  in  the  visible  part  of  the  spectrum. 

Accumulation  Mode 

A  radiative  transfer  model  used  to  calculate  the  dis¬ 
tribution  of  light  in  the  atmosphere  must  in  some  manner 
approximate  the  behavior  of  this  mode  as  a  function  of 
altitude.  The  characteristics  of  this  mode,  for  the  aerosol 
samples  collected  at  the  many  different  sites  and  under 
differing  meteorological  conditions,  were  compared  by 
plotting  these  characteristics  against  the  ratio  of  the  meas¬ 
urement  altitude  to  the  altitude  of  the  inversion  base  (alti¬ 
tude  ratio).  The  height  of  the  inversion  was  determined 
by  two  methods.  The  first  method  used  the  temperature 
profile  of  the  atmosphere  measured  using  the  airborne 
temperature  transducer  during  ascents  and  descents  of  the 
aircraft  The  second  was  based  on  estimates  of  the  top  of 
the  surface  haze  layer  made  by  an  on-board  observer  A 
comparison  of  the  inversion  height  estimated  from  both 
these  methods,  shown  in  Fig  3-4,  shows  that  the  two 
methods  are  in  general  agreement  with  one  another.  Data 
departures  from  the  expected  straight  line  agreement  may 
be  in  part  due  to  differences  of  two  to  three  hours 
between  the  times  the  on-board  observer  estimated  the 
inversion  height  and  the  aircraft  flew  a  vertical  profile. 

The  behavior  of  the  accumulation  mode  with  alti¬ 
tude  is  illustrated  in  Fig.  3-5  which  shows  values  of  the  fit 
parameter  A/,  versus  altitude  ratio  The  figure  shows  that 
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Fit  3-4.  A  comparison  between  the  thickness  of  the  surface  ha/e 
layer  as  estimated  by  the  on-board  observer  and  the  height  of  the 
inversion  base  as  determined  from  the  measured  temperature 
profile 

for  this  data  set  the  accumulation  mode  rarely  exists  above 
the  inversion  layer,  an  altitude  ratio  of  1.0  Of  the  nearly 
200  occurrences  of  this  mode  shown  in  the  figure,  only  9 
clearly  exist  above  the  mixing  layer  To  eliminate  modes 
artificially  created  by  the  fitting  routine,  cases  where  Vf, 
was  less  than  0.5pm’  im1  and  r,  was  less  than  0  25 >im  -  the 
fifth  Rovco  channel  -  were  ignored. 

Often  times  in  the  literature,  the  accumulation 
mode  is  presented  in  terms  of  its  total  particle  volume 
concentration  t  which  for  a  log-normal  distribution  is 
equal  to  If,  ir,(2tr>n5.  A  knowledge  of  <r,  as  a  function  of 
Vf,  would  simplify  determinations  of  l  and  would  reduce 
the  number  of  independent  variables  in  Eq.  (3.1).  The 
values  of  <r,  are  compared  to  the  corresponding  value  of 
Vf,  in  Fig.  3-6.  The  fairly  narrow  distribution  of  the  data 
points  about  the  average  value  of  .r,  for  the  data  of  the 
figure  indicates  that  .r,  may  be  a  constant,  independent  of 
the  value  Vf,.  The  average  value  of  .r,  is  Oil.  It  is  possi¬ 
ble.  though,  that  random  error  has  masked  a  relationship 
between  vf,  and  >r,.  The  five  data  points  existing  in  the 
region  of  Vf,  greater  than  45^m'  on'  and  .r,  -0  125 
represent  data  collected  at  Rodby  on  October  25  and  26  in 
a  very  heavy  ha/e  near  the  surface  of  the  Baltic  Sea  The 
Cambridge  dewpoint  hygrometer  measured  values  of  rela¬ 
tive  humidity  of  over  90%  when  these  samples  were  col¬ 
lected  These  five  points  might  be  representative  of  condi¬ 
tions  in  a  heavy  ha/e. 

A  comparison  between  the  values  of  If,  and  r,  is 
shown  in  Fig.  3-7.  The  data,  which  shows  increasing 
values  of  Vf,  with  increasing  values  of  r,.  are  fit  remark¬ 
ably  well  by  a  straight  line  as  shown  in  the  figure  The 
regression  line  (solid  line)  was  fit  to  (he  data  using  the 
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Fig-  3-5.  The  alinude  ratio,  i  t'  ralio  of  Ihe  aircraft  altitude  dur¬ 
ing  measurement  to  the  thickness  of  Ihe  mixing  layer,  as  a  func¬ 
tion  of  the  maximum  volume  concentration  in  the  accumulation 
mode.  Vf , 
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Fig.  3-6.  The  maximum  volume  concentration  Vf,  as  a  function 
of  (he  standard  deviation  <r ,  for  Ihe  accumulation  mode. 


least  squares  criterion  and  explains  89%  of  the  observed 
variation  in  the  values  of  r,  from  the  values  of  .Vf,.  The 
linearity  of  the  data  distribution  is  a  result  of  opposing 
trends  in  f,’  and  itsut  logf,  upon  which  Vf,  depends.  The 
dashed  line  in  Fig.  3-7  is  a  curve  of  the  values 
.Vf ,  —4/3  w  ?,5  70.  where  itsut  logr  has  been  fixed  at  70  par¬ 
ticles  per  cm'.  In  comparison,  the  trend  in  the  values  of 
its i<l  log  r  are  represented  in  Fig.  3-8.  which  shows  the 
measured  values  .Vf,/?,’  averaged  over  increments  of 
0  025Mm  in  radius  The  solid  line,  the  result  of  a  hand  fit 
to  the  averages,  shows  the  trend  in  dNul  logr.  It  appears 
that  the  opposing  trends  -  sign  of  the  second  derivative  - 
of  the  r,J  and  the  Vf,/r,'  dependence  shown  in  Figs.  3-7 
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Fig.  3-7.  The  maximum  volume  conceniraiion  Mt  as  a  function 
of  ihe  mode  radius  7,  for  (he  accumulation  mode  A  straight  line 
fit  to  the  data  (solid  line)  and  the  73  dependence  of  M  (dashed 
line)  are  shown  for  comparison 
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Fig.  3-8.  The  relative  change  in  number  density  at  the  mode 
radius  of  the  accumulation  mode  as  expressed  b>  the  ratio  \1  \! r 


and  3-8,  respectively,  is  responsible  for  the  linearity  of  the 
relationship  between  ,W,  and  f ,  for  this  data  set 

There  is  no  tendency  for  the  relative  humidity  to  be 
high  for  the  larger  values  of  f.  This  dependency  was 
explored  in  report  I  in  which  it  was  found  that  any  depen¬ 
dence  on  humidity  was  mainly  due  to  the  influence  of  the 
Rodby  data  collected  in  a  heavy  haze  When  the  entire 
data  set.  as  shown  in  Fig  3-7,  is  used  no  tendency  is 
found  for  the  relative  humidity  to  be  greater  than  or  equal 
to  70%  when  the  value  of  f,  is  greater  than  or  equal  to 
0.34 in*.  The  probability  is  almost  0.2  that  the  value  of 
chi-square  will  be  equaled  or  exceeded  by  a  random  sam¬ 
ple.  The  lack  of  an  apparent  dependence  of  f,  on  relative 


Fig.  3-9.  The  altitude  ratio  as  a  function  of  the  maximum 
volume  concentration  M2  for  the  middle  mode 


humidity  may  be  due  to  a  possible  drying  of  the  aerosol 
sample  by  this  particular  measuring  system. 

Middle  Mode 

As  shown  in  Fig.  3-3  there  is  a  mode  that  exists 
between  the  regions  generally  associated  with  the  accumu¬ 
lation  and  coarse  particle  modes.  A  plot  of  the  values  of 
the  fit  parameter  for  the  middle  mode  as  a  function  of 
altitude  ratio  is  shown  in  Fig.  3-9.  The  data  of  the  figure 
show  that  the  middle  mode  exists  only  within  the  mixing 
layer  and.  most  often,  in  the  lower  half  of  the  layer  The 
values  M}  for  the  middle  mode  never  get  as  large  as  the 
values  Mt  for  the  accumulation  mode  as  seen  by  compar¬ 
ing  Figs.  3-9  and  3-5.  To  reduce  the  influence  of  errors 
associated  with  low  particle  counts,  only  values  of 
Mi  ^  15p m'/tm3  are  presented.  False  middle  modes 
caused  by  deviations  between  the  log-normal  fit  to  either 
an  accumulation  or  a  coarse  particle  mode  were  reduced 
by  ignoring  cases  where  the  value  M:  was  less  than  20%  of 
the  value  M  for  the  accumulation  or  coarse  modes  when 
present. 

The  standard  deviation  for  the  middle  mode  >r2  is 
compared  to  M2  in  Fig.  3-10  for  all  of  the  cases  presented 
in  Fig.  3-9.  The  data  are  broadly  distributed  around  the 
average  value  of  <r2,  which  is  0.29  for  the  data  shown  in 
the  figure.  The  variation  in  the  data  is  not  associated  with 
any  one  set  of  aerosol  samples  but  appears  to  be  the  result 
of  random  errors.  The  discrepancy  in  the  data  can  be 
approximated  by  examining  the  differences  in  the  values 
it 2  for  two  middle  modes  measured  within  6  minutes  (0.1 
hrs.)  of  each  other  The  average  value  of  the  difference  in 
<T  is  0.08,  which  represents  errors  in  particle  counting, 
nonhomogenities  in  the  aerosol  distribution  at  a  given  alti¬ 
tude  and  errors  in  fitting  Eq.  (3.1)  to  the  data. 

The  values  of  f2  for  the  middle  mode  are  presented 
in  Fig  3-1 1  along  with  the  corresponding  values  of  as 
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Fig.  3- 1 II.  The  nmimum  volume  concern  ration  \/:  as  a  function  '-12.  Ihe  alliluUc  r.itio  .is  a  function  of  ihc  maximum 

ul  the  vlanvl.ircl  deviation  .r.  fur  ihc  nluMIc  mode  volume  concentration  l/j  lor  ihc  coarse  movie 


PARTICLE  RADII'S.  r:(Mm> 

Fig.  .1-11.  Ihc  maximum  volume  vniKcniralmn  1/.  as  ,i  limelion 
ol  ihc  movie  radius  /_.  lor  the  middle  movie  Ihc  /  dependence 
of  M  lx  shown  lor  comparison  Mashed  line  • 

determined  from  the  til  Shown  in  the  figure  fm  com¬ 
parison  is  ihe  i  dependence  of  U  for  constant  particle 
number  (dashed  line)  The  figure  shows  lhal  the  middle 
mode  most  frequently  occurs  in  the  si/e  range  0  5  to  l„m. 
The  data  points  existing  at  /  -  Ipm  were  not  collected  at 
one  particular  data  site  and  probably  represent  the  ambient 
aerosol,  even  though  they  are  based  on  relatively  small 
particle  counts  A  ./)  ./log/  value  of  5p V  , mx  at 
/  -  (Kvpm  represents  about  seven  lime'  the  number  of 
particles  than  at  /  I  In >«  When  the  data  of  the  figure 
are  plotted  as  If.  a  general  trend  of  decreasing  particle 
number  is  seen 

Coarse  Particle  Mode 

The  third  mode  appearing  in  the  volume  distribution 


plots  of  the  data  is  the  coarse  particle  mode  Its  source  is 
generally  considered  to  be  the  result  of  mechanical 
processes,  such  as  wind  blown  dust  A  representation  of 
its  behavior  with  altitude  is  demonstrated  in  Fig  .1-12 
which  shows  values  of  for  the  coarse  mode  as  a  func¬ 

tion  of  altitude  ratio.  To  reduce  the  influence  of  errors 
associated  with  low  particle  counts,  only  values  of 
U,  ii  have  been  displayed.  The  figure  clearly 

demonstrates  that  the  coarse  mode  is  not  confined  to  the 
mixing  layer  as  were  the  accumulation  and  middle  modes 
On  numerous  occasions  the  on-board  observer  reported 
'he  existence  of  ha/e  layers  above  the  surface  ha/e  layer 
Individual  examples  or  ha/e  layers  existing  at  altitude 
ratios  greater  than  1.0.  which  appear  to  contain  coarse 
mode  particles  but  no  accumulation  mode  particles,  were 
presented  in  report  1  The  set  of  seven  data  points  in  Fig 
1-12  centered  around  the  altitude  ratio  8  are  from  meas¬ 
urements  taken  at  Rodhy  and  Birkhof  The  Rodby  data 
were  collected  in  the  altitude  range  5000  to  6000m  above 
the  sea  surface,  while  the  Birkho1'  data  were  collected  at  an 
altitude  of  about  .1000m  above  the  ground  It  is  possible 
that  this  high  altitude  aerosol  was  transported  from  the 
south  since  a  majority  of  the  occurrences  of  the  coarse 
particle  mode  occurred  at  Trapani,  as  will  be  shown  in 
Table  .1  2  in  the  following  section  A  trajectory  analysis  of 
the  air  mass  containing  the  coarse  mode  particles  would  be 
necessary  to  confirm  this 

The  fit  parameter  .r ,  is  plotted  in  Fig.  5-13  as  a 
function  of  W<  for  all  of  the  data  presented  in  Fig.  3-12 
The  large  spread  in  the  data,  which  may  well  be  masking  a 
dependence  between  If,  and  ■»,.  is  probably  mostly  due  to 
noise.  The  average  deviation  in  .r,  is  0  09  for  data  sam¬ 
ples  collected  within  6  minutes  of  one  another.  The  data 
arc  fairly  well  distributed  about  the  value  0.21.  the  average 
value  of  .r,  The  data  points  near  the  upper  right  hand 
corner  of  the  figure,  with  <r  ,  in  the  range  0.3  to  0  4pm,  do 
not  represent  data  from  one  particular  station  or  month 
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Fig.  3*13.  The  maximum  volume  concentration  \/*  as  a  function 
of  the  standard  deviation  for  the  coarse  mode 


Fig.  3-15.  1  he  altitude  r.iuo  a>  a  function  of  i he  con-aunt  k  of 
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Fir.  5-14.  The  maximum  voli-mc  concentration  1/,  as  a  function 
of  (he  mode  radius  r ,  lor  ihc  coarse  mode 


The  values  of  r,  for  (he  coarse  panicle  mode  are 
shown  in  Fig  3-14  wnh  the  corresponding  values  of  U, 
The  distribution  of  data  points  in  ihc  figure  is  noticeably 
asymmetric  around  the  average  value  of  the  mode  radius 
for  the  data,  r,  --  5 ^m.  The  lack  of  any  data  points  in  the 
region  Vf,  >  5m m'  an'  and  r  •  y,,m  is  not  the  result  of  a 
deliberate  filtering  of  the  data  It  should  be  noted  that  the 
last  channel  of  the  fitted  Rovco  data  had  a  mean  radius  of 
5  9 nm.  which  is  probably  responsible  for  the  lack  of  sym¬ 
metry  in  the  data.  Modes  with  a  value  of  r,  in  the  region 
of  5.9 pm  or  greater  may,  therefore,  be  an  artifact  of  the 
model  fit  which  assumes  a  log-normal  shape  for  the  coarse 
mode.  A  comparison  between  the  values  of  r  for  the  mid¬ 
dle  and  coarse  particle  modes,  as  shown  in  Figs  3-11  and 
3-14.  respectively,  shows  a  distinct  break  between  these 


modes  at  a  particle  tadius  centered  around  2 ><m. 

The  Constant 

In  addition  to  the  three  log-normal  distributions  that 
were  used  to  describe  the  volume  modes,  a  constant  k  was 
included  in  Eq.  (3.1)  The  need  for  the  constant  is 
demonstrated  by  the  zero  offset  of  the  example  aerosol 
distributions  presented  in  Figs.  3-1  to  3-3  and  by  (he  many 
cases  where  there  is  no  discernible  mode  present  in  the 
volume  distribution.  The  value  of  the  constant  as  a  func¬ 
tion  of  altitude  ratio  is  shown  in  Fig.  3-15  for  all  of  the 
fitted  data.  As  seen  in  the  figure,  the  value  of  k  appears 
to  decrease  with  increasing  altitude  ratio,  reaching  a  max¬ 
imum  of  about  4m in'’ an'  near  the  ground  and  close  to 
0.2 um’  «'  at  the  highest  altitude  where  measurements 
were  made. 

Values  of  k  averaged  over  atmospheric  layers  of 
increasing  thickness  with  increasing  height  are  shown  in 
Fig  3-16  plotted  against  the  mean  value  of  altitude  ratio 
for  the  slab  layer  (circled  points).  At  altitude  ratios  of 
greater  than  I  the  average  values  of  the  constant  are  fit 
well  by  an  equation  of  the  form 


k  =  0  12  *  0.82(1  :i  (3.2) 


where  k  is  the  average  and  is  altitude  ratio  This  equa¬ 
tion  explains  96%  of  the  variation  in  k  A  straight  line 
was  fit  to  the  three  values  of  k  below  .-  =  1 

The  value  of  k  represents  an  estimation  of  the  resi¬ 
dual  total  particle  volume  of  the  aerosol  distribution  after 
the  log-normal  approximations  to  the  modes,  if  present, 
are  subtracted  out  No  clear  relationship  between  k  and 
the  values  of  \f  lor  the  three  modes  could  be  found  It 
appears  that  there  is  a  trend  of  increasing  k  with  increas¬ 
ing  total  particle  volume  of  the  three  modes,  but  the 
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Fig.  3-16.  The  altitude  ratio  .is  a  (unction  ol  k  averaged  over 
atmospheric  lavers  ol  increasing  thickness  with  increasing  height 
The  solid  line  represents  a  functional  lit  to  the  average  values  ol 
k 


scatter  in  the  data  is  great  Depending  nn  how  accurate 
the  aerosol  data  need  to  he  fit.  it  may  be  possible  to  drop 
the  constant  from  Fq.  <3.11  when  a  volume  mode  is 
present  For  example,  report  I  showed  that  the  measured 
volume  scattering  coefficient  at  the  wavelength  0  66/ou 
was  well  correlated  with  the  value  W,  of  the  accumulation 
mode,  independent  of  the  value  k. 

3.6  Combinations  of  the  Modes 

In  the  previous  section,  the  behavior  of  the  three 
modes  and  the  constant  with  altitude  was  studied,  along 
with  the  dependence  of  r  and  >r  on  U  for  each  of  the 
modes  It  is  of  interest,  now.  to  explore  how  frequently  in 
the  data  base  various  combinations  of  the  modes  occur 
and  if  the  frequency  of  occurren,  >1  the  combinations 
shows  a  preference  toward  altitude,  season,  relative  humi¬ 
dity.  or  site  of  measurement. 

The  first  two  columns  of  Table  3.1  present  the 
number  of  occurrences  of  the  different  combinations  of 
the  modes  within  and  above  (he  mixing  layer.  Each  of  the 
possible  combinations  are  labeled  using  a.  m.  or  <  to 
denote  the  accumulation,  middle,  or  coarse  particle 
modes,  respectively  The  label  none  refers  to  the  absence 
of  a  discernible  mode  The  volume  distributions  listed 
under  none  were  fit  with  only  the  constant  of  F.q  (3  1). 
which,  for  this  case  only,  is  equivalent  to  a  fit  of  the  aero¬ 
sol  number  distribution  using  a  Junge  distribution  with 
d".V  The  figures  in  parenthesis  give  the  distribution 
which  would  be  expected  if  there  were  no  tendency  for  the 
different  mode  combinations  to  exist  either  at  .  -  I  or 
:  >  I  The  table  shows  that  combinations  containing  the 
accumulation  or  middle  modes  rarely  exist  above  the  mix¬ 
ing  layer,  as  seen  in  Figs  3-5  and  3-d  In  contrast,  the 
coarse  mode  solo,  row  i.  is  equally  Itkelv  to  occur  above  as 
within  the  mixing  layer  as  seen  b\  the  small  differences 


Table  3.1.  The  occurence  of  different  combinations  of  the  modes 
as  a  function  of  altitude  ratio,  season,  and  relative  humidity. 


:>l 

Summer 

Winter 

RH<  70 

RH>70 

a 

85  (53) 

3  (35) 

58  (39) 

30  (49) 

35  (39) 

50  (46) 

a-m 

15  (9) 

0  (6) 

1  (7) 

>4  (8) 

4  (7) 

11  (8) 

a-c 

53  (34) 

4  (23) 

46  (25) 

II  (32) 

33  (23) 

18  (28) 

m 

19  (11) 

0  (8) 

4  (8) 

15  (II) 

4  (9) 

15  (10) 

m-c 

5  (3) 

0  (2) 

0  (2) 

5  (3) 

3  (2) 

2  (3) 

a-m-c 

18  HD 

0  (7) 

2  (8) 

16  (10) 

3  (8) 

15  (10) 

l' 

43  (48) 

r?  <32 > 

35  (35) 

45  (45) 

29  (17) 

8  (20) 

none 

88  11511 

175  (106) 

94  ( 1 16) 

169  (147) 

30  (35) 

47  (42) 

between  the  measured  number  of  occurrences  above  and 
within  the  mixing  layer  and  those  expected  from  chance. 
Samples  taken  within  10  minutes  of  one  another  at  the 
same  altitude  were  counted  as  one  sample  to  help  reduce 
the  bias  introduced  by  multiply  sampling  the  same 
ambient  aerosol  In  addition  to  showing  where  the 
different  mode  combinations  are  likely  to  occur  in  the 
atmosphere,  the  first  two  columns  of  Table  3.1  indicate 
what  combinations  are  least  prevalent.  For  :  >  1,  it  is 
likely  that  no  mode  will  be  observed  in  a  volume  distribu¬ 
tion  plot  as  seen  from  (he  large  number  of  occurrences  of 
none  compared  to  the  other  possible  combinations.  The 
combination  of  a  middle  mode  with  another  mode  at  :  $  1 
will  probably  be  measured  infrequently  according  to  this 
data  base.  It  is  more  likely  that  at  $  I  an  accumulation 
mode  will  be  measured  solo  or  in  combination  with  a 
coarse  mode. 

The  third  and  fourth  columns  of  Table  3.1  indicate 
the  seasonal  preference  of  the  combinations.  The  summer 
season  was  defined  for  this  investigation  to  include  the 
months  April  to  August,  the  other  months  of  the  year 
defining  the  winter  season.  The  table  shows  that  the  mid¬ 
dle  mode,  either  solo  or  in  a  combination,  tends  to  occur 
in  the  winter  months.  The  coarse  mode  by  itself  shows  no 
preference  toward  season,  while  the  combinations  a  and 
ii -i  are  likely  to  occur  in  summer.  Combinations  including 
a  and  m.  which  show  a  tendency  to  occur  in  the  summer 
and  winter  seasons,  respectively,  show  the  largest  values 
of  3/  in  those  respective  seasons.  The  combination  none 
indicates  a  preference  toward  winter. 

The  last  two  columns  of  Table  3  I  indicate  the 
preference  of  the  various  combinations  for  high  or  low 
humidity  The  data  presented  in  these  columns  represent 
aerosol  samples  collected  within  the  mixing  layer  only, 
ihus  avoiding  freezing  temperatures  that  can  occur  at 
'  'gher  altitudes  The  total  number  of  measured  samples 
in  the  last  two  columns  may  i.e  less  than  the  number  of 
samples  in  the  first  cc  <nr  a  given  combination,  since 
for  some  samples  no  of  relative  humidity  is  avail¬ 
able  As  shown  in  the  ..  ■  the  combinations  m,  a-m. 

and  a-m- c  which  contain  the  middle  mode  show  a  strong 
tendency  to  exist  at  the  higher  values  of  relative  humidity 
The  coarse  mode,  in  contrast,  shows  a  strong  tendency  to 
exist  at  the  lower  values  of  relative  humidity  when  in 
combinations  of  a- c  and  c.  The  apparent  preference  of  a, 
a-m.  and  none  toward  the  higher  values  of  humidity  is  not 
large  enough  to  be  statistically  significant 
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Table  3.2.  The  Occurrence  of  Different  Combinations 
of  the  Modes  as  a  Function  of  Site 
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The  number  of  times  the  different  mode  combina¬ 
tions  appear  in  the  volume  distributions  lor  the  different 
sites  of  measurement  is  shown  in  Table  3.2  The  percent 
of  the  total  number  of  flights  each  combination  represents 
within  and  above  the  mixing  layer,  separately,  is  shown 
for  each  site  The  two  most  southeastern  sites,  Trapani 
and  Birkhof.  show  the  highest  percent  of  occurrence  of 
the  coarse  particle  mode  solo,  row  c.  both  within  and 
above  the  mixing  layer.  Front  examination  of  synoptic 
charts  at  the  time  of  these  measurements  it  is  possible  that 
the  coarse  mode  material  had  its  origin  in  Northern 
Africa  It  is  interesting  that  the  solo  coarse  mode  exisis 
both  above  and  within  the  mixing  layer  for  the  two  most 
southeastern  sites,  but  its  existence  for  the  other  stations 
is  mainly  confined  to  altitudes  above  the  inversion  height 
Table  3  2  shows  that  the  accumulation  mode,  row  a. 
occurs  more  frequently  for  the  northern  most  stations 
The  middle  mode  in  combination  with  other  modes  and 
solo  appears  to  predominate  at  Rodby.  Meppen.  booster 
berg,  and  Yeovil  The  percent  of  occurrence  m  which 
there  was  no  discernible  volume  mode  is  large  at  altitudes 
above  the  mixing  layer,  as  seen  in  ihe  table  Within  the 
mixing  layer,  '.he  southern  most  stations  Irapam.  Birkhol 
and  Bruz  show  the  greater  occurrence  of  the  combination 
none.  The  number  of  measured  aerosol  samples  from  Mil  - 
denhall  are  too  small  to  see  a  pattern  for  this  site 

The  data  presented  in  Tables  3  I  ami  3  2  indicate 
that  the  volume  modes  are  lied  to  the  characteristics  >  I 
the  air  mass  containing  the  aerosol  sample  There  is  a 
definite  dependence  of  the  modes  on  altitude,  season,  and 
site.  The  dependence  on  relative  humidity  was  very 
strong  for  the  combinations  i  and  u-i.  which  were  fre¬ 
quently  measured  at  values  of  low  humidity  Ihe  depen¬ 
dence  was  not  as  definite  for  the  combinations  showing  a 
preference  for  the  higher  values  of  humidity,  a  possible 
result  of  a  drying  of  the  sample  in  ihe  intake  plumbmg 

3.7  Discussion  and  Conclusions 

An  extensive  set  ol  aerosol  si/e  distributions  meas¬ 


ured  m  ihe  lower  troposphere  over  Furope  and 
(Treat  Britain  were  analyzed  using  a  volume  (mass I 
display  As  has  been  recognized  in  the  past,  departures  of 
the  aerosol  number  distribution  from  a  power-law  distribu¬ 
tion  give  rise  to  volume  modes  The  analysis  ol  the 
characteristics  and  behavior  of  the  volume  modes  supports 
the  following  conclusions 

1  The  volume  distribution  can  be  described  well  by 
F.q  (3  II.  which  uses  a  log-normal  distribution  to 
describe  the  individual  volume  modes  The  data 
show  that  there  are  three  distinct  modes  accumula¬ 
tion.  middle,  and  coarse 

2  The  accumulation  mode  almost  always  exists  within 
the  mixing  layer  Its  value  ol  movie  radius  is  m  the 
range  0  25  to  0  47,,  m  and  it  is  well  described  as  a 
linear  lunction  ol  the  maximum  volume  ol  the 
mode.  U,  The  standard  deviation  a,  is  probably  a 
constant  ol  value  0  I  I 

3  The  middle  mode  is  confined  to  the  mixing  layer 
and.  generally,  to  the  lower  hall  ol  the  layer  I  Ik 
values  ol  i.  are  in  the  range  0  45  to  I  4i,.y  Most  ol 
(he  data  follows  a  trend  ol  increasing  l/.  with 
increasing  however,  all  the  data  show  a  general 
trend  of  decreasing  U  ■  with  increasing  Ihe 
values  of  it  appeared  lo  he  centered  around  a  mean 
value  ol  0  3 

4  The  data  indicate  that  Ihe  coarse  particle  mode  can 
exist  either  ahove  or  within  the  mixing  layer  Any 
relationships  between  Us.  . r ...  and  /,  that  may  exist 
lor  this  movie  are  masked  hv  what  appears  to  be  the 
influence  ol  random  noise  The  mean  values  ol  ,r. 
and  is  for  die  luted  data  set  are  02  and  5s. m. 
i  -speetivelv 

s  In  addition  to  the  three  log-normal  distributions 
used  to  describe  the  aerosol  volume  distribution,  a 
constant  was  used  in  l  q  (3  II.  The  value  o!  the 
constant  is  generally  less  than  the  value  ol  U  lor 
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any  of  the  modes  analyzed.  The  values  of  the  con¬ 
stant  averaged  over  intervals  of  altitude  ratio  are 
found  to  be  well  described  by  an  equation  of  the 
form  k  -  a+ft/c  where  a  and  ft  are  constants.  ;  is 
altitude  ratio,  and  k  is  the  average  value  of  the  con¬ 
stant. 

6  Within  the  mixing  layer  the  most  common  combina¬ 
tions  of  the  modes  are:  no  mode  at  all.  the  accumu¬ 
lation  mode  solo,  the  accumulation  and  coarse  mode 
together  (bimodal  distribution),  and  the  coarse 
mode  solo.  There  is  a  small  but  significant  number 
of  occurrences  of  the  combinations:  middle  mode 
solo:  accumulation,  middle,  and  coarse  modes 
together,  and  accumulation  and  middle  modes 
together.  The  combination  of  the  middle  and  coarse 
modes  rarely  occurs.  Above  the  mixing  layer  the 
combination  no  mode  predominates. 

7  The  middle  mode  lends  to  be  measured  in  the 
winter  season,  while  the  accumulation  mode  solo 
appeared  more  often  in  the  summer  data  than 
expeaed  from  random  chance.  The  occurrence  of 
the  coarse  particle  mode  solo  showed  no  preference 
toward  season. 

8  The  middle  mode  appears  more  often  when  the  rela¬ 
tive  humidity  is  high,  but  the  coarse  particle  mode 
shows  a  preference  toward  (he  low  values  of  humi¬ 
dity.  The  accumulation  mode  shows  no  significant 
dependence  on  relative  humidity.  The  relative 
humidity  dependence  may  have  been  affected  by  a 
drying  of  the  aerosol  in  the  intake  plumbing. 

9.  The  northern  most  stations.  Rodbv.  Meppen.  and 
Soesterberg.  show  a  higher  occurrence  of  the  accu¬ 
mulation  mode  than  the  southern  sites  The  south¬ 
ern  sites,  in  contrast,  show  a  larger  number  of 
measured  coarse  particle  modes  than  the  northern 
stations  The  middle  mode  seems  to  occur  most 
frequently  at  Rodby,  Meppen.  Soesterberg.  and 
Yeovil. 

There  have  been  a  number  of  aerosol  models 
developed  over  the  last  two  decades  that  have  been  used 
to  describe  the  aerosol  si/e  distribution,  as  mentioned  in 
the  introduction  Several  of  these  models,  presented  in  a 
volume  display,  are  shown  in  Fig.  3-1 7.  The  curves 
labeled  I  to  4  of  Fig  3- 17a  represent  the  aerosol  distribu¬ 
tion  predicted  from  the  Munn-Kai/  maritime  model  The 
model  consists  of  a  background  continental  component 
which  is  a  function  of  relative  humidity  only,  and  a  mari¬ 
time  contribution  which  is  considered  to  be  a  function  of 
relative  humidity,  wind  speed,  and  altitude  For  the  four 
curves  the  respective  values  of  relative  humidity,  wind 
speed,  and  altitude  are  1 )  70%.  6m/s.  150m.  2)  80%. 
4m/s,  150m.  3)  70%,  4m/s.  150m.  and  4)  70%.  4m/s. 
1200m.  The  continental  component  appears  as  a  straight 
horizontal  line  in  the  volume  display  While  (he  model 
does  depend  on  several  scaling  factors,  the  four  curves 
show  that  the  effects  of  the  maritime  component  predom¬ 
inate  at  particle  radii  greater  than  I  pm.  The  rural  and 
maritime  models  of  Shettle  and  Fcnn  (1979).  who  use  a 
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Fin-  3-17.  Volume  disinbution  plots  based  on  the  Munn-Kat/ 
maritime  aerosol,  curves  I  to  4  (Hughes  and  Richter.  1980).  a 
rural  and  maritime  model,  curves  5  and  b.  respectively  (Shelile 
and  Fenn.  1979).  and  the  Ha/e  M  and  lla/e  /.  models,  curves  7 
and  8.  respectively  (l)eirmendjian.  1969) 


log-normal  distribution,  are  shown  by  curves  5  and  6. 
respectively,  of  Fig.  3- 1 7b  for  a  relative  humidity  of  70% 
The  shape  of  curve  5  implies  a  bimodal  distribution.  The 
model  standard  deviations  are  constants  of  value  0.35  for 
the  maritime  and  the  subnncrometer  rural  modes  and  0.4 
for  the  rural  coarse  particle  mode  As  the  relative  humi¬ 
dity  changes  from  0  to  99%.  the  mode  radii  vary  from  0.19 
to  0.37,  5.5  to  15.  and  1.1  to  5.3 pm  for  the  two  modes  of 
the  rural  model  and  the  single  mode  of  the  maritime 
model,  respectively.  The  curves  labeled  7  and  8  represent 
the  Haze  A/  (marine  or  coastal)  and  Haze  L  (continental) 
models,  respectively,  of  Deirmendjian  (1969).  as  taken 
from  his  Table  5.  The  total  particle  volume  of  the  Haze 
M  model  was  adjusted  for  the  purpose  of  display.  The 
difference  between  'he  mode  radii  of  the  similar  curves  6 
and  7  and  the  curves  of  the  Munn-Kalz  model  might 
represent  differences  in  the  formation  history  of  the  aero¬ 
sols  stressed  by  these  maritime  models.  The  Munn-Katz 
model  may  stress  a  pure  wind  generated  aerosol  while 
curves  6  and  7  may  represent  a  more  mature  maritime 
aerosol  with  some  continental  influence. 

It  is  interesting  that  the  models  presented  in  Fig. 
3-17  not  only  give  rise  to  volume  modes  of  a  similar  shape 
to  those  seen  in  the  data  but  suggest  the  presence  of  three 
separate  aerosol  regimes  The  submicrometer  continental 
aerosol  modes  of  curves  5  and  8  form  the  first  region  with 
mode  radii  in  the  range  0.2  to  about  0.6p m,  the  second 
region  is  composed  of  the  maritime  modes  of  curves  1  to 
4.  6.  and  7  with  mode  radii  in  the  range  from  about  1  to 
3pm.  and  the  large  particle  mode  of  curve  5  forms  the 
third  region  in  the  range  5  to  15pm.  A  similar  division  is 
implied  by  the  behavior  and  characteristics  of  the  volume 
modes  of  the  fitted  data  set.  though  the  range  of  the 
regimes  is  slightly  different.  The  measured  accumulation 
mode  exists  from  the  lower  limit  of  the  Royco  to  about 
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0.47Mm,  the  middle  mode  exists  from  0.45  to  1.4m»i.  and 
the  coarse  particle  mode  is  in  the  region  3m m  to  the 
Royco’s  upper  limit,  5  9 nm.  Such  a  division  also  appears 
in  the  aerosol  volume  distributions  presented  in  Patterson 
and  Gillette  (1977)  and  Patterson  ei  al.  (1980)  These 
investigators  describe  a  submicrometer  mode  composed  of 
anthropogenic  aerosols,  a  maritime  aerosol  mode  centered 
at  a  radius  of  about  1m">,  and  an  aerosol  mode  at  about 
3m  m  composed  of  soil  particles. 

The  values  of  mode  radius  and  standard  deviation 
for  the  different  volume  modes  have  been  summarized  in 
Table  3.3.  The  symbol  --  indicates  that  no  value  was 
given  for  this  parameter.  The  mode  radii  for  the  accumu¬ 
lation  and  coarse  particle  modes  appear  to  agree  fairly  well 
The  standard  deviations  of  these  modes  for  the  aircraft 
data  are,  however,  significantly  less  than  those  found  in 
the  data  of  Whitby  (1978)  and  Patterson  and  Gillette 
(1977)  and  the  rural  model  of  Shetlle  and  Feint  <1979). 
For  the  coarse  mode,  this  may  be  a  result  of  low  particle 
counts  and  an  incomplete  sampling  of  the  entire  mode 
which  was  often  measured  near  the  upper  limit  of  the 
Royco  counter.  Investigations  of  the  discrepancy  existing 
for  the  accumulation  mode  are  continuing.  Preliminary 
analysis  of  the  submicrometer  mode  measured  near  ihe 
ground  at  Meppen,  with  a  separate  particle  counter,  show 
standard  deviations  of  about  the  same  values  as  those 
determined  from  the  aircraft  data.  The  mean  standard 
deviation  of  the  fitted  middle  mode  approximates  the 
value  used  in  the  maritime  model  of  Shetlle  and  Fenn 
(1979).  The  volume  mean  radius  of  1  y.m  measured  by 
Patterson  ei  al.  (1980)  for  a  Pacific  marine  boundary-layer 
aerosol  is  in  the  range  of  the  mode  radii  determined  from 
the  aircraft  data  for  the  middle  mode.  Care  must  be  exer¬ 
cised  in  interpreting  the  middle  mode  of  the  aircraft  data 
as  a  marine  aerosol,  since  the  source  of  the  aerosol  is  not 
known. 

Recommendation 


The  data  suggest  that  the  accumulation  mode  is  strongly 
associated  with  anthropogenic  sources  of  aerosol,  the  mid¬ 
dle  mode  with  a  maritime  environment,  the  coarse  mode 
with  wind  blown  dust,  and  the  constant  with  an  average 
background  aerosol.  Such  a  division  of  the  aerosol  size 
range  would  be  of  aid  in  predicting  the  optical  characteris¬ 
tics  of  an  aerosol. 


4.  SUMMARY 

This  Final  Report  under  contract  F19628-78-C-0200 
has  addressed  the  reduction  and  publication  of  previously 
unprocessed  raw  data  gathered  during  the  period  1970 
through  1978  by  the  airborne  instrument  system  and 
measurement  program  summarized  in  Duntley  ei  al. 
(1978).  The  Aircraft  Data  Reduction  effort,  task  one 
under  this  two  part  research  contract,  has  concentrated  on 
the  measurements  of  atmospheric  optical  properties,  pri¬ 
marily  scattering  coefficient  profiles  plus  sky  and  terrain 
radiances,  that  were  collected  during  the  Winter  and  Sum¬ 
mer  seasons  of  1978.  During  this  interval,  there  were 
approximately  45  data  flights  conducted  throughout 
western  and  northern  Europe.  The  atmospheric  optical 
data  from  these  and  other  flights  have  been  published  in 
six  of  the  ten  OPAQUE  related  AFGL  Technical  Reports 
listed  in  Table  2.3.  Over  670  aerosol  size  distributions 
covering  particle  radii  from  0.2  to  5.9  micrometers  were 
measured  contemporaneously  with  the  measurements  of 
sky  and  terrain  radiances  during  these  flights  As  noted  in 
Table  2.3.  portions  of  these  aerosol  data  have  also  been 
previously  reported  during  this  contract  interval. 

The  basic  data  types  which  have  been  processed  and 
reported  during  this  contract  interval,  and  thus  summar¬ 
ized  and  reviewed  in  this  Final  Report,  fall  into  the  fol¬ 
lowing  general  categories:  a)  altitude  profiles  of  atmos¬ 
pheric  volume  scattering  coefficient,  b)  fixed  altitude 
measurements  of  sky  and  terrain  radiances,  and  c)  fixed 
altitude  measurements  of  aerosol  size  distributions. 


In  future  aerosol  studies  it  is  recommended  that  a 
greater  effort  be  placed  on  determining  the  history  of  the 
air  mass  containing  the  aerosol  sample  The  dependencies 
on  altitude,  season,  and  site  shown  in  the  two  tables  for 
the  different  combinations  of  the  modes  could  be  better 
explained  if  the  history  of  the  aerosol  sample  were  known. 


Table  3.3.  Summary  of  (he  Volume  Mode  Parameters 


Accumulation 


Fitch  and  (  revs 

Whitby  M  W7*» 
average  continental  aerosol 

Pattenon  and  Gillette  1 1977* 
tropospheric  aerosol  with  dust 

Patterson  ei  al  t|9g0) 
maritime  aerosol 

Shettle  and  Fenn  H«T9t 
rural  model 
.  man  time  model 
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Several  general  characteristics  of  the  full  6  kilometer 
profiles  of  volume  scattering  coefficients  can  be  summar¬ 
ized  as  follows.  In  both  the  Winter  and  Summer  data 
sets,  the  Trapani  profiles  are  clearly  different  from  their 
companion  profiles  measured  in  northern  Europe.  In  the 
Winter  set.  the  Trapani  profiles  show  significantly  less  low 
altitude  structure  than  any  of  the  other  sets,  and  have  a 
complete  lack  of  any  abrupt  haze  top.  Similarly,  the  Sum¬ 
mer  Trapani  profiles  show  little  structure  and  additionally 
maintain  relatively  high  scattering  coefficient  values 
throughout  the  entire  6  km  altitude  band.  It  appears  that 
there  is  substantially  more  structural  irregularity  in  the 
scattering  coefficient  profiles  during  the  Winter  set  than 
during  the  Summer  This  is  particularly  true  for  the 
northern  European  stations.  Soesterberg.  Meppen  and 
Rodby. 

Several  illustrations  of  very  low  altitude  i.e.  ground 
level  to  700m  AGL.  photopic  volume  scattering  coefficient 
profiles  are  presented  from  the  fifty-four  profiles  previ¬ 
ously  reported  The  basic  question  addressed  by  these 
data  is  whether  or  not  the  scattering  coefficient  profile 
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remains  reasonably  eonstam  as  one  approaches  the  surface 
from  an  altitude  of  several  hundred  meters.  These  data, 
representing  all  four  temporal  seasons,  indicate  that  in 
forty-seven  out  of  fifty-four  (87%)  instances,  the  profile  is 
essentially  constant  in  value  within  this  very  low  altitude 
regime 

Sky  and  terrain  radiance  distributions  selected  from 
data  measured  during  all  five  OPAQUE  deployments  illus¬ 
trate  the  full  gamut  of  meteorological  conditions  from 
clear  and  cloud  free  to  fully  overcast.  A  remarkable 
degree  of  radiance  symmetry  and  uniformity  is  illustrated 
for  the  full  overcast  conditions.  However,  under  less 
extreme  coverages,  a  consistent  pattern  of  well  behaved 
radiance  patterns,  both  as  a  function  spectral  hand  and  of 
altitude,  is  easily  recognized  One  notes,  that  under  con¬ 
ditions  of  highly  variable  cloud  cover,  the  measured  sky 
and  terrain  radiances  tend  to  fluctuate  both  above  and 
below  the  radiance  distribution  of  a  clear  sky  by  a  factor  of 
two  or  more  in  some  instances,  but  in  general  remain 
within  about  ±50  percent  of  the  clear  sky  radiance  for 
corresponding  paths  of  sight 

Corrective  techniques  developed  during  this  contract 
interval  have  been  applied  to  selected  sub  sets  of  scatter¬ 
ing  coefficient  and  sky  radiance  measurements.  Analysis 
of  these  data,  plus  additional  measurements  of  directional 
scattering  function  at  scattering  angles  of  30  and  150 
degrees  led  to  the  conclusion  that,  in  most  cases,  the 
measurements  of  the  upper  altitude  total  volume  scatter¬ 
ing  coefficients  were  too  high  The  evidences  are  con¬ 
sistent  in  their  predictions  of  photopic  optical  scattering 
ratios  £)(:)  of  from  13  to  26  at  6  kilometers  altitude, 
roughly  half  the  measured  value  Additional  validations 
of  the  techniques  are  recommended 

Measured  near-sun  radiances  were  found  to  be  con¬ 
sistently  too  high  due  to  off-axis  stray  light  entering  the 
scanner's  optical  system  as  its  field  of  view  swept  through 
the  solar  disc.  The  implications  of  subsequent  testing 
were  that  high  altitude  or  clear  day  measurements  of  sky 
radiance  within  10°  of  the  solar  disc  will  have  stray  light 
ratios  of  as  much  as  20  to  35.  whereas  low  altitude  near¬ 
sun  measurements  through  more  turbid  atmospheres  will 
have  substantially  lower  stray  light  ratios  of  about  3  to  6  at 
the  same  angular  distance  from  the  disc  There  were  no 
detectable  stray  light  influences  upon  measurements  made 
beyond  25°  from  the  disc  under  any  of  the  conditions 
tested 

More  than  670  aerosol  si/e  distributions  covering 
particle  radii  for  02  to  59 /zm  were  collected  contem¬ 
poraneously  with  the  measurements  of  sky  and  terrain 
radiances.  The  measured  volume  (mode)  -  size  distribu¬ 
tions  were  characterized  by  three  additive  log-normal  dis¬ 
tributions,  one  for  each  of  the  three  volume  modes 
observed  in  the  data  The  behavior  of  the  modal  parame¬ 
ters  for  each  of  the  three  volume  modes  and  of  the  modes 
as  a  function  of  altitude,  humidity,  season,  and  geographic- 
site  of  measurement  were  examined  resulting  in  the  fol¬ 
lowing  conclusions 


1  The  accumulation  mode  is  almost  always  confined  to 
the  mixing  layer  The  mode  radius  r  ranges  from 
0.25  to  0.47/am  and  its  standard  deviation  cr  is  a 
constant  of  approximate  value  0.11.  The  mode 
radius  is  a  linear  function  of  the  value  dV/d  log r  af 
the  7.  The  accumulation  mode  is  the  most  pre¬ 
valent  volume  mode  occurring  as  the  only  mode 
more  often  in  the  summer  season  than  in  the  winter 
season  data.  It  shows  no  significant  dependence  on 
relative  humidity  The  accumulation  mode  appears 
more  frequently  in  the  data  collected  at  the  northern 
stations  than  the  southern  stations. 

2.  The  middle  mode  occurs  only  within  the  mixing 
layer,  most  often  in  the  lower  half  of  the  layer.  The 
mode  radius  ranges  from  0  45  to  1.4/xm  and  its 
standard  deviation  appears  to  be  a  constant  of  value 
0  3.  The  middle  mode  appears  more  often  in  the 
winter  than  in  the  summer  season  data  and  more 
often  when  the  relative  humidity  is  greater  than  70% 
than  when  it  is  less  than  70%. 

3.  The  coarse  particle  mode  can  exist  above  or  within 
the  mixing  layer.  The  mean  values  of  <r  and  r  for 
this  mode  are  0.2  and  5/am.  respectively.  The 
coarse  mode  shows  no  preference  toward  season, 
but  appears  more  often  in  the  low,  less  than  70%, 
humidity  data.  The  mode  appears  more  often  in  the 
data  collected  from  the  southern  most  sites. 

4.  At  altitudes  above  the  mixing  layer  the  volume 
display  typically  shows  no  discernible  modes  and  is 
well  fit  by  a  constant. 

5  The  data  and  comparisons  with  aerosol  models 
imply  that  the  volume  modes  are  associated  with 
specific  aerosol  sources  and.  therefore,  with  the  his¬ 
tory  of  the  airmass  in  which  the  aerosol  sample  was 
measured. 
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Abstract  The  ability  of  a  human  observer  to  perform  tasks  involving  visual  search  and  the  determination  ot 
visual  ranges  (i.e..  visibilities),  is  a  commonly  recognized  skill  that  is  presenlly  under  increasing  competition 
from  a  variety  of  sophisticated  electro-optical  systems  The  paper  reviews  a  few  of  ihe  pertinent 
psychophysical  characteristics  related  to  the  human  observer's  performance,  how  they  arc  related  to  some 
instrumental^  measured  quantities,  and  two  of  Ihe  more  commonly  encountered  sources  ol  error  in  visual 
determinations  of  visibility  Comparisons  between  several  observer  and  instrumental  determinations 
illustrate  the  link  between  the  two  experimental  regimes  and  introduce  some  appropriate  specialized  data 
obtained  with  an  airborne  nephelomcter  Since  the  determination  of  visibility  is  Ihe  simpler  special  case  ol 
determining  horizontal  path  contrast  transmittance,  several  characteristics  of  the  nephelomcter  data  which 
influence  its  use  in  the  more  general  determination  of  slant  path  contrast  transmittance  are  illustrated  In  this 
context,  the  implication  of  these  data  in  ones  ability  to  accurately  model  atmospheric  influences  upon 
contrast  transmittance  is  introduced  The  basic  conclusions  to  be  drawn  from  this  paper  are  twofold  I  irst. 
that  although  human  visual  performance  is  a  complex  psychophysical  task,  there  is  little  doubt  that  lot 
visibility  determinations  tl  can  be  adequately  related  to  instrumental  measurements  of  atmospheric  volume 
scattering  coefficient  and  second  that  a  reliable  specification  ol  the  directional  scattering  properties  ol  the 
atmospheric  aerosol  can  be  reasonably  deduced  from  this  measurement  of  total  volume  scattering  coetliciem 
alone 


I.  INTRODUCTION 

Since  visual  estimates  of  daytime  visibility  are  made 
and  recorded  routinely  by  professionals  and  amateurs 
for  a  variety  of  purposes,  it  ts  a  term  one  might  ex¬ 
pect  to  be  fully  and  clearly  understood  by  all 
Unfortunately,  this  ts  not  always  the  case,  even  among 
experienced  personnel  working  in  the  held  of  atmos¬ 
pheric  sciences.  Thus,  even  though  the  modern  litera¬ 
ture  on  Visibility  is  rich,  extending  from  the  179  entry 
annotated  bibliography  of  Grimes  (1969)  through  an 
additional  700  computerized  references  listed  at  the 
Visibility  Laboratory,  it  was  decided  (hat  yet  anothei 
note  might  be  helpful  if  it  could  provide  a  reasonably 
concise  juxtaposition  of  the  conceptual  variables  defin¬ 
ing  the  psychophysical  aspects  implied  in  the  term 
‘visibility’',  with  the  physical  characteristics  of  some 
equivalent  instrumental  measurements. 

This  note  then  will  first  present  a  brief  review  of  the 
cogent  concept u; '  variables  related  to  visual  de¬ 
terminations  and  their  relationship  to  instrumental 
determinations  of  visibility  It  will  then  illustrate  a 
series  of  nephelometer  measurements  which  are  com 
monly  associated  with  instrumental  determinations  of 
visibility,  but  which  in  this  case  can  be  extended  to 
address  the  more  general  subject  of  determining  slant 
path  directional  contrast  transmittance,  and  the 
modelling  of  the  atmospheric  properties  most  in¬ 
fluencing  these  determinations. 

•  Paper  presented  at  the  Symposium  on  Plumes  and 
Visibility:  Measurements  and  Model  C  omponents  Grand 
Canyon,  Arizona,  U  S  A  It)  14  November  1980 


2.  DAYTIME  VISIBILITY  FROM  V  ISl  At.  OBSFRV  AT  IONS 

The  determination  of  a  human  observ  er's  capability 
lo  discern  barely  perceptible  changes  in  ihe  apparent 
luminance  difference  between  an  object  and  its  sur¬ 
rounding  background  lie  the  luminance  contrast 
threshold  I.  has  been  a  fruitful  held  for  vision  research 
for  many  years.  An  extensive  experimental  effort  was 
initiated  at  the  Tiffany  foundation  in  New  York 
during  the  early  years  of  World  War  II.  the  data  lor 
which  were  tabulated  by  Duntlev  (1946)  and  the 
experimental  results  for  which  were  discussed  by 
Blackwell  11946)  An  important  point  well  illustrated 
by  these  data  is  that  the  contrast  threshold  for  the 
daylight  adapted  human  eye  changes  very  little  as  a 
function  of  background  luminance,  but  is  a  slrong 
function  of  the  angular  subtense  of  the  target  and  of 
the  time  spent  looking  at  the  target. 

To  illustrate  this  rather  fundamental  psychophysical 
relationship,  we  have  generated  the  composite  plots 
shown  in  Fig.  1.  In  this  figure,  the  contrast  values 
shown  along  the  horizontal  axis  are  defined  by  Equa¬ 
tion  (I),  and  in  each  case  have  been  adjusted  from 
the  experimental  liminal  values  (i.e..  50  prob  of 
detection)  to  values  appropriate  for  confident  sighting 
(99  confidence).  The  correction  factors  are  from 
Taylor  (1964).  Equation  (I)  thus  represents  the  most 
fundamental  relationship  involved  in  any  v  isual  search 
scenario 
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where  (  ,  represents  the  apparent  contrast,  at  distance 
r.  of  the  visual  target  t  against  its  background  b.  tH, 
represents  the  apparent  luminance  of  the  usual  target 
as  seen  from  distance  r:  and  hB,  represents  the  apparent 
luminance  of  the  background  against  which  the  target 
is  viewed  as  seen  from  distance  r. 

In  fig.  I.  the  curves  to  the  left,  labeled  Ai  =  x  arc 
from  the  I  iflany  data  lor  long  duration  viewing,  where 
the  viewing  times  At  were  24  s  and  longer.  The  curve  to 
the  right,  labeled  At  =  I  3  s  is  from  Taylor  1 19641.  As 
one  can  note  Irom  these  graphical  displays,  the  Tilfany 
long  duration  data  indicate  via  the  slowly  spreading 
shaded  area,  that  large  changes  in  the  luminance 
of  the  daylight  hori/on  background  (i.e  34(X)  to 
34<>luU  ‘m  : l  make  little  dillerence  to  the  discernible 
threshold  contrast  (  ,.  However,  as  illustrated  by  the 
horizontal  offset  between  the  A(  -•  x  and  the  At 
=  I  3s  plots,  the  time  used  in  viewing  the  object  does 
have  an  appreciable  influence  on  the  contrast 
threshold,  particularly  tor  targets  at  the  larger  angular 
sizes. 

These  data  also  clearly  illustrate  the  underlying 
rationale  for  the  WMO  (I97||  recommendation:  In 
order  that  daytime  visibility  measurements  should  be 
representative,  they  should  be  made  using  objects 
subtending  an  angle  of  not  less  than  0.5  degrees  at  the 
observers  eye  "  The  0.5  deg  angular  size  (eg  equiv¬ 
alent  to  the  solar  or  lunar  disci  and  the  commonlv 
recommended  threshold  contrast  of  0.025  are  clearly 
related  through  the  Taylor  At  13s  data 

3.1  Hi'turtnil  i  isutil  ttnirunu-ntiH  t  nmptiri\itn\ 

Armed  with  the  <j  priori  knowledge  that  target  size 
dcscernible  threshold  contrast,  and  observer  glimpse 


limes  are  all  interrelated,  and  olt  times  less  than  clearly 
detined.  one  might  fairly  ask:  "What  threshold  contrast 
is  actually  limiting  the  determination  of  the  visibility 
when  an  observer  just  sights  a  particular  target?" 

fortunately,  there  have  been  a  good  number  of  well 
done  experiments  wherein  the  visibility  range  was 
determined  visually  at  the  same  time  as  the  de¬ 
gradation  of  the  atmosphere  was  measured  by  physical 
instruments.  If  one  makes  a  few  assumptions  about  the 
performance  of  the  man  and  the  machine  in  these 
experiments,  some  encouragingly  consistent  results 
emerge.  Most  caveats  regarding  these  assumptions  are 
discussed  at  length  by  others  as  recently  as  Malm 
(1979)  and  need  not  be  pursued  further  here  Thus,  il 
one  assumes  that  the  visual  determinations  were  made 
using  a  truly  black  target  having  an  inherent  contrast 
of  -  I  against  a  clear,  cloud-free  horizon:  and  if  one 
also  assumes  that  the  machine  measures  an  attenuation 
coefficient  x  which  is  representative  ol  the  average 
coefficient  along  the  path  of  sight,  then  the  apparent 
contrast  is  given  by 

3  .  1 ,  e  ”,  (21 

where  /'.  is  the  transmittance  ot  the  path  of  sight  r 
along  which  ihe  determination  ol  t  ,  is  accomplished 
In  liquation  < 2 1.  one  notes  that  the  visibility  range  r 
is  defined  once  a  threshold  contrast.  (  ,  is  specified  It  is 
through  this  expression  that  the  perlormance  link 
between  man  and  machine  can  be  evaluated,  and  in  fact 
answer  the  initially  posed  question  regarding  the 
actual  limiting  threshold  contrasts 

Ihe  experimental  wotks  ol  Douglas  am  3oung 
1 1 945 1,  Punt  ley  1 1 94K i.  Middleton  i  I9s2i.  Hor  >ath  and 
Noll  iI9is9i  and  Mermg  i :  ,/l  i|9~||  all  relate  to 
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determining  the  efficacy  of  the  relationships  specified 
in  Equation  (2).  The  data  in  general  are  all  similar  in 
nature  to  those  of  Heringet  at.  (19711  illustrated  in  Fig- 
2,  which  compare  extinction  coefficients  determined 
from  transmissometer  measurements  with  visibilities 
determined  by  human  observers. 

If  there  is  a  reasonably  constant  threshold  contrast 
|C'r|.  then  data  of  this  sort  can  be  expected  to  fall  on  a 
straight  line  when  plotted  in  logarithmic  formats 
which  reflect  Equation  |2|  or  us  simplified  equivalent 

-  In  =  xr,  (?) 

which  is  merely  another  way  of  saying  that  if  the 
threshold  contrast  is  constant,  then  the  product  xr. 
the  threshold  sighting  range  times  the  attenuation 
coefficient,  is  also  constant 

In  Fig.  2,  the  straight  lines  drawn  through  the  data 
for  comparative  purposes  represent  the  conditions 
which  would  be  expected  to  exist  for  constant  thres¬ 
hold  contrasts  of  0  02  and  0.05 

There  is  a  satisfying  thread  of  constancy  running 
through  these  experimental  works  which  illustrates  the 
reasonableness  of  expecting  a  typical  threshold  con¬ 
trast  for  human  daylight  observations  to  exist  and 
that  it  be  analytically  relatable  to  instrumental 
measurements. 

Early  workers  in  the  field  of  visual  search  recom¬ 
mended  the  use  of  visibility  targets  at  least  0.5  deg  in 
size.  and  the  corresponding  use  of  a  threshold  contrast 
of  0025  which  would  thus  yield  the  attenuation 
coefficient-sighting  range  product  xr  of  3.7.  However, 
il  one  chooses  to  adopt  a  threshold  contrast  of  0.05. 
then  the  corresponding  xr  product  becomes  3.0  and 
additional  experimental  consistencies  emerge.  Accept- 
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mg  a  variability  of  ±20",,  in  visually  determined 
sighting  ranges  results  in 

xr  =  3.0  +  0.6.  (4| 

with  a  corresponding  variability  in  threshold  contrast 
of  0.027  to  0.091.  This  threshold  relationship  is 
illustrated  on  Fig  I  by  the  broken  lines  connecting  the 
horizontal  and  vertical  axes  via  the  Al  =  1  3  s  curve 
From  this  plot  it  can  be  seen  that  threshold  angular 
sizes  might  be  expected  to  vary  between  3.6  and 
350  min  of  arc 

An  examination  of  the  experimental  data  referenced 
in  the  preceeding  paragraphs  with  the  visibility  specifi¬ 
cation  implied  in  Equation  1 4 )  in  mind,  indicates  that 
by  far  the  bulk  of  these  experimental  data  fall  within 
the  bounds  ol  this  specification  In  a  way.  this  is  a 
remarkable  consistency  considering  the  variability  in 
the  mechanics  of  the  various  determinations,  and  the 
conditions  under  which  they  were  conducted 

it  is  thus  a  well  supported  contention  from  both 
psychophysical  and  physical  measurement  con¬ 
siderations.  that  usibiliiv  can  reasonably  be  de¬ 
fined  as  the  distance  that  a  black  target  of  suitable 
dimensions  can  be  seen  against  the  cloudless  horizon 
sky.  and  it  is  related  to  the  attenuation  coefficient  with 
an  accuracy  of  *  20' ,  by  the  relationship  indicated  in 
Equation  (4).  There  are.  however,  specific  penalties 
involved  if  one  chooses  to  ignore  the  caveats  associated 
with  visual  determinations  of  horizontal  visibility. 

2.2.  Siin-viutidurd  e/fi-ers 

2  2  1  Son-black  hiri/ft>  It 's  ol  interest  to  note  the 
influence  on  visibility  estimates  induced  through  the 
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use  of  either  non-black  targets,  or  backgrounds  other 
than  the  horizon  sky.  Derivations  to  conveniently 
evaluate  these  influences  are  given  in  Gordon  (1979) 
from  which  several  summarizing  illustrations  are  re¬ 
produced  herein.  As  a  convenience,  we  have  chosen  to 
use  the  reflectance  format  for  most  variables,  hence 
inherent  contrast  is  re-defined  as  a  function  of  the 
target  to  background  reflectance  ratio, 


Co 


-  —  1, 


(5) 


w here  C0  represents  the  preceived  contrast  at  distance 
zero  from  the  visual  target:  and  ,RU  and  bR0  represent 
the  respective  target  and  background  reflectances.  A 
similar  definition  is  of  course  appropriate  for  the 
apparent  contrast, 

The  effect  on  visual  range  determinations  of  using 
non-black  targets  is  illustrated  in  Fig.  3.  In  this  plot,  the 
horizontal  axis  is  conveniently  marked  in  both  in¬ 
herent  reflectance  ratio  and  the  equivalent  inherent 
contrast.  The  vertical  axis  is  marked  as  Relative  Visual 
Range,  which  is  defined  as  the  ratio  of  the  visual  range 
one  would  determine  observing  a  non-black  target,  to 
that  one  would  determine  if  observing  an  ideal  black 
target  having  the  optimum  inherent  contrast  of  minus 
one. 

It  is  apparent  that  for  the  relative  visual  range 
determination  to  be  within  5  of  the  ideal,  negative 
contrast  can  range  between  -  1  and  -0.86  However, 
for  varying  positive  contrasts,  the  relative  visual  range 
changes  quite  rapidly,  and  the  5  variance  from  ideal 
is  contained  in  a  very  narrow  band  between  +  1.9  and 
+  2.2. 


2.2.2.  Mon-horizon  sky.  For  visual  range  determi¬ 
nations  using  non-horizon  sky  backgrounds.  Mid¬ 
dleton  11952)  presented  a  series  of  expressions  for 
comparing  the  relative  visual  range  of  a  black  object  as 
seen  against  a  black  background  with  that  of  the  black 


object  as  seen  against  the  horizon  sky.  His  analysis 
indicated  that  for  a  threshold  contrast  of  0.02,  an 
acceptably  accurate  determination  of  visual  range 
would  result  if  the  black  background  was  at  least  1.5 
times  as  distant  from  the  observer  as  was  the  object. 
Since  this  distance  ratio  is  broadly  recognized  and  is 
referred  to  as  a  general  rule  of  thumb  by  the  WMO 
( 1 97 1 ),  it  is  appropriate  at  this  point  to  enlarge  slightly 
upon  the  previous  references. 

if  one  imagines  the  black  object  and  its  black 
background  initially  at  the  same  distance  from  the 
observer,  then  it  is  intuitively  clear  and  analytically 
correct,  that  in  this  orientation  there  is  zero  contrast 
between  them  However,  if  one  now  begins  to  move  the 
background  progressively  farther  and  farther  away, 
the  additional  path  radiance  generated  in  the  lengthen¬ 
ing  observer  to  background  path  will  begin  to  develop 
a  difference  in  apparent  luminance  between  the  two 
black  surfaces.  As  the  background  distance  lengthens, 
the  apparent  contrast  between  the  object  and  back¬ 
ground  will  increase  until  at  some  p>oint  it  reaches 
the  observers  threshold  contrast.  When  the  back¬ 
ground  has  reached  this  point,  the  determination  of 
the  black  objects  visual  range  will  be  the  same  as  if  it 
had  been  determined  against  the  horizon  sky. 

In  her  derivations  for  non-standard  effects.  Gordon 
(1979)  has  extended  this  general  class  of  analysis 
beyond  the  example  of  Middleton  to  include  a  selec¬ 
tion  of  backgrounds  other  than  the  black,  and  to 
establish  the  relative  visual  ranges  appropriate  for  a 
threshold  contrast  of  0.05.  A  key  parameter  in  her 
evaluation  is  the  ratio  between  the  background  re¬ 
flectance  and  the  equilibrium  reflectance  bR0/Rq.  The 
equilibrium  reflectance  R,  is  related  to  the  equilibrium 
radiance  Nr  i.e..  the  radiance  of  the  clear  horizon  sky, 
Duntley  (1948),  as  expressed  in  Equation  (6). 


R,(0.90  .</»  = 


\,(0,90  .<p\n 
H  (0  90  .0+180  ) 


(6) 
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fig  3  Relative  visual  range  as  a  function  of  inherent  contrast  and 
inherent  target  lo  background  reflectance  ratio  for  a  contrast  threshold 
of  005 
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where  N f  (0, 90 <j>)  represents  the  clear  horizon  sky 
radiance  as  viewed  at  zero  altitude  in  the  direction  <p 
from  the  sun;  H(0, 90c,  <#>  +  180“)  represents  the  2ir 
irradiance  falling  upon  a  vertically  oriented  flat  plate 
collector  whose  normal  is  directed  in  azimuth  toward 
<p  + 180°;  and  where  in  all  cases,  the  parenthetical 
notations  denote  the  altitude,  zenith  angle  and  azimuth 
angles  of  the  measured  path  of  sight. 

As  can  be  seen  in  Fig.  4  which  has  been  derived  from 
the  data  of  Gordon  (1979),  the  use  of  non-black 
backgrounds  which  have  reflectances  less  than  equilib¬ 
rium  (6R0/R,  <  1)  at  the  same  relative  distance  beyond 
the  target  yield  relative  visual  ranges  even  closer  to  one 
than  does  the  use  of  black  background.  For  back¬ 
grounds  brighten  than  equilibrium  (fcR„  R9  >  1). 
such  as  a  bright  cloud,  it  is  a  different  matter.  With  the 
use  of  these  bright  backgrounds  it  is  necessary  to  have 
the  background  approximately  twice  as  distant  as  the 
target  in  order  for  the  relative  visual  range  to  become 
reasonably  accurate. 

2.2.3.  Non-horizontal  paths  oj  siyfu.  As  one  departs 
from  the  special  case  of  determining  contrast  transmit- 
tances  along  horizontal  paths  through  homogeneous 
atmospheres  to  the  more  general  case  of  slant  path 
transmittances.  additional  complexities  emerge.  One 
must  now  consider  the  variations  in  at'enuation 
coefficient  a  as  a  function  of  altitude  along  the  inclined 
path,  as  well  as  the  directional  scattering  propei  ties  of 
the  aerosol  through  which  the  determination  is  to  be 
made.  It  is  at  this  point  that  instrumental  determi¬ 
nations  become  the  preferred  vehicle  since  they  can 
provide  this  additional  characterization  of  the  atmos¬ 
pheric  condition  in  a  manner  superior  to  attempts  by 
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visual  observation.  Experimental  measurements  of  the 
sort  required  for  these  slant  path  determinations  are 
presented  in  the  following  section. 

3.  PHYSICAL  MEASUREMENTS 

The  direct  physical  measurement  of  air  clarity  in  a 
manner  appropriate  for  realistic  comparison  with 
determinations  made  by  human  observers  has  long 
been,  and  still  is  an  endeavor  spiced  with  gentle  but 
lively  controversy.  Of  the  many  techniques,  the  most 
common  are  those  related  to  the  measurement  of 
lumina  ce  transmittance  over  some  specified  path 
length,  the  measurement  of  the  point  function  at¬ 
tenuation  or  scattering  coefficient,  or  the  measurement 
of  the  apparent  contrast  of  distant  targets.  In  each  case, 
the  particular  procedure  will  possess  inherent  short¬ 
comings  as  discussed  in  Malm  (1979),  and  may  or  may 
not  relate  well  to  the  total  human  perception  of  visual 
stimuli  as  discussed  in  Malm  ei  al.  (1980). 

Although  a  measured  value  of  visible  spectrum 
atmospheric  scattering  coefficient  is  often  used  with 
Equation  (4)  for  ihe  determination  of  daytime 
visibility,  if  one  is  addressing  the  more  general  case  of 
atmospheric  slant  path  contrast  transmittance,  more 
information  is  required.  One  must  obtain  a  specifi¬ 
cation  of  both  the  directional  scattering  properties  of 
the  atmospheric  aerosols,  and  their  vertical  and  horizon¬ 
tal  distnbuiion.  A  multi-channel  nephelometer  has 
been  used  by  the  Visibility  Laboratory  to  document 
these  specific  atmospheric  properties. 

3.1  Airborne  nephelometer  design 

There  are  a  variety  ofclectro-optical  devices  designed 
to  measure  the  energy  scattered  from  an  illuminated 


HAT  If)  or  BACKiiKOUNO  TO  IAROET  OISTANCE 

fig  4  Relative  visual  ranges  for  black  targets  viewed  against  ba¬ 
ckgrounds  olher  than  clear  horizon  sky  (threshold  conlrasl  =  0.05) 
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volume  of  aerosol  that  are  calibrated  in  terms  of 
risibility.  The  EG&G  Forward  Scatter  Visibility 
Meter  described  in  Hering  et  al.  1 1971 1  is  one  example. 
Integrating  nephelometers  similar  to  those  discussed 
in  C'harlson  el  al.  ( 19671  are  additional  examples  of 
devices  based  upon  one  or  the  other  of  the  two 
techniques  described  by  Beuttell  and  Brewer  (19491.  In 
all  cases,  these  scattering  meters  utilize  the  relationship 

s|c|  =  I  fllc./fldO.  (7) 

where  <r(c.  ji\  is  the  volume  scattering  function  at 
scattering  angle  /(  and  r  is  an  altitude  parameter. 

Both  of  the  procedures  described  by  Beuttell  and 
Brewer  (1949)  perform  the  integration  of  Equation  (7) 
optically  over  various  fractions  of  the  4 n  solid  angle  so 
as  to  include  measurements  o\er  as  broad  a  range  of 
scattering  angles  as  possible.  F  undamentally  the  two 
instrumental  procedures  are  optical  inverses  of  one 
another.  The  rtrst  uses  a  diffusely  emitting  light  source 
to  illuminate  a  sample  volume  of  aerosol  w  hich  is  then 
observed  via  a  small  held  telescope  assembly  as  in  the 
work  of  C  rosby  and  Koerber  11963).  C'harlson  ei  ul. 
I1967|  and  Horvath  and  Noll  (1969)  The  second 
reverses  the  optical  collector  emitter  relationship  by 
using  a  well  collimated  projector  as  the  illuminating 
source  and  a  cosine  corrected  irradiomeier  to  observe 
the  scattered  flux.  This  second  design  has  been  used  by 
Duntlev  et  ill.  (1970.  I97X) 

The  Visibility  Laboratory  nephelometer  func¬ 
tionally  illustrated  in  Fig.  5.  is  unique  in  several  design 
aspects  which  make  it  particularly  attractive  experi¬ 


mental  tool.  First,  the  optical  system  used  a  cylindncally 
limited  projector  beam,  stopped  to  provide  a  rec¬ 
tangular  beam  cross  section  which  provides  good 
geometrical  definition  of  the  illuminated  volume.  The 
observed  scattering  volume  is  approximately  5  cm 
x  9  cm  in  cross  section,  and  107  cm  in  length  for  a  total 
sample  volume  of  48 1 5  cm  Second,  the  detector 
assembly's  primary  optica)  channel  uses  a  cosine 
corrected  and  mechanically  shielded  irradiometer  head 
w  hich  performs  the  integration  of  Equation  (7)  over  all 
scattering  angles  between  5  and  172  minimizing 
truncation  losses  of  the  sort  discussed  in  Fitzgerald 
(1977).  Rabinoff  and  Herman  (1973)  and  Ensor  and 
Waggoner  (1970).  Third,  the  detector  assembly's  second¬ 
ary  optical  channels  measure  relatively  narrow  angle, 
i.e.,  2  held  of  view,  directional  scattering  at  both  30 
and  150  which  allows  the  determination  of  a  forward 
to  backward  scattering  ratio  that  can  be  used  to  more 
thoroughly  characterize  the  sample  aerosol.  Fourth, 
the  device  is  multi-spectral  w  ithin  the  visible  spectrum 
Thus  visible  spectrum  blue-red  scattering  ratios  in 
addition  to  photopically  corrected  measurements  can 
be  accomplished.  And  filth,  the  entire  device  is  de¬ 
signed  to  operate  as  a  ratio  measuring  system  which 
tends  to  minimize  calibration  errors  and  stability 
problems 

The  calibration  constants  for  the  device  relate  the 
ratio  of  the  flux  scattered  from  the  sample  aerosol  to 
the  flux  reflected  from  a  calibrated  reflective  plaque 
with  the  physical  properties  of  the  overall  optical 
system,  and  the  most  probable  shape  of  the  scattering 
function  associated  w  ith  the  sample  aerosol  Constants 
have  been  derived  for  scattering  functions  ranging 
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from  the  Rayleigh  molecular  through  ten  increasingly 
steep  functions  representing  a  broad  range  of  real 
atmospheric  aerosols  as  determined  from  Barteneva 
(I960).  Compensation  for  the  approx.  7"„  variation  in 
calibration  due  to  aerosol  variations  is  achieved  by 
selecting  the  most  appropriate  constant  on  the  basis  of 
the  nephelometers  directional  scattering  measure¬ 
ments. 

Inasmuch  as  the  size  and  configuration  of  the 
Visibility  Laboratory  nephelometer  preclude  regular 
calibration  by  purging  the  sample  volume  with  known 
gases  as  in  Charlson  et  ul.  (1967)  and  Bodhaine  (1979). 
secondary  checks  to  validate  the  relative  integrity  of 
the  measurements  have  been  an  important  and  con¬ 
tinuous  facet  of  our  experimental  procedure. 

Coincident  ground  based  measurements  made  by  an 
integrating  nephelometer  and  a  long  path  transmisso- 
meter  offei  a  convenient  cross  check  under  conditions 
of  no  atmospheric  absorption  and  reasonably  uniform 
horizontal  stratification.  The  data  in  Fig.  6.  illustrate  a 
set  of  comparison  measurements  made  between  the 
Visibility  Laboratory  integrating  nephelometer  and  an 
Eltro  transmtssometer  similar  to  that  illustrated  in 
Douglas  et  ul.  (1977)  during  a  joint  operation  tn 
northern  Germany.  Duntley  el  at.  (1976).  These  data 
represent  two  pairs  of  measurements  made  on  each  of 
two  separate  days,  and  show  good  comparison  be¬ 
tween  the  two  system  performances. 

3.2.  Experimental  measurements 

3.2.1.  Comparisons  with  visual  observations.  One 
must  keep  in  mind  that  the  use  of  stable  and  reliable 


instrumentation  does  not  imply  that  the  measurements 
one  gathers  will  describe  a  stable  optical  environment. 
Particularly  in  the  lower  troposphere,  there  are  sig¬ 
nificant  spatial  and  temporal  variations  in  the  optical 
state  of  the  atmosphere  that  make  machine  vs.  machine, 
or  man  vs.  machine  comparisons  difficult.  A  simple 
illustration  of  this  variability  is  shown  in  Fig.  7.  In  this 
illustration  each  rectangle  represents  the  maximum 
variation  in  derived  and  reported  visibilities  that 
occurred  during  seventeen  different  data  missions.  The 
rectangular  elements  are  typified  by  time  intervals  of 
approx.  2h,  and  spatial  intervals  of  approx.  60  km. 
This  time-space  interval  is  typical  of  that  used  during 
our  experimental  flight  program  for  the  accumulation 
of  one  set  of  ground  level  to  6  km  measurements.  Each 
rectangle  represents  one  of  17  separate  flights.  Thus  it 
is  clear  that  even  though  the  machine  and  the  observer 
may  both  be  doing  their  jobs  properly  as  the  trend  of 
the  data  in  Fig.  7  might  indicate,  it  is  essential  that  they 
be  exactly  coincident  in  the  time  and  site  of  their 
performance  if  one  is  to  expect  precise  correlation. 

3.2.2.  Vertical  variation  oj  atmospheric  aerosols.  If 
one  is  to  address  the  generalized  problem  of  slant  path 
contrast  transmittance  through  the  atmosphere,  then 
in  addition  to  determining  a  single  value  of  scattering 
coefficient  which  is  suitable  for  approximating  hori¬ 
zontal  visibility  via  Equation  (4).  the  spatial  variation  in 
scattering  coefficient  with  changing  altitude  is  also 
essential.  These  variations  have  been  routinely  docum¬ 
ented  as  reported  in  Duntley  et  at.  (1978)  over  large 
portions  of  the  U  S  A.  and  Europe. 

Two  examples  illustrating  the  broad  variety  in 
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Fig  6  Comparison  of  nephelometer  n-d  Itansmissometet  data.  Meppen. 
Germany  (Duntley  it  ul.  19721 
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VISIBILITY  (km)  REPORTED  BV  METEOROLOGICAL  STATIONS 
WITHIN  57  km  OF  TRACK  OR  VIS  LAB  GROUND  STATION 


Fig.  7  Comparison  of  visibilities  derived  from  nephelometer  data  (from 
Equation  |4|  where  visibility''  r  =  J,j|  with  those  reported  by  nearby 
meteorological  stations 


vertical  structure  that  one  might  anticipate  are  shown 
in  Figs  8a  and  8b.  In  these  multi-spectral  displays  the 
filter  identification  code  is  as  follows:  Filter  2,  mean 
wavelength  at  447  nm;  Filter  3,  mean  wavelength  at 
663  nm:  Filter  5,  mean  wavelength  at  750  nm,  and  Filter 
4,  mean  wavelength  at  550  nm,  representing  the  pseudo 
photopic  response  which  relates  directly  to  the  sensi¬ 
tivity  of  the  human  eye. 

Flight  C-372  illustrates  the  commonly  found  con¬ 
dition  of  a  sharp  naze  top  at  about  2  km  above  ground 
level  separating  the  relatively  heavy  low  altitude  haze 
from  the  clear  air  above.  It  also  illustrates  an  additional 
feature  not  commonly  shown  as  clearly  as  in  this 
example.  Whereas,  one  normally  expects  the  scattering 
coefficient  to  decrease  as  altitude  increases,  it  is  not 
uncommon  to  find  instances  of  increasing  haziness  and 
thus  increasing  coefficients  as  one  climbs  from  the 
surface  to  near  the  haze  top.  This  characteristic  is 
clearly  illustrated  in  the  data  for  Filter  3  and  is 
normally  one  found  in  profiles  exhibiting  relatively 
strong  haze  top  demarcations. 

Flight  C-155  illustrates  the  extreme  opposite  con¬ 
dition  which  is  far  more  appropriate  to  those  con¬ 
cerned  with  monitoring  and  preserving  our  western 
desert  environment.  These  data  were  obtained  south  of 
Albuquerque,  New  Mexico  during  a  morning  flight 
under  clear  conditions,  with  the  flight  crew  reporting 
light  haze  on  the  horizons.  These  are  probably  as 
idealized  a  set  of  profiles  as  one  might  ever  find.  There 
is  little  or  no  structure  with  altitude,  and  as  indicated 
by  the  high  degree  of  parallelism  between  the  data  in 
each  filter,  there  was  exceptional  temporal  stability 
within  the  measurement  interval. 


Examination  of  profiles  of  the  sort  illustrated  in  Fig. 
8(a)  indicates  that  it  is  not  uncommon  for  low  altitude 
changes  in  scattering  coefficient  to  exceed  a  factor  of 
ten  within  the  last  IGOOm  above  the  terrain.  The 
magnitude  and  altitude  of  this  often  abrupt  haze 
boundary  are  obviously  major  influences  upon  the 
determination  of  slant  path  transmittances  through 
this  lower  altitude  regime.  It  is  equally  apparent  that 
few  atmospheres  are  as  uniform  and  well  behaved  as 
that  represented  by  flight  C-155,  and  that  to  assume  so, 
as  is  often  done  in  simple  modelling  exercises,  may  lead 
to  significant  error. 

3.2.3.  Directional  scattering  of  atmospheric 
aerosols.  The  measurement  of  two  discrete  values  of 
atmospheric  directional  scattering  function  coin- 
cidently  with  the  measurement  of  total  volume  scatter¬ 
ing  coefficient  is  a  powerful  capability  that  is  built  into 
the  Visibility  Laboratory  integrating  nephelometer.  It 
is  essential,  if  one  is  to  describe  the  optical  properties  of 
an  atmosphere,  either  for  a  specific  contrast  transmit¬ 
tance  determination,  or  for  a  general  purpose  model¬ 
ling  application,  that  a  reliable  entry  into  the  directional 
characteristics  of  those  properties  be  available.  The 
work  of  Barteneva  (I960)  has  been  heavily  utilized  by 
the  Visibility  Laboratory  as  this  entry,  but  there  have 
been  continuous  queries  as  to  the  complete  generality 
one  could  assign  to  the  Barteneva  data  (Johnson  et  al., 
1979).  Measurements  made  with  our  instrument  have, 
I  believe,  answered  that  question  beyond  reasonable 
doubt. 

Barteneva  (1960)  provides  a  catalog  of  proportional 
directional  scattering  functions  <r(/))/s  for  the  photopic 
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Fig  8  Airborne  measurements  of  volume  scattering  coefficient 
profiles  under  typically  ha/>  conditions  and  under  clear,  cloud 
tree  conditions. 


sensor  based  upon  624  measurements  made  during 
1 955  1958  at  various  locations  in  the  U.S.S.R.  and  at 
sea.  For  each  volume  scattering  function  she  gives  a 
range  of  total  volume  scattering  coefficients  appro¬ 
priate  to  that  function.  Median  values  of  scattering 
coefficient  were  derived  for  each  of  the  gradual  scatter¬ 
ing  function  classes,  Nos  110.  This  resultant  catalog 
of  gradual  functions  and  their  associated  median 
scattering  coefficients  is  appropriate  for  describing 
atmospheres  characterized  as  very  clear  to  heavy  haze, 
Nos  I  to  10  respectively. 

It  became  apparent  during  the  analysis  of  the 
measurements  from  the  directional  channels  of  Visi¬ 
bility  Laboratory  nephetometer  that  its  data  compared 
well  to  the  Barteneva  data.  When  its  median  scattering 


coefficient  values  and  its  gradual  class  ir(/f)/s  curves 
were  put  into  the  form  of  their  ratio  to  the  Rayleigh 
ti  e.,  molecular  scattering,  as  was  the  nephelometer 
data),  then  there  was  exoellent  correlation.  Two  sets  of 
these  comparative  data  are  illustrated  in  Figs  9  and  10. 

It  is  of  significance  to  note  that  whereas  the  original 
Barteneva  data  represent  photopic  measurements  at 
ground  level,  the  Visibility  Laboratory  measurements 
include  data  from  several  separate  bands  within  the 
visible  spectrum  and  from  within  a  broad  range  of 
attitudes. 

In  Figs  9  and  10.  the  horizontal  axis  is  Optical 
Scattering  Ratio  Q(r)  which  is  defined  as  the  ratio  of 
the  total  volume  scattering  coefficient  s(z)  to  the 
Rayleigh  scattering  coefficient.  „s(r)  where  r  is  the 
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Fig.  9.  Comparison  of  multi-spectral  directional  seal* 
tenng  ratios  measured  b\  Visibility  Laboratory  air¬ 
borne  nephclometcr  with  the  photopic  ratios  I'rom 
Barteneva  (19601 

OPAQUE  V  GROUND  DATA 


I  ig  Id  tompariNon  ol  multi-spcctral  directional 
scattering  ratios  measured  by  Visibility  Laboratory 
ground- based  nephclometcr  with  the  photopic  ratios 
from  Barteneva  (I960! 


altitude  designator: 
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The  vertical  axis  is  Volume  Scattering  Function  Ratio, 
q(z)  which  is  the  ratio  of  the  total  directional  scattering 
function  at  angle  beta,  a{:,  fi )  to  the  Rayleigh  scattering 
function  at  the  same  angle  and  altitude.  R<x(r, 


aU.fi) 

aU.fi) 


Figure  9  represents  multi-spectral  directional  scat¬ 
tering  function  data  measured  during  fourteen  flights 
made  in  the  vicinity  of  St.  Louis.  Ml.,  during  the 
summer  of  1971.  In  each  case  the  measurements  made 
at  scattering  angles  of  30  and  150  are  plotted  as  i/(rl 
against  the  coincidently  measured  values  of  total 
scattering  coefficient  converted  to  values  of  £>(;).  The 
solid  lines  are  the  equivalent  data  extracted  from  the 
Barteneva  catalog.  The  correlation  is  good  over  the 
entire  range  of  G(-) 

Figure  10  represents  additional  multi-spectral  data 
of  the  same  sort  as  in  Fig.  9.  that  was  measured  at 
ground  level  in  conjunction  with  19  flights  made  in 
western  Europe  during  the  summer  of  1978.  Again,  the 
correlation  is  strong  over  a  broad  range  of  Gl-I 

The  excellent  comparisons  of  both  the  in-flight  and 
surface  measurements  with  the  Barteneva  directional 
data  as  illustrated  in  Figs  9  and  10  argue  powerfully  for 
accepting  the  general  applicability  of  the  Barteneva 
scattering  function  catalog.  It  is  this  experimental 
verification  of  genera)  applicability  that  establishes  the 
practicality  of  developing  simple  and  reliable  models 
of  atmospheric  directional  scattering  properties. 


4.  C  OM  l.l  SKINS 

Although  many  early  experimentalists  established 
the  viability  of  visual  determinations  of  daytime 
contrast  transmittance  along  horizontal  paths  of  sight 
and  related  this  determination  to  the  concept  of 
"visibility",  there  are  equally  well  documented  and 
sometimes  not  as  fully  understood  associated  caveats 
However,  within  the  boundaries  established  by  these 
caveats,  horizontal  visibility  determinations  derived 
from  instrumental  measurements  of  atmospheric  ex¬ 
tinction  coelficienl  are  adequately  relatable  to  visuai 
observations. 

The  extension  of  these  concepts  to  the  more  general 
case  of  contrast  transmittance  determinations  along 
vertically  inclined  paths  of  sight  implies  knowledga- 
bility  of  the  vertical  variability  in  atmospheric  extinc¬ 
tion  coelficienl  and  the  directional  scattering  proper¬ 
ties  along  the  path  These  data  have  been  obtained 
through  the  use  of  a  multi-channel  nephelometer  and 
indicate  a  strong  correlation  between  the  measured 
scattering  coefficient  and  the  aerosol  directional  scat¬ 
tering  function. 

The  modelling  implications  of  these  data  are  indeed 
enticing.  The  evidence  to  date  strongly  supports  the 
contention  that  with  an  adequate  specification  of  total 
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volume  scattering  coefficient  s(r)  and  its  Rayleigh 
component,  one  can  enter  the  Barteneva  catalog  or  its 
parameterized  equivalent,  and  obtain  a  truly  repre¬ 
sentative  expression  of  the  directional  scattering  fun¬ 
ction  o{z,  jj)  for  any  spectral  band  within  the  visible 
spectrum.  Thus,  a  reliable  specification  of  the  direc¬ 
tional  scattering  properties  of  the  atmospheric  aerosol 
can  be  reasonably  deduced  from  a  measurement  of  the 
total  volume  scattering  coefficient  alone. 
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ABSTRACT 

The  radiation  emerging  from  the  top  of  the  earth's  atmosphere  is  affected  by  the  reflection  characteristics 
of  the  underlying  surface.  Laboratory-gathered  bidirectional  reflectance  data  were  used  to  characterize 
the  reflection  matrix  for  three  different  types  of  natural  surfaces.  The  radiation  emerging  from  the  top 
of  a  pure  molecular  and  an  aerosol-laden  plane  parallel  atmosphere  was  calculated  for  these  three 
diverse  types  of  surfaces  and  their  Lambert  model  equivalents.  The  calculations  indicate  that  the  radiation 
emerging  from  the  top  of  the  atmospheres  is  significantly  different  for  a  real  surface  compared  to  a 
Lambert  model.  For  the  molecular  atmosphere  the  difference  is  most  pronounced  at  a  wavelength  of 
0.6  pm.  this  is  mainly  due  to  the  light  directly  transmitted  downward  through  the  atmosphere,  reflected, 
and  then  directly  transmitted  outward  through  the  top  of  the  atmosphere.  For  a  turbid  atmosphere  the 
difference  also  is  due  to  light  with  a  history  of  diffuse  transmission  because  the  strong  forward  scattering 
of  the  aerosols  partially  offsets  the  smoothing  effect  of  the  diffuse  multiple  scattering.  The  calculations  indi¬ 
cate  that  the  polarization  of  the  light  emerging  from  the  top  of  the  atmosphere  strongly  mirrors  the 
polarization  reflected  from  the  underlying  surface.  Thus,  satellite-measured  polarization  would  be  an 
excellent  parameter  for  distinguishing  ground  features  with  similar  reflectance  but  different  polariza¬ 
tion  characteristics. 


1.  Introduction 

To  facilitate  the  interpretation  of  data  gathered 
from  satellites  and  aircraft,  many  investigators  have 
made  computations  of  the  radiation  emerging  from 
the  earth's  atmosphere  overlying  various  types  of 
surfaces.  In  most  cases,  the  radiation  reflected  from 
the  surface  was  considered  isotropic  and  unun- 
polarized,  corresponding  to  the  Lambert  approxima¬ 
tion.  The  radiation  emerging  from  the  top  of  an  at¬ 
mosphere  overlying  different  natural  surfaces  has 
been  investigated  by  Coulson  el  ul.  (1966)  using  a 
Rayleigh  atmosphere  and  by  Keopke  and  Kriebel 
( 1978)  using  a  turbid  atmosphere  overlying  vegetated 
surfaces  reflecting  light  of  an  unknown  degree  of 
polarization.  These  investigators  found  significant 
differences  in  this  emergent  radiation  compared  to 
the  radiation  emerging  from  an  atmosphere  over- 
lying  a  Lambert  surface,  even  though  the  albedo  of 
the  Lambert  surface  was  equal  to  the  albedo  of  the 
natural  surface.  Because  surface  reflection  influences 
the  emergent  radiation,  the  characteristics  of  the 
radiation  emerging  from  a  pure  molecular  (Rayleigh) 
and  an  aerosol-laden  atmosphere  overlying  different 
types  of  surface  materials  is  investigated  here. 


‘  Present  affiliation:  Visibility  Laboratory,  Scripps  Institution 
of  Oceanography.  University  of  California,  San  Diego  9209 V 


In  this  paper,  the  radiation  emerging  from  the  top 
of  a  Rayleigh  and  a  turbid  plane  parallel  non¬ 
absorbing  atmospheric  model  overlying  three  types 
of  natural  surfaces  of  known  directional  reflectance 
and  polarization  is  described. 


2.  Atmosphere-surface  model 

The  solar  radiation  emerging  from  the  top  of  an 
atmosphere  illuminated  by  parallel  radiation  7tF  can 
be  expressed  as  an  infinite  series,  with  each  term  in 
the  series  having  its  own  history  of  transmission, 
scattering  and  surface  reflection.  A  few  of  the  terms 
depicting  the  radiation  emerging  from  the  atmosphere 
are  shown  in  Fig.  1 .  By  summing  the  terms  in  the 
series,  the  radiance  I(r;  /x,  <f>)  emerging  from  the  top 
of  a  planetary  atmosphere  of  normal  optical  thick¬ 
ness  t,  can  be  written  as 

I(r  =  0;  n ,  <t>) 

=  SttF  +  e  r|,wERe"T|"1“En-F  +  e  t'"ERT7tF 
first  second  third 

+  T*Re  T'"‘«EjrF  +  T*RTttF  +  <e  ''"“E  +  T*)R 
fourth  fifth 

X  (S*R  4-  S*RS*R  +  •  •  •)(<>  ',U"E  +  T)jtF.  (1) 
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Fig.  I.  Diagrammatical  representation  of  some  of  the  components  of  Ihe  total  emergent  light,  nh  is  the  solar  radiation 
incident  on  the  atmosphere  having  a  total  optical  thickness  t,.  and  S  and  T  are  the  diffuse  reflection  and  transmission 
matrices,  respectively.  The  termr  E  represents  direct  transmission  through  the  atmosphere,  where  E  is  the  unit  matrix. 
The  reflection  characteristics  of  the  surface  are  described  by  R. 


where  0  =  cos "V  is  the  nadir  viewing  angle,  <f>  the 
azimuthal  angle,  I(r;  jz;  <b)  and  F  are  vectors,  and  E 
is  the  3  x  3  unit  matrix.  The  position  of  the  sun,  or 
source,  is  given  by  the  zenith  angle  00  =  cos'Vo 
and  the  azimuthal  angle  <t>0  =  0°.  At  the  top  of  the 
atmosphere  r  =  0. 

For  convenience,  parallel  radiation  of  net  flux 
n-F  per  unit  area  normal  to  the  beam  is  assumed 
incident  on  the  top  of  the  atmosphere.  Using  the 
Stokes  parameter  representation  ( F, ,  FQ,  Fv,  F i ), 
F  is  ( 1 , 0, 0, 0).  Experiments  show  that  the  ellipticity 
of  light  reflected  from  the  natural  surfaces  is  suffi¬ 
ciently  small  to  be  neglected  (Coulson,  1978)2,  in 
which  case  Fv  -  V  =  0. 

The  first  term  of  Eq.  ( 1)  represents  the  solar  radia¬ 
tion  diffusely  reflected  from  the  top  of  the  atmos¬ 
phere.  The  second,  third,  fourth  and  fifth  terms 
represent  the  solar  radiation  that  is  directly  or  dif¬ 
fusely  transmitted  through  the  atmosphere,  reflected 
from  the  surface,  and  then  directly  or  diffusely 
transmitted  back  out  from  the  top  of  the  atmosphere. 
The  last  term  of  Eq.  (1)  describes  the  radiation  re¬ 
flected  more  than  once  from  the  surface.  For  clarity, 
the  factor  is  included  as  part  of  the  diffuse 

reflection  matrix  S  and  the  diffuse  transmission 
matrix  T  first  defined  by  Chandrasekhar  (1960). 
The  3x3  reflection  matrix  R  for  a  real  surface  is 
obtained  from  laboratory  measurements  of  the  bi¬ 
directional  reflectance,  degree  of  linear  polarization, 
and  angle  \  between  the  vertical  plane  and  the  plane 
of  polarization  of  the  reflected  light.  For  a  Lambert 
surface  all  the  elements  of  R  are  zero  except  Rn 


’  Coulson,  K.  L.,  1978:  Private  communications. 


which  is  a  function  of  the  total  reflectance  of  the 
surface.  Although  it  is  not  specifically  shown,  the 
product  of  two  bidirectional  matrices  also  implies 
hemispherical  integration.  For  example, 

S*R  =  I  I  S*( Tj ;  /a.  d>;  /a',  <t>’) 

Jo  Jo 

x  R(/x’,  d>';  /x0,  Mdn '<J>b ' . 

The  matrix  S*  is  related  to  the  diffuse  reflection 
matrix  S  by 

S*  =  QSQ, 

where 

/  1  0  0  \ 

Q  =  |  0  1  0  . 

\  0  0  - 1  / 

The  matrix  Q  takes  into  account  a  change  in  sym¬ 
metry,  discussed  by  Hovenier  (1%9),  which  occurs 
when  the  atmosphere  is  illuminated  from  the  bottom 
rather  than  from  the  top.  The  definition  of  T*  is 
similar. 

The  first  atmosphere  used  in  the  atmosphere- 
surface  model  was  a  Rayleigh  atmosphere.  The 
values  of  the  wavelength-dependent  t,(R)  used  were 
taken  from  Elterman  ( 1 968 )3  and  are  presented  in 
Table  1.  The  diffuse  reflection  and  transmission 
matrices  describing  the  pure  molecular  atmosphere 


3  Elterman.  L.,  I%8:  UV,  visible,  and  IR  attenuation  for  alti¬ 
tudes  to  50  km.  I%8.  AFCRL-68-0153.  49  pp. 
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were  evaluated  by  utilizing  a  computer  program 
developed  by  Dave  and  Warten  (1968). 4 

An  aerosol-laden  atmosphere  was  introduced  into 
the  atmosphere-surface  model  to  investigate  the 
effects  of  the  strong  forward  scattering  characteris¬ 
tic  of  aerosols.  The  total  normal  optical  thickness 
r,{  / )  of  the  turbid  atmosphere  is  the  sum  of  the 
Rayleigh  optical  thickness  r ,(/?)  and  the  normal  opti¬ 
cal  thickness  of  the  aerosol  content  T,(a).  The 
values  of  r,(«).  shown  in  Table  1,  were  taken  from 
Elterman  ( 1 968). 1  The  refractive  index  of  the  parti¬ 
cles  is  assumed  to  be  m  =  1.33  -  0/  with  a  size 
distribution  described  by  the  Haze  L  model  of 
Deirmendjian  ( 1969).  The  elements  of  the  diffuse 
reflection  and  transmission  matrices  describing  the 
turbid  atmosphere  were  calculated  using  the  adding 
method  with  a  program  graciously  provided  by  Dr. 
T.  Takashima  and  de>:ribed  in  Takashima  (1974). 
Because  of  the  method  oi  •:umerical  integration  used 
in  the  adding  method,  the  matrices  S  and  T  could  not 
be  calculated  at  the  exact  angles  tt0  used  to  illuminate 
the  real  surfaces  in  the  laboratory  experiments.  For 
surface  illumination  0„  =  53  and  78.5°,  the  differ¬ 
ence,  being  less  than  2.5°,  was  ignored.  For  ()„  =  0° 
the  difference  is  1 1°  but  its  influence  is  small  because 
S  and  T  as  a  function  of  0„  change  slowly  in  the  re¬ 
gion  l)„  =  0°. 

Once  the  diffuse  reflection  and  transmission 
matrices  describing  the  different  atmospheres  have 
been  computed,  only  the  surface  reflection  matrix, 
describing  the  natural  surfaces,  is  needed  in  order 
to  evaluate  the  terms  of  Eq.  (1)  for  the  different 
atmosphere-surface  combinations.  The  Stokes  param¬ 
eters  /,  Q  and  U.  together  with  the  degree  of  polari¬ 
zation.  were  computed  for  each  of  the  first  five  terms 
of  the  equation  by  using  the  Gaussian  quadrature 
method  of  integration. 

3.  Reflection  characteristics  of  the  natural  surfaces 

The  surfaces  used  in  this  study,  from  the  Lam¬ 
bert-like  gypsum  sand  to  the  highly  anisotropic- 
polarizing  Yolo  loam  (a  dark,  loamy  soil)  represent 
the  reflection  characteristics  of  a  variety  of  natural 
surfaces.  The  reflection  data  for  Yolo  loam,  gypsum 
sand  and  Mojave  Desert  sand  are  taken  from  the 
laboratory  measurements  by  Ghen  cl  til.  ( 1967). For 
these  measurements  the  surfaces  of  the  three 
samples  were  level.  The  bidirectional  reflectance 
and  polarization  of  the  reflected  light  were  measured 


1  Dave.  J.  V  and  K  M.  Warten.  |%8:  Program  for  computing 
l he  Stokes  parameters  of  the  radiation  emerging  from  a  plane 
parallel,  non-absorbing.  Rayleigh  atmosphere.  Rep  No.  *20 
*248.  IBM  Scientific  Center.  Palo  Alto.  (A.  80  pp 

Chen.  Hsi-shu.  (  R  Nagaraja  Rao  and  A  Sekera. 
Investigation  of  the  polarization  of  light  reflected  by  natural 
surfaces.  So.  Rep.  No  2.  Contract  AF  I9ffi28>- *850.  University 
of  California.  I  os  Angeles.  ^  pp 


Table  1.  Values  of  the  wavelength-dependent  optical  thick¬ 
nesses  used  to  construct  the  atmospheres  at  each  of  the  three 
wavelengths. 


Wavelength 
(m  m) 

r,  (*) 

Ti  (a) 

(1) 

0.4 

0.36 

0.32 

0.68 

0.5 

0.15 

0.26 

0.41 

0.6 

0.07 

0.24 

0.31 

for  three  types  of  illumination:  1)  natural  light,  2) 
light  polarized  90°  and  3)  light  polarized  45°  to  the 
vertical  plane.  The  use  of  the  different  types  of 
illumination  enabled  the  determination  of  nine  ele¬ 
ments  of  R  for  a  specific  surface,  wavelength  and 
angle  of  illumination.  The  angle  x  was  assumed  to 
be  perpendicular  to  the  scattering  plane;  this  as¬ 
sumption  was  found  by  observation  to  be  rea¬ 
sonable. 

The  reflected  radiances  measured  by  Chen  et  al. 
were  normalized  with  the  value  of  the  radiance  re¬ 
flected  from  the  given  surface  into  the  nadir  (when 
=  0°  normalization  was  performed  with  respect 
to  0  =  15°).  The  correction  factors  used  to  re¬ 
normalize  the  data  to  the  radiance  reflected  from  a 
standard  surface  of  known  reflectance  were  obtained 
from  Coulson  (1978)2  and  Walraven  and  Coulson 
(1972)6  who  measured  the  reflectance  of  the  same 
surface  samples  at  the  same  wavelengths.  For  Yolo 
loam  and  Mojave  Desert  sand  data  were  available 
for  all  three  types  of  illumination  except  for  ti„  =  0° 
where  only  data  for  illumination  by  natural  light 
were  available.  For  gypsum  sand  it  was  necessary 
to  renormalize  the  directional  reflectance  inde¬ 
pendent  of  illumination  type.  The  measured  polari¬ 
zation  of  the  light  reflected  from  the  surfaces  was 
not  affected  by  the  renormalization. 

The  bidirectional  reflectance  of  the  three  surfaces 
in  the  principal  plane  d>  =  0°(I80°)  for  the  wave¬ 
length  \  =  0.6  ^m  and  =  0“  is  shown  in  Fig.  2. 
The  reflectance  of  a  Lambert  surface  would  be  char¬ 
acterized  by  a  straight  horizontal  line  in  this  figure. 
The  reflectance  of  the  three  surfaces  is  independent 
of  the  azimuthal  angle,  so  that  the  hemispheric 
pattern  of  reflectance  can  be  generated  by  a  rotation 
around  the  nadir.  The  highly  reflecting  gypsum  sand 
shows  a  distinct  darkening  at  the  large  nadir  angles 
..nd  a  fairly  constant  value  of  reflectance  at  the 
nadir  angles  H  50°.  The  absorbent  Yolo  loam, 
in  contrast,  has  a  lower  reflectance  as  well  as  exhibit¬ 
ing  a  backward  maximum  centered  about  H  =  0°. 
Moreover.  Yolo  loam  shows  no  distinct  darkening 
at  large  nadir  angles  as  does  gypsum  sand.  The 
strong  backward  maximum  and  the  lack  of  darkening 

’  Walraven.  R.  I ...  and  k  1  (.Vu  Ivon.  W72:  Measurements  of 
Ihc  light  reflecting  properticv  of  white  gypsum  sand.  Contrih. 
ton. i\  Si  i  .  No  ~t.  I’nivcrvity  of  California.  Davis. 
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h  -  0°  4,  -  180“ 

Fic.  2.  Bidirectional  reflectance  of  gypsum  sand.  Mojave 
Desert  sand  and  Yolo  loam  in  the  principal  plane  for  A  =  0.6  /im 
and  t)0  =  0°. 

at  large  angles,  characteristic  of  Yolo  loam  but  not 
of  gypsum  sand,  is  probably  a  result  of  the  particles 
of  gypsum  sand  being  fairly  translucent  and  those  of 
Yolo  loam  opaque.  The  features  of  the  curve  char¬ 
acterizing  the  reflectance  of  Mojave  Desert  sand 
seem  to  be  a  blend  of  those  found  in  gypsum  sand 
and  in  Yolo  loam. 

The  effect  of  increasing  the  zenith  angle  of  the 
source  is  shown  in  Fig.  3.  The  measurements  are 
for  A  =  0.6  nm  and  t)0  =  53°.  The  curves  for  all 
three  surfaces  show  a  marked  increase  in  total  re¬ 
flectance  and  an  increase  in  the  magnitude  of  the 
backward  maximum.  The  strong  forward  reflection 
in  the  <|>  =  0°  plane  increases  with  increasing  9„. 


Fio.  3.  Bidirectional  reflectance  of  gypsum  sand.  Mojave 
Desert  sand  and  Yolo  loam  in  the  principal  plane  for  A  =  0.6  jim 
and  6„  =  53°,  and  of  Mojave  Desert  sand  for  A  *  0.4  ftm  and 
»a  =  53°. 


4>  =  0”  <1.  =  180' 

Fig.  4,  Degree  of  polarization  of  radiation  reflected  from 
gypsum  sand,  Mojave  Desert  sand  and  Yolo  loam  in  the  principal 
plane  for  A  =  0.6  /xm  and  H„  =  0°. 


Also  shown  in  the  figure  is  the  bidirectional  re¬ 
flectance  of  Mojave  Desert  sand  for  A  =  0.4  jim 
and  tt„  =  53°.  The  decrease  in  reflectance  with  de¬ 
creasing  wavelength  is  characteristic  of  all  three 
surfaces. 

The  degree  of  linear  polarization  of  the  radiation 
reflected  from  the  three  surfaces  in  the  principal 
plane  for  A  =  0.6  /am  and  0o  =  0°  is  presented  in 
Fig.  4.  A  comparison  of  the  profiles  in  Fig.  2  with 
those  of  Fig.  4  shows  that  the  reflectance  increases 
inversely  to  the  maximum  degree  of  polarization 
(Omov's  principle).  In  the  principal  plane,  positive 
polarization  indicates  that  the  plane  of  polarization 
is  perpendicular  to  the  plane:  negative  polarization 
indicates  that  the  plane  of  polarization  is  parallel 
to  this  vertical  plane.  The  positively  polarized  light 
is  considered  a  property  of  singly  reflected  light: 
the  negatively  polarized  light  occurring  in  the  re¬ 
gion  of  the  anti-source  is  considered  a  result  of 
light  being  doubly  reflected  (Wolff,  1975). 

Large  values  of  polarization  for  the  three  surfaces 
are  seen  in  Fig.  5  for  f)„  =  53°  and  A  =  0.4  /am.  The 
large  values  of  polarization  in  the  <t>  =  0°  plane  are 
probably  a  result  of  the  dominance  of  singly  re¬ 
flected  light  in  this  region.  The  maximum  degree  of 
polarization  for  Mojave  Desert  sand  is  larger  when 
the  wavelength  is  0.4  /am  than  when  it  is  0.6  /am  as 
shown  in  the  figure.  This  agrees  with  Omov's 
principle  since  the  reflectance  of  the  surface  in¬ 
creases  as  the  wavelength  increases. 

4.  Radiance  outward  from  a  Rayleigh  atmosphere 

For  these  computations,  each  real  surface  was 
assumed  to  be  an  infinite  plane  underlying  a  plane 
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Fk..  5.  Degree  of  polarization  of  radiation  reflected  from 
gypsum  sand,  Mojave  Desert  sand  and  Yolo  loam  in  the  principal 
plane  for  A  0.4pmandl)„  -  s.T,  and  for  Mojave  Desert  sand  for 
A  0.6  nm  and  =  53”. 


90  60  30  0  30  60  90 

NADIR  ANGLE  (Degrees) 

d>  =  0°  d>  =  180° 

Fit,  6  Relative  radiance  of  the  major  surface  reflection  terms 
of  the  total  emergent  radiation  as  defined  by  Eq.  (1)  for  a  pure 
molecular  atmosphere  overlying  Mojave  Desert  sand  (principal 
plane,  A  =  0.4  pm,  N  -  53°).  Solid  lines  indicate  real  surface 
reflection;  dashed  lines  indicate  Lambert  surface  reflection. 


parallel  pure  molecular  model  of  the  earth's 
atmosphere. 

a.  Result s  for  a  Mojave  Desert  sand  surface 

The  relative  radiances  of  the  first  three  surface 
reflection  terms  of  Eq.  ( 1 1  for  the  light  emerging  from 
the  top  of  a  Rayleigh  atmosphere  overlying  Mojave 
Desert  sand  for  A  =  0.4  pirn  and  H0  =  53°  are  shown 
by  the  profiles  in  Fig.  6.  The  radiance  of  the  second 
term  of  Eq.  ( 1)  as  a  function  of  nadir  angle  for  real 
surface  reflection  (solid  line)  is  compared  to  the  radi¬ 
ance  of  the  same  term  for  reflection  by  a  Lambert 
surface  (dashed  line)  of  the  same  total  mono¬ 
chromatic  reflectance,  A  ^  0.12.  At  large  nadir 
angles  the  differences  between  the  two  curves  show 
the  effect  of  the  forward  and  backscattered  peaks 
in  the  reflected  radiance  for  the  desert  sand  de¬ 
picted  in  Fig.  3.  Though  not  as  great,  these 
differences  are  still  present  for  the  third  term  which 
represents  diffuse  illumination  of  the  surface.  The 
curves  of  both  the  second  and  third  terms  show  that 
the  L.ambert  approximation  overestimates  the  real 
surface  reflection  in  the  region  of  the  nadir.  The 
strong  smoothing  influence  of  diffuse  transmission 
out  from  the  top  of  the  atmosphere  is  illustrated  by 
the  lack  of  any  observable  difference  between  the 
fourth  term  of  Eq.  (1)  and  its  Lambert  counterpart. 
The  differences  between  the  remaining  surface  re¬ 
flection  terms  of  Eq.  ( 1 1  and  their  Lambert  counter¬ 
parts  are  negligible. 

The  influence  on  the  emergent  radiance  of  reflec¬ 
tion  by  Mojave  Desert  sand  instead  of  its  Lambert 
equivalent  is  shown  in  Fig.  7.  The  figure  shows  the 
percent  change  or  needed  correction  in  the  total 


NADIR  ANGLE  (Degrees) 

<6  =  0°  <6  -  180° 


Fi<;.  7.  Correction  to  the  total  radiance  emerging  from  a 
Rayleigh  atmosphere  overlying  a  Lambert  surface  of  total  re¬ 
flectance  A  that  occurs  when  the  reflection  matrix  of  Mojave 
Desert  sand  is  used  in  term  two  of  Eq.  (I). 
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Fig.  8.  Correction  >o  the  total  radiance  emerging  from  a 
Rayleigh  atmosphere  overlying  a  Lambert  surface  of  total  re¬ 
flectance  A  that  occurs  when  the  reflection  matrix  of  gypsum 
sand  is  used  in  term  two  of  Eq.  ( I)  or  in  terms  two  and  three  as 
shown  by  the  dashed  curve. 

emergent  radiance  calculated  using  the  Lambert 
approximation  in  Eq.  (1)  as  compared  to  the  results 
obtained  when  real  surface  reflection  is  included  in 
the  dominant  second  term  and  the  Lambert  ap¬ 
proximation  in  the  remaining  surface  reflection 
terms.  A  positive  correction  indicates  that  Lambert 
surface  reflection  underestimates  the  real  surface 
reflection  and  vice  versa.  The  positive  values  of  cor¬ 
rection  in  the  region  of  the  nadir  and  the  negative 
values  of  correction  at  large  0  for  A  =  0.4  and  0.6 
fxm  when  0„  =  0°  show  the  influence  of  the  back- 
scatter  maximum  and  the  decrease  in  bidirectional 
reflectance  with  increasing  0  for  the  desert  sand  seen 
in  Fig.  2.  The  small  values  of  correction  for  the 
wavelength  0.4  jtm  are  a  result  of  the  large  optical 
thickness  at  this  wavelength  which  decreases  the 
contribution  of  the  surface  reflection  terms  to  the 
total  emergent  light.  The  increasing  importance  of 
the  forward  scattered  peak  of  the  bidirectional 
reflectance  as  the  solar  zenith  angle  increases  is 
also  shown  in  Fig.  7.  When  A  =  0.5  /xm  the  maxi¬ 
mum  change  in  the  <f>  =  0°  plane  is  —14%  for  0„ 
=  78.5°  and  7%  for  0O  =  53°.  As  the  wavelength  in¬ 
creases  to  0.6  pm,  decreasing  the  influence  of  the 


first  term  of  Eq.  (1),  the  effects  of  the  real  surface 
become  very  large  for  0O  =  78.5°.  It  is  interesting 
that  the  Lambert  surface  approximation  works  well 
in  the  region  6  =  40-45°  and  <f>  =  0°  for  all  the 
cases  studied  in  Fig.  7. 

The  third  term  of  Eq.  (I)  is  of  next  importance, 
after  the  second  term,  in  introducing  the  effects 
of  real  surface  reflection.  The  influence  on  the  cor¬ 
rection  of  the  inclusion  of  real  surface  reflection  by 
the  third  term  which  represents  diffuse  illumination 
of  the  surface  is  largest  for  the  larger  solar  zenith 
angles.  This  term  has  its  maximum  influence  in  the 
region  of  the  nadir  where,  for  A  =  0.5  /xm,  it  con¬ 
tributes  -1.4  and  -2.6%  to  the  correction  when 
0O  =  53  and  78.5°,  respectively.  When  0O  =  0°  the 
cumulative  influence  of  the  surface  reflection  terms, 
other  than  the  second,  on  introducing  the  effects  of 
real  surface  reflection  is  negligible. 

b.  Results  for  a  gypsum  sand  surface 

The  correction  introduced  by  gypsum  sand  is  pre¬ 
sented  in  Fig.  8.  The  change  introduced  by  the 
second  term  for  A  =  0.6  pm  and  0„  =  78.5°  again 
shows  the  importance  of  the  forward  and  back- 
scattered  maxima  on  the  total  emerging  light.  When 
A  =  0.4  pm,  the  change  is  much  smaller  due  to  a 
decrease  in  the  surface  reflection  and  an  increase  in 
the  optical  thickness  of  the  atmosphere.  The  cumula¬ 
tive  effect  of  real  surface  reflection  by  the  second 
and  third  terms,  shown  by  the  dashed  curve,  shows 
that  in  the  region  of  the  nadir  the  effect  of  these  two 
terms  is  fairly  equal  and  their  total  effect  on  the 
emergent  radiance  is  significant.  The  influence  of  the 
darkening  at  large  0  and  the  backscattered  maximum 
present  in  the  curves  of  bidirectional  reflectance 
for  gypsum  sand  is  clearly  visible  in  the  curve  for 
A  =  0.6  and  0o  =  0°.  For  all  the  cases  shown  in 
Fig.  8  the  Lambert  approximation  works  well  in 
the  region  of  0  =  40°. 

c.  Results  for  a  Yolo  loam  surface 

To  investigate  the  effects  of  a  surface  with  a  low 
total  reflectance,  a  Yolo  loam  surface  was  included 
in  the  atmosphere-surface  model.  Fig.  9  presents 
the  curves  of  correction  which  demonstrate  the 
influence  of  real  surface  reflection  by  the  second 
term  only.  The  pronounced  backscattered  maximum 
characteristic  of  Yolo  loam  is  clearly  shown  by  the 
positive  values  of  correction  in  the  region  of  the 
nadir  for  A  =  0.4  and  0.6  pm  when  0O  =  0°.  A  change 
in  solar  zenith  angle  to  78.5°  for  A  =  0.6  /xm  pro¬ 
duces  large  corrections  in  the  region  of  the  nadir  and 
the  forward  and  backscatter  maxima.  These  correc¬ 
tions  are  effectively  reduced  by  the  larger  at¬ 
mospheric  reflection  and  smaller  surface  reflection 
present  when  A  =  0.4  /xm.  Including  the  effect  of 
real  reflection  by  the  third  term,  however,  approxi¬ 
mately  doubles  the  correction  in  the  region  of  the 
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<l>  -  0°  </>  =  180” 

Fig.  9.  Correction  to  the  total  radiance  emerging  from  a 
Rayleigh  atmosphere  overlying  a  Lambert  surface  of  total  re¬ 
flectance  A  that  occurs  when  the  reflection  matrix  of  Yolo  loam 
is  used  in  term  two  of  Eq.  (I). 

nadir.  The  Lambert  approximation  works  well  in  the 
region  of  0  ~  45°  for  the  cases  presented  in  Fig.  9. 

The  previous  discussions  were  confined  to  the 
<f>  =  0°  (180°)  plane.  The  difference  in  radiance  of 


the  second  term  for  reflection  by  Mojave  Desert 
sand  and  its  Lambert  equivalent  is  presented  in  the 
hemisphere  map  (Fig.  10).  The  source  zenith  angle 
is  53°  and  the  wavelength  is  0.5  /im.  The  large 
differences  are  centered  on  the  principal  plane, 
indicating  the  importance  of  this  plane.  The  positive 
differences  at  large  nadir  angles  are  a  result  of  the 
forward  and  backscattered  maxima.  The  negative 
values  in  the  region  of  the  nadir  indicate  the  Lambert 
surface  approximation  overestimates  the  reflected 
radiance  in  this  region.  Because  the  surface  re¬ 
flection  characteristics  are  symmetrical  about  the 
plane  <t>  =  0°  (180°)  the  values  in  Fig.  10  are  repre¬ 
sentative  of  the  values  occurring  in  the  opposite 
hemisphere. 

5.  Radiance  outward  from  a  turbid  atmosphere 

A  turbid  atmosphere  was  used  in  the  atmosphere- 
surface  model  with  Mojave  Desert  sand  to  study  the 
effects  of  aerosol  scattering  on  the  emergent  light. 
The  effect  of  real  surface  reflection  by  the  second 
term  on  the  total  radiance  is  shown  for  both  the 
turbid  and  pure  molecular  atmospheres  in  Fig.  1 1 . 
The  decrease  in  correction  caused  by  the  increased 
multiple  scattering  of  the  turbid  atmosphere  is  evi¬ 
dent  by  comparing  curves  representing  the  different 
atmospheres  but  the  same  wavelength  and  solar 
zenith  angle.  In  addition,  the  figure  compares  the 
effect  of  the  second  and  third  terms.  The  combined 
influence  of  real  surface  reflection  by  the  second  and 
third  terms  for  \  =  0.6  /j.m  and  0O  =  78.5°  is  shown 
by  the  dashed  line  in  Fig.  II.  The  effect  of  the  third 
term,  depicted  by  the  large  differences  between  the 
dashed  line  and  the  line  symbolized  by  plus  marks 
indicate  that  the  strong  forward  scattering  of  the 
aerosols  have  significantly  offset  the  smoothing 
effect  of  multiple  scattering  by  the  aerosols.  The 
influence  of  the  third  term  is  also  significant  at 


Fig.  10.  Hemispheric  map  of  the  difference  in  relative  radiance  of  the  second 
term  of  Eq.  (I)  for  real  and  Lambert  surface  reflection  (Mojave  Desert  sand, 
X  =  0.5  |xm.  #0  =  53°). 
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Flo.  11.  Correction  to  the  total  radiance  emerging  from  Ray¬ 
leigh  and  turbid  atmospheres  overlying  a  Lambert  surface  of  total 
reflectance  A  that  occurs  when  the  reflection  matrix  of  Mojave 
Desert  sand  is  used  in  term  two  of  Eq.  ( 1 )  or  in  terms  two  and  three 
as  shown  by  the  dashed  curve. 


K  =  0.4  /Lin,  though  not  shown  in  the  figure.  When 
0O  =  78.5°  and  X.  =  0.4  /xm  the  maximum  correction 
in  the  region  of  the  nadir  is  increased  from  less  than 
0.3%  to  slightly  over  2%  when  the  effect  of  this 
term  is  included.  When  0„  =  0°  the  cumulative  in¬ 
fluence  of  the  surface  reflection  terms,  other  than 
the  second,  is  negligible,  except  at  A  =  0.6  /im.  In 
this  case,  the  second  and  fourth  terms  can  combine 
at  large  0  to  give  a  total  correction  of  -3%.  In 
general,  the  corrections  occurring  for  the  turbid 
atmosphere  increase  with  increasing  wavelength 
but  remain  below  those  found  in  the  pure  molecular 
atmosphere.  The  Lambert  surface  approximation 
works  well  in  the  region  of  0  =  45°  for  the  Rayleigh 
anti  turbid  atmospheres  presented  in  Fig.  II. 

Including  the  effects  of  aerosol  scattering  in  the 
atmosphere-surface  model  tends  to  decrease  the  im¬ 
portance  of  using  a  real  surface  for  calculating  the 
emergent  radiance.  Nevertheless,  the  influence  of 
real  surface  reflection  for  an  aerosol-laden  atmos¬ 
phere  is  still  significant.  The  apparent  importance  of 
the  fourth  term  at  large  nadir  angles  when  0O  =  0°  is 


probably  the  result  of  a  relatively  large  direct 
transmission  of  light  to  the  surface,  followed  by  a 
large  diffuse  transmission  out  from  the  top  of  the 
atmosphere  at  large  nadir  angles.  When  0O  =  78.5° 
the  importance  of  the  third  term  at  small  nadir 
angles  is  likely  due  to  a  highly  anisotropic  beam  of 
diffuse  light  incident  on  the  surface  which,  after 
surface  reflection,  is  directly  transmitted  out  through 
the  atmosphere  at  the  smaller  optical  thickness 
occurring  in  the  region  of  the  nadir.  The  strong 
forward  scattering  of  the  aerosols  preserves  some  of 
the  anisotropy  of  the  surface  reflected  light  even  in 
an  optically  thick  atmosphere  with  large  multiple 
scattering. 


6.  Polarization  outward  from  a  Rayleigh  atmosphere 

Another  interesting  feature  of  reflection  by  natural 
surfaces  is  their  effect  on  the  behavior  of  the  polari¬ 
zation  of  the  emergent  light. 

a.  Results  for  a  Mojave  Desert  sand  surface 

The  degree  of  polarization  of  the  emergent  radia¬ 
tion  in  the  principal  plane  for  Mojave  Desert  sand 
is  shown  in  Fig.  12  for  \  =  0.5  /xm  and  0O  =  53°. 
Included  for  comparison  are  the  curves  of 
polarization  for  a  Rayleigh  atmosphere  overlying  a 
Lambert  surface  of  equivalent  total  reflectance  and 
the  polarization  of  the  light  reflected  from  the  desert 
sand  surface  itself.  The  maximum  polarization  of  the 
outward  directed  light  for  the  sand  surface  is  —8% 
more  than  that  for  the  Lambert  surface,  and  is  shifted 
~5°  toward  the  nadir.  There  is  no  apparent  shift  in 
the  neutral  points  (points  of  zero  polarization)  which 
may  be  a  result  of  the  small  change  in  polarization 
in  this  region  for  the  sand-reflected  light. 

The  difference  between  the  degree  of  linear 
polarization  emerging  from  the  top  of  a  Rayleigh 
atmosphere  overlying  Mojave  Desert  sand  and  a 
Lambert  surface  of  the  same  total  reflectance  as  a 
function  of  nadir  angle  is  shown  in  Fig.  13  for 
several  wavelength  and  solar  zenith  angle  combina¬ 
tions.  The  increasing  values  of  the  difference  with 
increasing  0  for  \  =  0.6  /xm  and  0„  =  0°  agree  with 
the  trend  in  polarization  shown  in  Fig.  4  for  the 
desert  sand.  In  a  similar  manner,  the  other  curves 
in  the  figure  follow  the  trend  of  the  measured  degree 
of  polarization  of  light  reflected  from  the  desert 
sand.  To  place  the  differences  of  Fig.  13  in  perspec¬ 
tive  for  sensing  from  space,  ground-based  polarizing 
radiometers  have  been  constructed  that  measure  the 
degree  of  polarization  to  ±0.1%,  but  a  satellite- 
borne  instrument  will  probably  measure  the  degree 
of  polarization  to  ±  1.0%  (Walraven,  1980). 

The  section  on  emergent  radiance  showed  that  the 
second  term  of  Eq.  (1)  was  dominant  in  transmitting 
the  characteristics  of  the  natural  surfaces.  The 
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second  term  and  any  term  representing  direct 
transmission  out  from  the  top  of  the  atmosphere  will 
not  be  polarized  when  reflection  is  from  a  Lambert 
surface.  Approximating  the  polarization  of  the 
emergent  light  using  the  total  emergent  radiance  for 
real  surface  reflection  but  only  using  the  polariza¬ 
tion  introduced  by  the  first  and  second  terms  gives 
close  agreement  with  the  results  obtained  when  the 
polarization  introduced  by  all  the  terms  of  Eq.  (I) 
is  used.  For  the  three  wavelengths  the  maximum 
difference  in  degree  of  polarization  between  the 
approximation  and  the  rigorous  case  for  Mojave 
Desert  sand  when  tf0  =  0°  occurs  at  large  0  and  is 
1.4'; .  If  the  polarization  of  the  next  important 
term,  the  fourth,  is  included  in  the  approximation. 


Fit.  12  A  comparison  of  (he  degree  of  polarization  of  the 
light  emerging  from  the  top  of  a  Rayleigh  atmosphere  overlying 
a  Mojave  Desert  sand  surface  and  its  Lambert  equivalent.  Also 
shown  is  the  degree  of  polarization  of  light  reflected  from  the 
desert  sand  < principal  plane,  K  =  0.5  pm.  H„  =  53"). 


90  60  30  0  30  60  90 

NADIR  ANGLE  (Degrees) 


<*.  =  0° 


<t>~  180° 


FtC.  13.  The  difference  in  the  degree  of  polarization  of  the 
radiation  emerging  from  the  lop  of  a  Rayleigh  atmosphere  that 
occurs  for  reflection  by  Mojave  Desert  sand  and  its  Lambert 
equivalent. 


the  maximum  difference,  as  shown  in  column  II 
of  Table  2,  is  reduced  to  0.2%.  For  solar  zenith 
angles  of  53  and  78.5°  the  maximum  difference 
between  the  approximation  and  the  rigorous  case  is 
reduced  to  a  maximum  of  0.3%  if  the  third  term  of 
Eq.  (1)  is  included  in  the  approximation. 

b.  Results  for  a  gypsum  sand  surface 

Gypsum  or  white  quartz  sand  is  a  highly  reflect¬ 
ing  surface,  but  the  light  reflected  from  the  sand  is 
not  highly  polarized.  However,  as  shown  in  Fig.  14, 
there  are  significant  departures  of  the  degree  of 
polarization  from  the  Lambert  case  when  the  effects 
of  the  real  polarizing  surface  are  taken  into  account. 
For  a  wavelength  of  0.4  and  0.6  jim  when  0O  =  0° 
the  differences  in  degree  of  polarization  increase 
with  increasing  values  of  nadir  angle,  except  for 
0  >  75°  and  A  =  0.4  gm  where  the  large  optical 
thicknesses  decrease  the  effect  of  the  real  surface. 
The  importance  of  the  fourth  term  is  demonstrated 
by  the  decrease  of  the  maximum  difference  from 
1.8  to  0.2%.  as  shown  in  Table  2,  when  real  surface 
reflection  by  this  term  is  included  in  the  approxima- 
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Table  2.  Hie  maximum  difference  (I)  between  the  approximated  and  rigorously  calculated  values  of  polarization  for  the  different 
atmosphere- surf ace  models,  and  the  maximum  difference  (II)  when  the  next  important  term  of  Eq.  (1)  is  included  in  the  spproximation. 


(I) 

(ID 

Maximum 

Maximum 

difference 

difference 

Added 

Atmosphere 

Surface 

to  (deg) 

(%) 

(%) 

term 

Rayleigh 

desert  sand 

0 

1.4 

0.2 

fourth 

53,  78.5 

0.8 

0.3 

third 

gypsum  sand 

0 

1.8 

0.2 

fourth 

78.5 

0.5 

0.2 

third 

Yolo  loam 

0 

0.8 

0.05 

fourth 

78.5 

0.5 

0.3 

third 

Turbid 

desert  sand 

0 

1.2 

0.5 

fourth 

78.5 

1.4 

0.3 

third 

tion.  The  influence  of  large  optical  depths  is  seen  by 
comparing  the  curves  in  the  <t>  =  0°  plane  for 
A  =  0.4  /rm  and  A  =  0.6  /tm  when  00  =  78.5°.  The 
increased  multiple  scattering  occurring  for  A  =  0.4 
/zm,  compared  to  A  =  0.6  pirn,  masks  the  contribu¬ 
tion  of  the  real  surface  to  the  total  polarization.  For 
both  wavelengths,  however,  the  approximation 
works  well,  as  shown  in  Table  2,  showing  the 
dominance  of  the  first  and  second  terms  of  Eq.  (1) 
in  contributing  to  the  polarization. 

c.  Results  for  a  Yolo  loam  surface 

The  polarization  emerging  from  a  Rayleigh  at¬ 
mosphere  overlying  the  dark,  but  highly  polarizing, 
Yolo  loam  surface  is  shown  in  Fig.  15  for  A  =  0.6 
ftm  and  0O  =  0°.  The  higher  values  of  polarization 
for  the  atmosphere-loam  model  at  0  >  30°  compared 
to  the  atmosphere-Lambert  model  reflect  the  polari¬ 
zation  characteristics  of  the  loam  in  this  region.  In 
the  region  of  the  nadir,  the  negatively  polarized 
light  reflected  by  the  surface  introduces  a  neutral 
point  for  the  atmosphere-loam  model,  whereas 
there  was  apparently  none  for  the  atmosphere- 
Lambert  model. 

The  departures  from  the  Lambert  case  caused  by 
the  highly  polarizing  Yolo  loam  surface  are  shown  in 
Fig.  16  for  a  number  of  wavelength  and  solar  zenith 
angle  combinations.  The  shape  and  orientation  of  the 
curves  with  respect  to  one  another  are  similar  to 
the  case  of  reflection  by  Mojave  Desert  sand  (Fig. 
13).  However,  unlike  the  other  surfaces,  Yolo  loam 
shows  a  large  negative  difference  in  the  region  of  the 
nadir  for  0O  =  0°.  This  is  due  to  the  presence  of  large 
values  of  negative  polarization  of  the  reflected  light 
in  this  region,  as  shown  in  Fig.  4.  The  approxima¬ 
tion  again,  as  shown  in  Table  2,  works  well  in  esti¬ 
mating  the  polarization  of  the  light  outward  from  the 
atmosphere  overlying  the  loam. 

7.  Polarization  outward  from  a  turbid  atmosphere 

The  polarization  of  the  total  radiance  emerging 
from  a  turbid  atmosphere  is  generally  smaller  than 


for  a  Rayleigh  atmosphere  as  shown  in  Fig.  17  for 
the  Mojave  Desert  sand  surface.  The  wavelength  is 
0.6  /im  and  the  zenith  angle  of  the  sun  is  78.5°.  Also 
compared  in  the  figure  are  the  polarization  curves  of 
the  light  emerging  from  the  turbid  atmosphere  over- 
lying  the  desert  sand  and  its  equivalent  Lambert 
surface.  Both  curves  show  the  existence  of  the  neu¬ 
tral  point  (Arago)  at  <f>  =  0°.  The  point  is  predicted 
for  Rayleigh  scattering  alone  when  values  of  0„  are 
large  but  its  occurrence  was  probably  suppressed 


Fig.  14.  The  difference  in  the  degree  of  polarization  of  the 
radiation  emerging  from  the  top  of  a  Rayleigh  atmosphere  that 
occurs  for  reflection  by  gypsum  sand  and  its  Lambert  equivalent. 
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because  of  the  large  contribution  of  positive  polari¬ 
zation  by  the  desert  sand  in  the  d>  =  0°  plane.  It  is 
interesting  that  the  differences  in  polarization  be¬ 
tween  the  curves  for  the  turbid  and  Rayleigh  atmos¬ 
pheres  overlying  desert  sand  can  be  as  large  as  the 
differences  occurring  between  the  curves  for  the 
turbid  atmosphere  overlying  the  sand  and  the  Lam¬ 
bert  surface. 

The  difference  in  degree  of  polarization  as  a  func¬ 
tion  of  nadir  angle  between  real  surface  and  Lambert 
surface  reflection  is  presented  in  Fig.  18  for  turbid 
and  Rayleigh  atmospheres  overlying  Mojave  Desert 
sand.  For  angles  «70°  and  0o  =  0°  the  differences  in 
polarization  for  the  turbid  atmosphere  are  close  to 
those  for  the  molecular  atmosphere.  At  larger  nadir 
angles  the  two  cases  diverge  from  one  another,  a 
result  of  an  increase  in  multiple  scattering.  At  a 
solar  zenith  angle  of  78.5°  the  effect  of  the  turbid 
atmosphere  in  reducing  the  departures  from  Lambert 


90  60  30  0  30  60  90 

NADIR  ANGLE  (Degrees) 

<t>  -  0°  <b-  180° 

Fig.  15.  A  comparison  of  the  degree  of  polarization  of  the 
light  emerging  from  the  top  of  a  Rayleigh  atmosphere  overlying 
a  Yolo  loam  surface  and  its  Lambert  equivalent.  Also  shown  is 
the  degree  of  polarization  of  light  reflected  from  the  loam 
(principal  plane,  A  =  0.6  jim.  9„  -  0°). 


Fig.  16.  The  difference  in  the  degree  of  polarization  of  the 
radiation  emerging  from  the  top  of  a  Rayleigh  atmosphere  that 
occurs  for  reflection  by  Yolo  loam  and  its  Lambert  equivalent. 


are  large  compared  to  the  Rayleigh  atmosphere.  The 
departures  from  Lambert  for  the  turbid  atmosphere, 
however,  are  still  significantly  large.  The  increased 
importance  of  the  terms  of  Eq.  (1)  with  histories  of 
diffuse  transmission  are  evident  by  the  large  differ¬ 
ences  occurring  (shown  in  Table  2)  for  a  turbid 
atmosphere  overlying  desert  sand.  With  the  addition 
of  either  the  third  or  fourth  terms,  however,  the  ap¬ 
proximation  of  the  polarization  is  greatly  improved. 

8.  Conclusions 

The  surfaces  used  in  this  study  represent  a  variety 
of  reflection  characteristics  from  the  Lambert-like 
gypsum  sand  to  the  highly  anisotropic  polarizing 
Yolo  loam.  The  calculations  support  the  following 
conclusions: 

1)  The  radiance  emerging  from  the  lop  of  a 
Rayleigh  atmosphere  can  be  significantly  different 
fora  real  surface  compared  to  its  Lambert  equivalent 
and  the  difference,  in  general,  increases  as  the 
solar  zenith  angle  and  the  wavelength  increase. 

2)  The  surface  reflection  characteristics  are  car: 
ried  mainly  by  the  second  term  of  Eq.  (1).  This  fact 
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Fic.  17.  A  comparison  of  (he  degree  of  polarization  of  the 
light  merging  from  the  top  of  the  turbid  atmosphere  plus  desert 
sand,  Kayleigh  atmosphere  plus  desert  sand,  and  turbid  atmos¬ 
phere  plus  Lambert  surface  models.  Also  shown  is  the  degree  of 
polarization  of  light  reflected  from  the  desert  sand  (principal 
plane,  k  =  0.6  *im,  6a  =  78.5°). 


has  been  used  by  Coulson  et  al.  (1966)  and  Tanre 
et  ul.  (1979)  to  simplify  their  calculations. 

3)  The  third  term  is  of  next  importance;  its  influ¬ 
ence  in  the  region  of  the  nadir  can  exceed  that  of  the 
second  term  at  the  shorter  wavelengths  and  larger 
solar  zenith  angles. 

4)  The  increased  multiple  scattering  of  the  turbid 
atmosphere  decreases  the  importance  of  including 
the  characteristics  of  real  surface  reflection  in  cal¬ 
culating  the  emergent  light.  The  maximum  correc¬ 
tion,  however,  may  reach  3.0%  for  A  =  0.4  /am  and 
0O  =  0°  or  78.5°  and  more  than  10%  for  A  =  0.6  /am 
and  0O  =  78.5°. 

5)  The  Lambert  approximation  works  well  in  the 
range  0  =  40-45°,  <f>  =  0°  for  allpf  the  atmosphere- 
surface  combinations  studied.  It  may  be  possible  to 
use  the  radiance  measured  in  this  range  to  normalize 
the  entire  radiance  field  producing  a  representation 


of  the  non-Lambertian  reflection  characteristics  of 
the  surface. 

6)  Large  differences  can  occur  between  the  degree 
of  polarization  of  radiation  emerging  from  the  top 
of  a  pure  molecular  atmosphere  overlying  a  natural 
as  compared  to  a  Lambert  surface.  The  profiles  of 
these  differences  are  directly  related  to  the  curves 
of  polarization  of  the  light  reflected  by  the  surfaces. 
As  was  the  case  with  radiance,  the  departures  of  the 
polarization  from  the  Lambert  case  are  largest  for 
the  larger  solar  zenith  angles  and  wavelengths. 

7)  The  degree  of  polarization  of  the  light  emerging 
from  the  top  of  the  Rayleigh  atmosphere  overlying 
a  real  surface  can  be  approximated  to  an  uncertainty 
of  0.3%  or  less  by  including  real  surface  reflection 
in  only  the  second  and  either  the  third  or  fourth 
terms  of  Eq.  (1). 

8)  Compared  to  a  Rayleigh  atmosphere  the  differ¬ 
ences  in  polarization  are  smaller  for  the  turbid 
atmosphere,  but  they  are  still  significant.  Again,  the 
polarization  of  the  emergent  light  can  be  accurately 
approximated,  in  this  case  to  0.5%  or  less,  by  using 
only  the  second  and  either  the  third  or  fourth  terms 
of  Eq.  (1). 


Fig.  18.  The  difference  in  the  degree  of  polarization  of  the 
radiation  emerging  from  the  top  of  Rayleigh  and  turbid  atmos¬ 
pheres  that  occurs  for  reflection  by  Mojave  Desert  sand  and  its 
Lambert  equivalent. 
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There  are  surprisingly  large  departures  in  the 
polarization  of  radiation  emerging  from  the  top  of 
the  Rayleigh  and  turbid  atmospheres  for  a  real  as 
compared  to  a  Lambert  surface.  These  departures 
show  that  there  is  a  significant  amount  of  infor¬ 
mation  about  the  real  polarizing  surface  in  the 
polarization  of  the  emergent  light.  Polarization  ap¬ 
pears  to  be  a  very  good  parameter  for  obtaining 
additional  information  about  the  features  of  a  natural 
ground  surface. 

This  work  needs  to  be  extended  into  several  areas. 
One  area  is  the  study  of  how  much  large-scale 
topographic  relief  will  reduce  the  effect  of  real  sur¬ 
face  reflection  on  the  emergent  radiation.  In  addi¬ 
tion,  the  study  of  the  effects  of  a  surface  composed 
of  two  or  more  distinct  surface  types  is  of  interest. 

Acknowledgments.  It  is  mv  pleasure  to  thank  Prof. 
Kinsell  L.  Coulson  for  then  any  constructive  discus¬ 
sions  we  had.  1  also  wish  to  thank  Dr.  Robert 
Walraven  for  his  helpful  comments  and  suggestions. 
Dr.  Tsutomu  Takashima  generously  made  his  pro¬ 
gram  available  to  model  the  inhomogeneous  atmos¬ 
phere.  The  phase  functions  describing  the  aerosol 
were  graciously  provided  by  Dr.  Robert  Fraser.  The 
aid  of  Mr.  Keith  Grant  is  gratefully  acknowledged. 
Computer  time  was  kindly  provided  by  the  Atmos¬ 
pheric  and  Geophysics  Division  of  Lawrence  Liver¬ 


more  Laboratories.  The  support  for  this  investiga¬ 
tion,  which  was  provided  by  the  National  Science 
Foundation  under  Grants  DV40893,  DPP76-22260 
and  DPP77- 19362,  and  by  the  Air  Force  Geophysics 
Laboratory  under  Contract  F19628-78-C-0200,  is 
gratefully  acknowledged. 


REFERENCES 

Chandrasekhar,  S.,  1960:  Radiative  Transfer.  Dover,  393  pp. 

Coulson,  K.  L..  E.  L.  Gray  and  G.  M.  B.  Bouricius,  1966:  Effect 
of  surface  reflection  on  planetary  albedo.  Icarus ,  5, 1 39- 148. 

Deirmendjian,  D.,  1969:  Electromagnetic  Scattering  on  Spherical 
Polydispersions.  Elsevier,  290  pp. 

Hovenier,  J.  W.,  1969:  Symmetry  relationships  for  scattering  of 
polarized  light  in  a  slab  of  randomly  oriented  particles. 
J.  Atmos.  Sci.,  26,  488-499. 

Koepke,  P  ,  and  K.  T.  Kriebel,  1978:  Influence  of  measured 
reflection  properties  of  vegetated  surfaces  on  atmospheric 
radiance  and  its  polarization.  Appl.  Opt.,  17,  260-264. 

Takashima,  T..  1974:  Method  of  computing  the  effect  of  surface 
reflection  on  the  atmospheric  radiation.  Puhl.  Astron.  Soc. 
Japan.  26  ,  361-366. 

Tanre,  D.,  M.  Herman,  P.  Y.  Deschamps  and  A.  de  Leffe,  1979: 
Atmospheric  modeling  for  space  measurements  of  ground 
reflectances,  including  bidirectional  properties.  Appl.  Opt., 
18,  3587-3594 

Walraven.  R.  L.,  1980:  Polarization  imagery.  Opt.  Eng.,  20, 
14-18. 

Wolff.  M.,  1975:  Polarization  of  light  reflected  from  rough  plane¬ 
tary  surface.  Appl.  Opt..  14,  1395-  1405. 


Reprinted  from  Jovxnal  or  AfpUB)  Meteoeolocy,  Vol.  20,  No.  10,  October  1981 
American  Meteorological  Society 
Printed  in  U.  S.  A. 


APPENDIX  D 


Measurements  of  Aerosol  Size  Distributions  in  the 
Lower  Troposphere  over  Northern  Europe 

Bruce  W.  Fitch 

Visibility  Laboratory,  Scripps  Institution  of  Oceanography,  University  of  California,  San  Diego,  92093 

Ted  S.  Cress' 

Air  Force  Geophysics  Laboratory,  Hanscom  AFB,  MA  01731 
(Manuscript  received  II  November  1980,  in  final  form  23  June  1981) 

ABSTRACT 

Airborne  measurements  of  particle  size  distributions  were  made  at  several  altitudes  within  and  above 
the  mixing  layer  at  sites  near  Ahlhom  and  Meppen  (West  Germany),  Rodby  (Denmark),  and  Bruz  (France). 
The  distributions  were  measured  over  the  range  0.2- 3. 9  jtm  in  particle  radius  using  a  Royco  model 
220  particle  counter.  The  experimental  data  gathered  at  the  sites  were  analyzed  in  terms  of  volume-size 
distributions  instead  of  the  commonly  used  particle  number  distributions.  The  volume  distribution  plots 
were  found  to  be  bimodal  with  the  accumulation  mode  centered  at  a  mean  radius  in  the  range  0.26-0.49 
fim  and  the  coarse  particle  mode  existing  beyond  1.0  /im.  The  data  show  the  accumulation  mode  is 
wefl  deftned  by  a  log-normal  distribution.  The  values  of  the  measured  volume  scattering  coefficient 
and  mean  particle  radius  of  the  accumulation  mode  increase  as  the  maximum  particle  volume  of  the 
mode  increases.  The  existence  of  an  accumulation  mode  was  almost  always  confined  to  the  mixing  layer. 
It  is  interesting  that  haze  layers  above  the  mixing  layer  were  found  to  have  a  distinct  coarse  particle  mode 
but,  generally,  no  distinct  accumulation  mode.  The  total  concentration  of  particles  in  each  of  the  two 
modes  appear  related  for  larger  values  of  concentration  in  the  coarse  particle  mode. 


1.  Introduction 

Sophisticated  numerical  techniques  have  been  de¬ 
veloped  by  a  number  of  investigators,  including  Dave 
(1978)  and  Takashima  (1975),  to  solve  for  the  radia¬ 
tion  emerging  from  an  atmosphere.  These  investi¬ 
gators  and  others,  however,  must  in  some  way  rely 
upon  a  knowledge  of  the  size  distribution  of  the 
atmospheric  aerosol  as  a  function  of  altitude,  season, 
and  location.  The  size  distributions  frequently  used 
are  the  haze  models  of  Deirmendjian  (1969),  who 
used  a  modified  gamma  function;  the  power-law 
size  distribution  function  of  Junge  (1963);  or  a  log¬ 
normal  distribution  as  described  by  Aitchison  and 
Brown  (1957).  Generally,  the  shape  of  the  distribu¬ 
tion  is  held  constant  while  only  the  particle  concen¬ 
tration  changes  with  altitude. 

There  are  few  measurements  to  support  the  use  of 
these  model  size  distributions.  A  series  of  airborne 
measurements  were  made  by  Blifford  and  Ringer 
(1969),  who  collected  aerosol  samples  in  the  central 
United  States  during  1966  at  altitudes  from  0.3  to  9 
km  above  the  ground.  A  majority  of  the  samples,  22 
of  them,  were  acquired  near  Scottsbluff,  Nebraska. 
Sampling  of  the  troposphere  has  also  been  con- 


1  Present  affiliation:  Air  Force  Office  of  Scientific  Research, 
Bolling  AFB,  Washington,  DC  20332. 


ducted  by  Laulainen  el  al.  (1978).  To  augment  these 
data  and  increase  our  understanding  of  the  atmos¬ 
pheric  aerosol  content,  an  extensive  measurement 
program  was  conducted  in  the  lower  troposphere  at 
seven  sites  in  northern  Europe  and  Great  Britain 
during  each  of  the  four  seasons.  This  report  presents 
the  results  of  analysis  of  some  of  these  data. 

2.  Method  of  measurement 

Aircraft  measurements  of  the  number  and  size  dis¬ 
tribution  of  the  atmospheric  aerosol  were  conducted 
along  tracks  of  constant  altitude  above  the  ground 
near  Bruz  (France),  Rodby  (Denmark),  Ahlhom 
and  Meppen  (W.  Germany)  during  the  spring  and  fall 
of  1976  and  the  summer  of  1977.  All  of  the  50  km  long 
sampling  tracks,  the  locations  of  which  are  shown  in 
Fig.  1,  were  over  land,  except  for  the  Rodby  track 
which  was  over  the  Baltic  Sea.  These  airborne  aerosol 
measurements,  described  by  Cress  (1980),*  were 
part  of  a  larger  atmospheric  optical  measurement 
program  (Duntley  et  al.,  1978)*  conducted  under 


*  Cress,  T.  S.,  1980:  Airborne  measurements  of  aerosol  size 
distributions  over  northern  Europe,  Vol.  1.  Rep.  AFGL-TR-80- 
0178,  Air  Force  Geophysics  Laboratory,  124  pp. 

1  Duntley,  S.  Q.,  R.  W.  Johnson  and  J.  I.  Gordon,  1978:  Air¬ 
borne  measurements  of  optical  atmospheric  properties,  summary 
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Fig.  I .  Location  of  the  ground-base  data  sites  (circles)  and  the  aircraft  tracks  (line  segments). 


the  sponsorship  of  the  Air  Force  Geophysics 
Laboratory. 

Sampling  was  performed  in  cloud-free  air  at  alti¬ 
tudes  from  — 150-6100  m  above  the  surface.  The 
altitudes  varied  from  one  flight  to  another  because 
of  changes  in  cloud  height  and  mixing  layer  depth. 
Whenever  possible,  two  samples  were  made  at  dif¬ 
ferent  altitudes  within  the  mixing  layer  and  two 
more  above.  The  mixing  layer  was  generally  capped 
by  an  inversion  layer.  The  sampling  time  of  the  Royco 
was  typically  4  min,  but  occasionally  samples  were 
taken  for  I  min  or  as  long  as  10  min,  depending 
on  the  situation  and  the  density  of  particles  being 
counted.  A  single  flight  required  from  30  min  to 
220  min  depending  on  the  number  of  altitudes  flown. 

The  instrumentation  on  the  C- 1 30  aircraft  included 
a  Royco  220  single  particle  counter  and  an  integrating 
nephelometer.  Sample  air  for  the  Royco  was  pro¬ 
vided  by  a  ram  air  isokinetic  sampling  probe  located 
on  top  of  the  C-130  fuselage  with  the  probe  entrance 
33  cm  above  the  aircraft  skin  and  45  cm  to  the  right  of 
the  radome  which  housed  the  nephelometer.  The 
entire  Royco  system  consisted  of  the  Royco  counter, 
a  Technical  Measurement  Corporation  Model  102 
Gamma  Scope  II  Pulse  Height  Analyzer,  and  a  Hew¬ 
lett-Packard  Clock/Printer.  The  electronics  of  the 
analyzer  were  adjusted  to  position  the  counting  of 
0.4  |xm  radius  size  particles  in  channel  18  of  the  total 


and  review  III.  Rep.  AFGL-TR-78-0286,  Visibility  Laboratory. 
University  of  California,  San  Diego,  Scripps  Institution  of  Ocean¬ 
ography,  La  Jolla,  108  pp.  (NTIS  AD  A073I2IJ. 


89  channel  range  covering  particle  sizes  0.2  to 
5.9  )tm. 

The  Royco  is  a  right-angle  scattering  counter  that 
was  calibrated  with  submicrometer  size  polystyrene 
spheres.  The  operating  characteristics  of  this  instru¬ 
ment  have  been  the  subject  of  a  number  of  in¬ 
vestigations  including  those  of  Quenzel  (1969),  Liu 
el  al.  (1974),  and  Cooke  and  Kerker  (1975).  The  two 
fundamental  sources  of  error  result  from  an  in¬ 
adequate  knowledge  of  the  index  of  refraction  of 
the  ambient  aerosol  and  modifications  to  the  particle 
size  distribution  which  may  occur  in  the  plumbing 
between  the  ambient  air  and  the  counter. 

Investigations  by  Cooke  and  Kerker  (1975)  and 
Quenzel  (1969)  have  demonstrated  that  significant 
errors  in  the  determination  of  particle  size  can 
occur  when  the  refractive  index  m  of  the  particles 
is  not  known.  Assuming  a  calibration  with  poly¬ 
styrene-latex  spheres  (m  =  1.58-0/),  Quenzel 
=  1.33-0/)  or  an  average  atmospheric  aerosol 
(m  =  1 .50—0/)  is  used,  the  sizing  error  is  less  than  a 
(m  =  1.95-0.66/)  over  the  range  of  particle  radius 
0.1 -3.0  //m.  If  the  refractive  index  of  water  (m 
=  1 .33— Of)  or  an  average  atmospheric  aerosol 
(m  =  1 .50— Of)  is  used,  the  sizing  error  is  less  than  a 
factor  of  2.0.  A  value  of  m  =  1.50-0.02/',  that 
includes  absorption,  produced  a  sizing  error  that  in¬ 
creased  with  increasing  particle  radius.  The  mag¬ 
nitude  of  the  error  was  difficult  to  determine  from 
Fig.  2  of  Quenzel  (1969)  but  its  maximum  is  esti¬ 
mated  to  be  in  the  range  2  to  3.  A  ground  based 
intercomparison  study  between  a  Royco  220  and 
active  and  classical  Particle  Measurement  Systems 
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Fig.  2.  The  aerosol  volume  per  increment  d  logr  as  a  function  of  particle 
radius  for  data  taken  at  Rodby  (Denmark)  on  25  October  1976  at  an  altitude  of  305 
m  above  the  sea  surface  (solid  line).  The  mode  radius  r  and  maximum  M,  of  the 
volume  distribution  accumulation  mode  are  shown.  A  log-normal  distribution 
with  r  =  0.45  gm  and  standard  deviation  n  -  0. 13  is  shown  (dashed  line)  for 
comparison. 


aerosol  counters  showed  agreement  to  within  a  factor 
of  2  to  5  in  particle  concentration  (Cress  and  Fenn, 
1978)/ 

The  physical  properties  of  the  ambient  aerosol 
sample  may  have  been  modified  in  its  passage  from 
outside  the  aircraft  to  the  Royco  scattering  chamber. 
Comparison  of  consecutive  I  min  data  samples, 
however,  indicate  that  the  repeatability  of  the  airborne 
Royco  was  good.  The  flow  rate  was  approximately 
47  cubic  centimeters  per  second.  Possible  modifica¬ 
tions  to  the  aerosol  sample  may  be  a  result  of  a 
drying  of  the  particles  or  a  selective  loss  of  particles. 
The  small  number  of  changes  in  the  direction  of 
flow  within  the  plumbing  probably  minimized  the 
effects  of  particle  loss.  This  is  supported  by  the  good 
agreement  between  this  data  set  and  the  impactor 
data  of  Blifford  and  Ringer  (1969)  as  shown  by 
Cress/  The  aerosol  measurements  most  likely  repre¬ 
sent  a  low  humidity  or  dry  particle  distribution 
because  some  evaporation  probably  occurred  in  the 
plumbing.  The  shape  of  the  distribution  was  probably 
preserved  due  to  the  small  size  range  measured. 
These  sources  of  error  are  discussed  in  greater 
detail  by  Cress/ 

The  integrating  nephelometer,  described  by  Duntley 
et  al..3  measures  the  volume  scattering  coefficient 
v  over  the  range  of  scattering  angle  5°  to  170°  and 
the  scattering  coefficient  at  30°  and  150°.  The 
nephelometer  operated  at  the  photopic  and  in  nar- 


*  Cress,  T.  S.  and  R.  W.  Fenn.  1978:  OPAQUE  aerosol  counter 
intercomparison  25  April  1977-4  May  1977.  Rep.  AFGL-TR- 
78004,  Air  Force  Geophysics  Laboratory,  56  pp. 


row  wavelength  bands  centered  at  0.478,  0.664  and 
0.765  /am.  Calibration  using  dry  nitrogen  and  com¬ 
parisons  with  an  Eltro  transmissometer  indicate  that 
the  accuracy  of  the  volume  scattering  coefficient 
is  within  20%  (Duntley  et  til..  1976)/ 

3.  Results  of  measurements 

Traditionally,  size  distribution  data  have  been 
presented  in  figures  of  \og(dN/d  logr)  versus  logr, 
where  r  is  the  particle  radius  and  N  is  the  total 
number  density  of  particles  smaller  than  r.  A 
presentation  of  this  type  has  been  found  to  mask  the 
multimode  nature  of  distributions  that  appears  when 
the  data  are  plotted  as  volume  (or  mass)  versus 
particle  size,  as  pointed  out  by  Whitby  et  al.  (1972) 
and  Shettle  (1975).  For  this  reason  it  was  decided 
to  use  the  volume  size  display  in  this  analysis. 

A  plot  of  dVId  logr  versus  logr  for  the  meas¬ 
urements  taken  at  Rodby  on  25  October,  1976  at 
an  altitude  of  305  m  above  the  sea  surface  is  shown 
in  Fig.  2.  Values  of  the  total  particulate  volume, 
dVId  logr,  over  an  increment  d  logr  are  shown  in  the 
figure  by  circles  for  a  few  of  the  89  data  channels. 
The  solid  curve  represents  a  fit  to  the  data  of  all  the 
channels.  The  inversion  height  was  about  700  m  in 
this  case  as  determined  from  temperature  profile 
measurements  made  concurrently  with  the  Royco 


5  Duntley,  S.  Q.,  R.  W.  Johnson  and  J.  I.  Gordon,  1976:  Air¬ 
borne  measurements  of  optical  atmospheric  properties  in  Northern 
Germany.  Rep.  AFGL-TR-76-0188.  Visibility  Laboratory,  Uni¬ 
versity  of  California,  San  Diego.  Scripps  Institution  of  Ocean¬ 
ography,  La  Jolla.  221  pp.  (NTIS  AD  A03I496). 
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Fig.  3.  The  mode  radius  r  as  a  function  of  the  maximum  volume  Mv  for  all  the 
data  collected  within  the  mixing  layer.  Curves  of  dVId  logr  for  constant  particle 
number  density  dNId  logr  (solid  lines)  show  the  r3  dependence  of  Mr.  Data  for 
which  the  relative  humidity  is  greater  than  or  equal  to  70%  are  indicated  by  the 
open  circles. 


measurements  and  from  estimates  of  the  surface 
haze  layer  thickness  made  by  an  on-board  observer. 
The  curve,  shown  in  Fig.  2  centered  about  the 
radius  0.45  fan,  is  closely  approximated  by  a  log¬ 
normal  curve  (dashed  line).  The  symmetry  and 
shape  of  the  data  curve  is  typical  of  the  analytical 
data.  Whitby  (1978)  termed  this  the  accumulation 
mode  and  stated  that  its  main  source  of  mass  is  from 
the  coagulation  of  the  smaller  size  aerosols  and  from 
gas  to  particle  conversion.  The  peak  in  the  curve  at 
about  3  fim  indicates  the  presence  of  the  coarse 
particle  mode,  which  according  to  Whitby  (1978) 
generally  results  from  mechanical  processes  such  as 
those  that  create  wind-blown  dust  or  sea  spray. 

a.  Accumulation  mode 

For  this  study  the  accumulation  mode  is  modeled 
by  the  value  of  the  mode  radius  f  the  standard 
deviation  a,  and  Mv — the  value  of  dVId  logr 
occurring  at  the  mode  radius.  Values  of  Mv  are 
plotted  against  those  of  f  in  Fig.  3  for  all  of  the 
seasons,  stations,  and  altitudes  below  the  inversion 
height  that  were  available  for  this  study.  Values  of 
Mv  >  40  /im1  cm'3  represent  data  collected  only  at 
Rodby  during  the  winter.  The  figure  shows  values 


of  M r  increasing  with  values  of  r.  Increasing  values 
of  A/,,  do  not  necessarily  imply,  however,  that  the 
particle  number  is  increasing.  Also  shown  in  the 
figure  are  curves  of  dVId  logr  for  constant  particle 
number,  dNId  logr,  which  illustrates  the  r3  depend¬ 
ence  inherent  in  the  volume  size  distribution.  For 
values  of  r  <  0.34  /xm  there  is  an  increase  in  the 
measured  particle  number  with  increasing  r.  For 
values  off  >  0.34  ^im,  however,  the  particle  number 
appears  to  be  constant  with  increasing  values  of 
r.  The  data  presented  in  Fig.  3  suggest  that  there  is  a 
functional  relationship  between  f  and  Mc.  This  rela¬ 
tionship  will  be  more  fully  explored  for  different  air 
mass  types  using  the  complete  set  of  Royco  meas¬ 
urements.  Plots  comparing  a,  which  range  in  value 
from  0.08  to  0.20,  with  values  of  both  Af„  and  tot&l 
particle  volume  V  of  the  mode  show  no  discernible 
relationship  between  these  parameters.  This  is  prob¬ 
ably  due  to  the  sensitivity  of  the  determination  of 
(7  to  the  fitting  procedure.  In  order  to  more  easily 
visualize  the  character  of  the  accumulation  mode 
and  to  simplify  its  log-normal  description,  the  value 
Mt  =  VM2ir)1,s  is  used  for  the  analysis  of  this  mode. 

There  does  not  appear  to  be  a  dear  dependence 
on  relative  humidity  as  measured  with  a  Cambridge 
Model  137-C3  aircraft  hygrometer  system.  The  open 
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circles,  representing  cases  where  the  relative  humidity 
was  greater  than  or  equal  to  70%,  predominate  at 
the  larger  radii  as  shown  in  Fig.  3.  This  apparent 
dependence  on  relative  humidity,  however,  is  mainly 
due  to  a  small  data  sample  taken  at  Rodby  in  an 
unusually  heavy  haze.  If  this  Rodby  data  is  neglected 
from  the  figure,  then  any  tendency  for  the  rela¬ 
tive  humidity  to  be  *70%  when  the  value  of  the 
mode  radius  is  *0.34  ^m  is  probably  not  significant. 
There  is  a  probability  level  of  ~0.05  that  the  value 
of  chi  square  will  be  equalled  or  exceeded  by  a 
random  sample.  If  values  greater  than  or  equal  to 
80%  are  considered  there  again  appears  to  be  no 
clear  dependence  on  relative  humidity.  This  lack 
of  a  clear  dependence  on  relative  humidity  may 
be  a  result  of  drying  of  the  aerosol  in  the  intake 
plumbing. 

In  addition  to  measurements  of  the  aerosol  size 
distribution,  the  volume  scattering  coefficient  s  at 
X  =  0.664  fim  was  also  measured.  The  scattering  co¬ 
efficient  may  be  expressed  as 

rlogr, 

s  =  7T  Q.(m,x)r3  — - d  log/-, 

Jtogri  d  lOgr 

where  x  -  2nr/A  is  the  size  parameter  of  a  spherical 
particle,  Q,(m,  x)  is  the  efficiency  factor  for  scatter¬ 
ing  (Van  de  Hulst,  1957),  and  r,  and  r2  are  the  lower 
and  upper  limits  of  the  particle  radii  under  con¬ 
sideration.  The  relationship  between  s  and  the  model 
parameters  can  be  shown  by  using  a  log-normal 
description  of  the  aerosol  volume  distribution  and 
the  expression  dVId  log/-  =  %  nr3  (dNId  logr).  After 
simple  manipulation,  s  can  be  written  in  the  form 

3  rloirt 

J  =  -  Q,(m,x)r'Mv 

4  Jlocri 

x  exp{-0.5((logr  -  logr)/<rjJ}d  logr.  (i) 

Even  though  the  Royco  is  sensitive  to  a  narrower 
particle  size  range  than  the  nephelometer,  cal¬ 
culations  show  that  the  distribution  below  0.4  /um, 
typically,  adds  <15%  to  the  total  scattering  at  X 
=  0.455  //.m.  At  the  wavelength  0.664  fim  this 
percent  contribution  is  even  smaller,  so  the  results  of 
measurements  by  the  Royco  and  the  nephelometer 
should,  as  seen  in  Eq.  (1),  be  strongly  related. 

Values  of  the  volume  scattering  coefficient  s  for 
the  wavelength  0.664  fim  are  plotted  in  Fig.  4  as  a 
function  of  Mv.  Though  values  of  j  are  not  uniquely 
dependent  upon  values  of  Af„,  as  seen  in  Eq.  (1), 
the  figure  does  show  a  trend  of  increasing  s  with 
increasing  Mc.  This  trend  can  be  explained  by  a  shift 
in  the  volume  distribution  to  larger  radii  for  constant 
particle  number,  which  increases  M,.,  and  by  increases 
in  particle  number  which  occur  at  the  smaller  values 
of  Mc  as  shown  in  Fig.  3.  The  scatter  in  the  data,  in 
part,  is  a  result  of  the  fact  that  s  was  not  always 
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Fic.  4.  The  volume  scattering  coefficient  s  for  the  wavelength 
0.664  fim  as  a  function  of  the  maximum  volume  Mr.  Cases  were 
s  and  Mc  were  measured  simultaneously  are  circled. 

measured  simultaneously  with  M„.  At  times  there 
was  a  two  to  three  hour  difference  between  the  meas¬ 
urements.  There  is  an  apparent  decrease  in  the  scatter 
of  the  data  when  simultaneous  measurements  of  s 
and  Mv  are  plotted  (circled  data  points).  Some 
residual  scatter  can  be  expected  due  to  the  effects 
of  variations  in  relative  humidity  which  were 
detected  by  the  nephelometer  but  probably  not  by 
the  Royco  due  to  drying  in  the  intake  plumbing. 

In  five  of  the  flights  studied,  the  air  was  so  clear 
and  the  particle  counts  so  low  that  there  was  no 
distinct  accumulation  mode  measured  at  altitudes 
below  the  inversion  height.  All  of  these  cases  oc¬ 
curred  during  the  winter  season  at  Bruz,  Rodby  and 
Meppen  and,  in  all  but  one  case,  the  air  mass  type 
was  maritime.  The  average  meteorological  range, 
as  estimated  by  the  on-board  observer,  was  always 
>14  km.  The  values  of  dVId  logr  at  0.3  /um  were 
rarely  >2  fim3  cm*3  (corresponding  to  a  measured 
V  over  the  range  of  the  mode  of  <  1  jim3  cm-3)  and 
s  measured  for  the  0.664  ^tm  band  was  <6.8  x  10~3 
m*'.  These  values  of  dV/d  log r  and  s  are  in  agree¬ 
ment  with  the  general  trend  of  the  data  in  Fig.  4, 
even  though  no  value  of  Mr  could  be  determined 
for  these  five  cases. 

To  investigate  the  variation  of  M,  with  respect  to 
altitude,  values  of  M,.  were  plotted  against  the  ratio 
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Fig.  5.  The  maximum  volume  in  the  accumulation  mode  Mr  measured  at  Rodby 
on  25  October  1976  (circles)  and  on  26  October  1976  (crosses)  as  a  function  of 
altitude  ratio,  the  ratio  of  the  aircraft  altitude  during  measurement  to  the  thick¬ 
ness  of  the  mixing  layer. 

of  the  measurement  altitude  to  the  altitude  of  the  aerosol  measurements.  From  the  altitude  0.4  to  1.0 
inversion  base.  Data  taken  at  Rodby  on  25  and  26  the  temperature  profile  for  both  days  was  close  to  the 
October  are  presented  in  Fig.  5.  In  the  mixing  dry  adiabatic  lapse  rate.  In  addition,  below  the 
layer — an  altitude  ratio  of  <1.0 — the  values  of  Mv  altitude  ratio  0.4  the  temperatures  for  4  July  were 
decrease  with  increasing  altitude  for  both  days.  For  larger  than  those  for  7  July.  This  is  probably  because 
ratios  >  1.0  (note  the  change  of  scale)  the  values  the  7  July  data  was  taken  in  the  morning  and  the 
of  Mv  are  near  zero.  4  July  data  was  collected  in  the  afternoon.  The 

All  of  the  Rights  investigated  show  small  values  scatter  in  the  7  July  data  at  the  low  altitude  ratios 
of  Mr  at  altitude  ratios  of  >1.0.  The  decrease  in  is  in  part  a  result  of  the  time  required  to  take  these 
Mr  with  increasing  altitude  in  the  mixing  layer,  data,  ~3. 5  h.  In  contrast,  three  measurements  made 
however,  is  not  always  present.  The  results  of  at  Bruz  on  6  July  over  a  4  min  period  at  an  altitude 
measurements  made  at  Bruz  on  4  and  7  July  are  ratio  of  0.3  show  Mv  ranging  from  only  15  -18 
presented  in  Fig.  6.  There  is  a  decrease  in  the  values  /urn3  cm-3, 
of  Mv  with  increasing  altitude  in  the  mixing  layer 

for  the  7  July  data  but  not  for  the  4  July  data,  which  h.  Coarse  particle  mode 
appears  invariant  with  altitude.  The  temperature 

profile  to  an  altitude  ratio  of  0.4  shows  7  July  as  more  Very  often,  low  particle  counts  were  measured  in 
thermally  stable  than  4  July  at  the  time  of  the  the  size  range  of  the  coarse  particle  mode.  These  low 
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Fig.  6.  The  maximum  volume  in  the  accumulation  mode  M,,  as  a  function  of  alti¬ 
tude  ratio  for  data  collected  at  Bruz  on  7  July  1977  (circles)  and  on  4  July  1977 
(crosses). 


counts  made  it  difficult  to  determine  values  of  the 
model  parameters  M„,  a  and  r  for  this  mode. 
Consequently,  the  coarse  particle  mode  was 
represented  by  the  integrated  volume  concentration 
V  of  particles  from  particle  radius  2.0  to  5.7  /xm.  A 
plot  of  V  versus  the  altitude  ratio  for  Bruz  on  4  and  7 
July  is  shown  in  Fig.  7.  Within  the  mixing  layer,  the 
values  of  V  decrease  with  increasing  altitude  for  both 
days.  Above  the  mixing  layer,  the  increase  in  the 
values  of  volume  concentration  with  increasing 
altitude  ratio  indicate  the  presence  of  a  second  haze 
layer.  The  on-board  observer  confirmed  the  exist¬ 
ence  of  this  second  haze  layer.  It  is  interesting  that 
there  is  no  indication  of  a  second  haze  layer  in  the 
accumulation  mode  data  of  Fig.  6.^  When  there  was 
no  second  haze  layer  above  the  mixing  layer,  the 
trend  of  V  and  A/„  with  altitude  was  similar. 


The  relationship  between  Mr  of  the  accumulation 
mode  and  V  of  the  coarse  particle  mode  is  repre¬ 
sented  in  Fig.  8,  which  includes  all  the  analyzed  data 
from  the  mixing  layer.  The  total  volume  concentra¬ 
tion  based  on  a  log-normal  fit  to  the  accumulation 
mode  can  be  determined  by  using  the  expression 
V  =  Mv&(2n)'12  where  the  average  &  for  the  data  set 
is  0.11.  Values  of  volume  concentration  greater 
than  3.5  jxm3  cm-3  represent  data  from  Rodby  (25 
and  26  October)  taken  in  a  saturated  layer  caused 
by  the  warm  water  of  the  Baltic  and  strong  winds. 
The  connected  values  in  Fig.  8  show  the  variations 
that  can  occur  over  2-3  h  at  a  given  site  and  alti¬ 
tude.  The  set  of  data  points  occurring  at  Mr  >  15 
/xm3  cm~3  and  V  <  1 .0 /xm3  cm-3  are  from  two  flights 
conducted  on  4  August  at  the  Meppen  site.  During 
the  measurements,  the  air  mass  was  classified  as 
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Fig.  7.  The  volume  concentration  of  particles  in  the  coarse  particle  mode  V  as  a 
function  of  altitude  ratio  for  data  collected  at  Bruz  on  7  July  1977  (circles)  and  on 
4  July  1977 (crosses). 


stagnant  with  ground  stations  in  the  area  reporting 
haze,  and  winds  <  4  m  s  '.  The  data  in  Fig.  8  indi¬ 
cate  that  values  of  Mv  generally  increase  with  values 
of  V,  except  for  the  4  August  data.  The  increase, 
however,  is  strongly  dependent  on  only  two  days 
of  data  from  Rodby.  It  probably  represents  the 
effect  of  different  air  mass  types  on  values  of  V  and 
Mr  and  not  the  effect  of  a  correlation  between  the 
two  parameters.  According  to  Whitby  (1978),  the 
two  modes  have  different  source  and  production 
mechanisms,  so  no  strong  correlation  need  exist  be¬ 
tween  the  two  modes. 

4.  Discussion  and  conclusions 

Aircraft  measurements  of  particle  size  distribu¬ 
tions  were  made  at  several  altitudes  within  and 
above  the  mixing  layer  at  sites  near  Ahlhom  and 
Meppen  (W.  Germany),  Rodby  (Denmark)  and  Bruz 
(France).  The  experimental  information  gathered  at 


the  sites  was  analyzed  in  terms  of  volume-size 
distributions  instead  of  the  standard  particle  number 
distributions.  The  data  support  the  following  conclu¬ 
sions  for  the  aerosol  characteristics  at  the  measure¬ 
ment  sites: 

1)  There  are  at  least  two  distinct  modes  com¬ 
prising  the  size  distribution.  The  accumulation  mode 
is  well  defined  by  a  log-normal  distribution  char¬ 
acterized  by  <r,  r  and  Mv.  Values  of  Mv  appear  to 
increase  with  values  of  r.  This  is  caused  by  an 
increase  in  particle  number  with  increasing  F  for  F 
less  than  0.34  jtm,  and  it  is  a  result  of  the  f*  depend¬ 
ence  of  Mv.  No  clear  relationship  could  be  found 
between  <r  and  the  other  model  variables.  An  ap¬ 
parent  relationship  between  Mv  and  the  volume  con¬ 
centration  of  particles  in  the  coarse  particle  mode 
V  is  probably  the  result  of  the  different  air  mass 
types  present  during  the  study. 

2)  The  values  of  Mv  measured  in  the  mixing  layer 
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Fig.  8.  The  volume  concentration  of  particles  in  the  coarse  particle  mode  V 
as  a  function  of  maximum  volume  in  the  accumulation  mode  M,  for  data  collected 
within  the  mixing  layer.  The  connected  points  show  typical  variations  occurring 
over  a  period  of  2-3  h  at  a  particular  site  and  altitude. 


appear  well  correlated  with  the  measured  values  of 
the  volume  scattering  coefficient. 

3)  The  existence  of  an  accumulation  mode  was 
confined  to  the  mixing  layer  for  the  majority  of  the 
cases  studied.  During  the  winter  season  for  condi¬ 
tions  of  good  visibility  in  a  maritime  airmass,  how¬ 
ever,  the  particle  counts  were  generally  so  low  in  the 
mixing  layer  that  a  distinct  accumulation  mode  was 
not  evident. 

4)  Within  the  mixing  layer  the  value  of  Mr  is  not 
always  constant  with  altitude,  but  may  depend  on 
the  stability  of  the  layer  and  the  strength  of  the 
capping  inversion. 

5)  Size  distribution  models  for  a  specific  air  mass 
type  which  change  only  the  total  number  of  particles 
may  not  always  be  correct.  Significant  changes  in 
f  with  changing  Mr  for  a  given  air  mass  were  ob¬ 
served.  Work  is  currently  underway  to  improve  the 
determination  of  air  mass  type  and  to  better  under¬ 
stand  its  influence  on  the  properties  of  the  two  modes. 

6)  The  coarse  particle  mode  was  not  confined  to 
the  mixing  layer.  It  is  interesting  that  haze  layers 
above  the  mixing  layer  were  found  to  have  a  distinct 
coarse  particle  mode  but,  generally,  no  accumula¬ 


tion  mode.  Its  difficult  to  determine  the  origin  of 
these  layers  which  may  be  a  specific  source,  or 
they  may  simply  be  aged  aerosols. 

7)  The  aerosol  measurements  probably  represent 
a  low  humidity  or  dry  particle  distribution  due  to 
evaporation  in  the  intake  plumbing.  Its  likely  that 
the  shape  of  the  distribution  of  each  mode  was  pre¬ 
served  because  of  the  small  range  in  radius  spanned 
by  each  of  the  two  modes.  Each  of  the  modes  may 
have  been  affected  differently,  however,  due  to  dif¬ 
ferences  in  particle  size  and  possible  differences 
in  chemical  composition. 

The  rest  of  the  aerosol  data  collected  in  northern 
Europe  are  being  studied  at  present.  With  the 
added  data  set,  which  is  approximately  twice  as  large 
as  was  used  in  this  investigation,  attention  will  be 
focused  on  how  the  characteristics  of  the  accumula¬ 
tion  and  coarse  particle  modes  change  with  air  mass 
type  and  season.  Techniques  to  better  determine 
o-  for  the  accumulation  mode  and  M,  ,  r  and  a-  for 
the  coarse  particle  mode  are  being  developed.  In 
addition,  the  modeling  of  these  modes  at  different 
altitudes  in  the  lower  troposphere  based  on  a  knowl- 
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edge  of  the  inversion  height,  volume  scattering 
coefficient  and  air  mass  type  is  being  pursued. 
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